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Hadronic r decays provide a clean laboratory for the precise study of quantum chromodynamics 
(QCD). Observables based on the spectral functions of hadronic r decays can be related to QCD 
quark-level calculations to determine fundamental quantities like the strong coupling constant, 
parameters of the chiral Lagrangian, |Vu S |, the mass of the strange quark, and to simultaneously 
test the concept of quark-hadron duality. Using the best available measurements and a revisited 
analysis of the theoretical framework, the value as(mj) = 0.345 ± 0.004 cxp ±0.009 t h is obtained. 
Taken together with the determination of as(M^) from the global electroweak fit, this result leads 
to the most accurate test of asymptotic freedom: the value of the logarithmic slope of aj 1 (s) is 
found to agree with QCD at a precision of 4%. The r spectral functions can also be used to 
determine hadronic quantities that, due to the nonperturbative nature of long-distance QCD, 
cannot be computed from first principles. An example for this is the contribution from hadronic 
vacuum polarization to loop-dominated processes like the anomalous magnetic moment of the 
muon. This article reviews the measurements of nonstrange and strange r spectral functions and 
their phenomenological applications. 
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I. INTRODUCTION 

Because of its relatively large mass and the simplicity of its decay mechanism, the r lepton offers many 
interesting, and sometimes unique, possibilities for testing and improving the Standard Model (SM). These 
studies involve the leptonic and hadronic sectors and encompass a large range of topics, from the measurement 
of the leptonic couplings in the weak charged current and the search for lepton flavor violation, providing 
precise lepton universality tests, to a complete investigation of hadronic production from the Quantum Chro- 
modynamics (QCD) vacuum. For the latter case, the r decay results have proven to be complementary to 
those from e + e~ data, allowing us to perform detailed studies at the fundamental level through the deter- 
mination of the spectral functions, which embody both the rich hadronic structure seen at low energy, and 
the quark behavior relevant in the higher energy regime. The spectral functions play an important role in 
the understanding of hadronic dynamics in the intermediate energy range. They represent the basic input for 
QCD studies and for evaluating low-energy contributions from hadronic vacuum polarization. This is required 
for the calculation of the muon anomalous magnetic moment and the running of the electromagnetic cou- 
pling constant. The topics of interest also include: studying isospin violation between the weak charged and 
electromagnetic hadronic currents, evaluation of chiral sum rules making use of the separate determination of 
vector and axial-vector components, and a global QCD analysis including perturbative and nonperturbative 
contributions, offering the possibility for a precise extraction of the strong coupling at a relatively low energy 
scale, and hence providing a sensitive test of the running when compared to the value obtained at the Mz 
scale. The spectral functions into strange (|A5| = 1) final states, though Cabibbo-suppressed, can be used to 
determine the s-quark mass and the CKM element \V US \. In addition, the study of t + t~ pair production at 
LEP has proven to be an invaluable tool for investigating the neutral current sector of the electroweak theory, 
particularly through the measurement of the r polarization. 

This article reviews the measurements of nonstrange and strange r spectral functions and their phenomeno- 
logical applications. It is organized as follows. We begin in Section II with an overview on the experimental 
conditions leading to the measurement of the t branching fractions (Section III) and hadronic spectral func- 
tions (Section IV) by the ALEPH and OPAL experiments at LEP as well as CLEO at CESR. Using isospin 
symmetry, the r vector spectral functions are compared to the corresponding e+e~ annihilation data in Sec- 
tion V. We discuss in this context isospin-breaking and radiative corrections. The use of spectral functions 
in dispersion relations to compute hadronic vacuum polarization is reviewed in Section VI. A comprehensive 
review of hadronic r decays within QCD is presented in Sections VII and VIII. The main topics here are the 
determinations of the strong coupling constant and the mass of the strange quark at the scale of the r mass 
together with \V US \. We discuss nonperturbative contributions to the r hadronic width and present results 
from chiral QCD sum rules. 

For all these topics we critically examine experimental and theoretical limitations, and provide detailed 
comparisons between experiments as well as perspectives for the future. 

II. THE EXPERIMENTAL STAGE 

A. Different experimental conditions 

Tau physics in the nineties has been dominated by two rather different and in many views complementary 
experimental facilities: on one hand the LEP experiments ALEPH, DELPHI, L3 and OPAL, operating at the 
Z resonance at center-or-mass (CM) energies of 91.2 GeV, and on the other hand CLEO at CESR running at 
the T(4S) resonance (10.6 GeV). The r pair production and decay characteristics at LEP were: 

• a mean r flight length of 2.2 mm. 

• a strong boost of the produced r's leading to a collimated back-to-back event topology; hence tracks and 
photon clusters overlap in the detector, requiring excellent double hit resolution in the tracking system 
and a highly granular electromagnetic calorimeter. 

• a large track reconstruction efficiency and low nuclear interactions and multiple scattering due to the 
high energetic particles in the final state. 

• since hadronic background from Z decays occurs with significantly larger particle multiplicity, the selec- 
tion of t pairs is almost background free (< 1 — 2%), and the efficiency is large (up to 96% at ALEPH, 
depending on the final state). 
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• limited statistics: up to 1.6 x 10 5 r pairs have been recorded by each of the LEP experiments between 
1990 and 1995. 

For CLEO, one has 

• a mean t flight length of 0.25 mm. 

• well separated r decay products in the detector. 

• significant hadronic background from e + e~ — > qq continuum events due to the relatively low average 
multiplicity in the final state of these events. This requires r tagging, i.e., one r is reconstructed 
("tagged") in a leptonic decay or a hadronic mode with less than three hadrons in the final state. The 
recoil side of the tagged event is then used as an unbiased sample for physics studies. 

• large statistics: for example, the t~ — > ir~ir°v T spectral function analysis (CLEO Coll., 2000d) uses a 
data sample that contains 3.2 x 10 6 produced r pairs. 

Due to its superior statistics, CLEO discovered a number of rare modes, and performed studies of their 
hadronic substructures. Examples for these are t~ — ► r}K~v T (CLEO Coll., 1996) and t~ — ► r]ir~ir°v T (CLEO 
Coll., 1992), the latter one proceeding through a vector current that is due to the Wess-Zumino chiral 
anomaly (Decker and Mirkes, 1993; Pich, 1987). The anomaly violates the rule that the weak vector and 
axial- vector currents produce an even and odd number of pseudoscalars, respectively. The t~ — > rj-K^u T mode 
is a G-parity suppressed second-class current, which has not been observed yet. 

For the r decays with a branching fraction greater than 0.1%, the cleaner LEP events provided the more 
precise branching fraction measurements, many of them with a sub-percent relative accuracy, which are a 
crucial ingredient to the QCD study of hadronic r decays. The Cabibbo-suppressed r decays with |A5| = 1 
have been studied by all experiments, albeit with limited success due to the restricted particle identification 
capabilities (except for DELPHI and CLEO III). Much progress is expected in this area from the B-Factory 
experiments BABAR and Belle, which have dedicated Cherenkov particle-identification devices. 

B. Investigating the resonance structure in r decays 

Most of this review will be concerned with the interplay between measurements involving hadronic (or 
more precisely semileptonic) r decays and QCD. This may appear a bit paradoxical in view of the relatively 
low energy scale of r physics. Indeed, r decays to hadrons (+ v T ) reveal a rich structure of resonances, the 
study of which opens an interesting field of research on meson dynamics. The leptonic environment provides 
a significant advantage for isolating clean hadronic systems and measuring their parameters. Many r decay 
modes have been investigated for their resonance structure. The phenomenologically important 7r~7r° vector 
state is analyzed in detail in Section V. 

Taking advantage of its large statistics sample, the CLEO Collaboration was able to perform partial-wave 
analyses of multi-hadron final states. In particular, the (37r)~7r° mode is expected to proceed through the 
vector current dominated by p, p(1450), p(1700) isovector mesons. Resonances are also produced in two- 
body decays, such as ivn, rjir, ain. CLEO measured (CLEO Coll., 2000c) the ivn contribution, which they 
successfully modeled with interfering p, p(1450), p(1700) resonances as shown in Fig. 1. At least the first two 
contributions are needed to reproduce the experimental mass spectrum. The helicity angle distribution in the 
lutt~ center-of-mass is found to be consistent with the expected / = 1 wave (Fig. 1), whereas a wrong G-parity 
contribution arising from a weak second-class current, such as t — > b^v T , b± — > luit~ , would lead to an even 
I value. Such a contribution has been bound to be less than 5.4% of the ojit~ rate. Besides ujtt~ , the total 
(37r) _ 7r° mode is shown to be dominated by pirir, consistent with originating from a\-K, a\ — > pir. 

The decay r — ► v T (Zit)~ is the cleanest mode to study axial- vector resonance structure. The spectrum is 
dominated by the J p = 1 + a\ state, known to decay essentially through pir. The hadronic structure of this 
decay has been studied in a model-independent analysis by the OPAL collaboration, where no evidence for 
vector or scalar currents has been found (OPAL Coll., 1995a, 1997). A comprehensive analysis of the 7r~27r° 
channel has been presented by CLEO. First, a model-independent determination of the hadronic structure 
functions gave no evidence for non-axial- vector contributions (< 17% at 90% CL) (CLEO Coll., 2000b). 
Second, a partial- wave amplitude analysis was performed (CLEO Coll., 2000a): while the dominant pir mode 
has been confirmed, it came as a surprise that an important contribution (<~ 20%) from scalars (<r, /o(1470), 
/2(1270)) was found in the 2tt subsystem. The precisely determined — ► 7r~27r° line shape shows the opening 
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FIG. 1 Left: the spectral function (see Section IV for the definition) for the decay r~ — > u>n~v T vs. the hadronic 
mass from the CLEO analysis (CLEO Coll., 2000c). The curves are the results of fits to various combinations of p, 
p(1450), p(1700) resonance contributions. Right: the angular distribution in the um center-of-mass frame compared to 
predictions for different (I) partial waves. Second-class currents would reveal themselves as I — 0, 2 waves. 



of the K*K decay channel. Some interference with a higher mass state (a^) was advocated in (DELPHI Coll., 
1998). However, no conclusive evidence was found in the CLEO analysis, with a sensitivity about tenfold 
higher. 

Evidence for the decay r — > fi(1285)n~ i/ T has been obtained, with /i decay into i]ir + ir~ and ?77r 7r (CLEO 
Coll., 1997a), and recently into 2tt + 2it~ (BABAR Coll., 2005b), the latter decay mode being consistent with 
/9°7r + 7T~. Thus many spectroscopy issues can be tackled with r leptons in a clean environment for single 
resonance studies. 



III. TAU BRANCHING FRACTIONS 

Tau hadronic spectral functions (Section IV) are obtained from the measured vector and axial-vector in- 
variant mass spectra, normalized to the corresponding inclusive branching fractions. These are derived from 
the sum of exclusive branching fractions. Before discussing the spectral functions, we shall thus briefly review 
the measurement of the r branching fractions. 



A. The measurement of branching fractions 

In principle, the measurement of the branching fractions of exclusive r decays relies on the identification of 
the corresponding final state and the knowledge of the number of the produced r leptons in the sample under 
study. The latter number can be derived from the production cross section, known to high precision in e + e~ 
annihilation, and the measured luminosity. The measurement of the final state however requires excellence in 
many aspects of detector technology and analysis techniques: charged particle tracking and identification for 
e/ii/w/K separation, photon detection and ir° (77) reconstruction. 

Since very clean and high-efficiency t + t~ samples have been produced at LEP, another approach can 
be considered where all decay channels are simultaneously measured. This global method was pioneered by 
CELLO at PETRA (CELLO Coll., 1990) with much smaller statistics and less favorable conditions than those 
encountered at LEP. The detector performance achieved by ALEPH, in particular the excellent granularity 
of its electromagnetic calorimeter, made it possible to implement this method and take full advantage of the 
LEP conditions on the Z peak (ALEPH Coll., 1992, 1996a,b, 2005). 

The ALEPH global analysis of r branching fractions proceeds from the selection of the r + r~ sample, 
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based on the rejection of non-r backgrounds: Z decays into e + e~, and qq pairs, 77-induced processes 

e + e~ — > e + e~ X (where X can be any fermion-antifermion pair), and cosmic-ray events. The kinematic 
properties of these backgrounds and their characteristic detector response permit to achieve a high t+t~ 
selection efficiency of 78.9% (91.7% when the t + t~ direction points into the active area of the detector) and 
a low background fraction of 1.2%. Each r decay is classified on the basis of charged particle and photon/-^ 
multiplicities, and the nature of the charged particles. Charged particle identification is achieved through 
the response of electromagnetic and hadronic calorimeters, and the ionization loss measurement. Photons 
need to be identified among the neutral electromagnetic clusters, some of which are induced by charged 
hadron interactions or by spatial fluctuations of nearby showers (fake photons). Fourteen classes are thus 
defined in (ALEPH Coll., 2005) for reconstructed decays, where the last one collects the fraction (3.6%) of 
rejected candidates corresponding to the cases where particle identification is not reliable. The first 13 classes 
correspond to the following final states produced with v T : e~~P e , pTV^ h~ ', h~ir°, h~2ir°, h'Sir , ft,~47r°, 
2h~h+, 2h-h + Tr°, 2/j-/i+27r°, 2h-h+3TT°, 3h~2h+, 3/i"2/i+7r°, where h stands for a pion or kaon. 
The branching fractions are determined using 

nf»-n** = 5>>Af° d , (1) 
3 

^yprod 

Bi = " a > ( 2 ) 

Z-'j 3 

where n° hs (n^ ks ) is the number of observed t candidates (the non-r background events) that are reconstructed 
in class i, the cross-talk efficiency of events that are produced in class j and reconstructed in class i, and 
_/V? ro the number of produced events in class j. The efficiency matrix sji is determined from Monte Carlo 
simulation, based on the event generator KORALZ07 (Jadach et at, 1994) and including the full detector 
response. The simulation has been calibrated using data for a large number of observables. Among these are 
particle identification efficiencies, t + t~ selection efficiency and the number of fake photon candidates. The 
selection efficiency remains constant within ± <~ 5% for all 13 topologies considered. The linear system of 
equations (1) is solved for the 7Vj 3rod unknown, hence the branching fractions Bj. 

Extensive systematic studies are performed to assess possible biases in the decay assignment. These studies 
are conducted on data samples in order not to rely on the detector simulation, however detailed it can be. 
They include: event selection, non-r background, particle identification, photon detection efficiency, converted 
photons, photon identification, fake photons, ir° reconstruction, bremsstrahlung and radiative photons, 7r° 
Dalitz decays, secondary nuclear interactions, tracking of charged particles, dynamics in the Monte Carlo 
generator. The dominant systematic uncertainties for the hadronic modes stem from the photon/71" treatment 
and secondary interactions. 

While measurements of branching fractions larger than 0.1% are more precise at LEP because of its small 
systematic uncertainties, the larger statistics accumulated by CLEO makes it more competitive for smaller 
branching fractions. The different experimental conditions also provide complementary information and valu- 
able cross checks. 

The ALEPH global analysis assumes a standard t decay description. However, unknown decay modes not 
included in the simulation may exist. Since large detection efficiencies are achieved in the t~t + selection, 
which is therefore robust, these decays would be difficult to pass unnoticed. An independent measurement 
of the branching fraction for undetected decay modes, i.e., modes not passing the selection cuts, using a 
direct search with a one-sided t tag, was done by ALEPH (Snow, 1993), limiting this branching fraction 
to less than 0.11% at 95% CL. This result justifies the assumption that the sum of the branching fractions 
for visible r decays is equal to one. The consistency between the results and the r decay description in the 
simulation can be investigated further by verifying that the "branching fraction" for a 14 th class in the decay 
classification, which collects rejected events, turns out to be zero within the measurement errors. The result, 
$14 = (0.066 ± (0.042)%, provides a nontrivial test of the overall procedure at the 0.1% level for branching 
fractions. 

The global method at LEP reaches its limit for decay modes with small branching fractions, which compete 
with leading channels with fractions that could be more than two orders of magnitude larger. This limit is 
also determined by the overall r statistics. It turns out that the two factors contribute at about the same 
level. Thus final states with more than five hadrons are poorly determined at LEP, whereas CLEO can use 
its larger statistics to apply stricter cuts for these higher-multiplicity decays with low branching ratios. 
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A final check of the global analysis can be performed by comparing the sums of the relevant exclusive 
decay fractions to the topological branching fractions Bi, where i is the charged particle multiplicity in 
the decay. Even though the latter branching fractions have only limited direct physics implications, their 
determination represents a valuable cross check as the measurement does not depend on photon reconstruction 
in a calorimeter. Assuming a negligible contribution from charged-particle multiplicities larger than five, in 
agreement with the 90% CL limit achieved by CLEO (B 7 < 2.4 x 1(T 6 (CLEO Coll., 1997c)), ALEPH finds: 

B 3 = (14.652 ± 0.067 ± 0.086)% , (3) 
B 5 = (0.093 ± 0.009 ± 0.012)% , (4) 

in agreement with the values from a dedicated topological analysis performed by DELPHI (DELPHI Coll., 
2001), 

#3° p = (14.569 ± 0.093 ± 0.048)% , (5) 
2?5° p = (0.115 ± 0.013 ± 0.006)% . (6) 



B. Leptonic branching fractions 

The measurement of the two leptonic decay modes of the r is of great importance in our context, since they 
are complementary to the hadronic modes, the fraction of which can be determined by a simple difference. 
The leptonic branching ratios also turn out to be the most precisely determined r decay fractions, not affected 
by systematic uncertainties that are specific to the presence of hadrons in the detector. 

Before examining the experimental situation, we shall recall the theoretical expectations within the standard 
V — A theory with leptonic coupling gi at the WtP( vertex. The L — fi, t leptonic [l = e, \i) partial widths 
can be computed, including radiative corrections (Marciano and Sirlin, 1988) and safely neglecting neutrino 
masses: 

G L G e m 5 L J m 2 e \ L cL 

, nn .s / I IT I °W°-y 

where 



r ( ^^ (7 )) = ^^/^j^, (7) 



4V2M 



2 

W 

rL , 3 m l 

5w = 1 + - 1 



W 

s l = 1 + o^)/25 2 
7 2?r V 4 

f(x) = 1 - 8x + 8x 3 - x 4 - I2x 2 \nx (8) 

Numerically, the radiative and W propagator corrections are small: 5^ = 1 + 2.9 x 10~ 4 , 5l f = 1 — 43.2 x 10~ 4 , 
= 1 + 1.0 x 10- 6 , and 5Hj = l- 42.4 x lO" 4 . 
The experimental situation on the electronic and muonic branching fractions, B e and B^, is given in Fig. 2. 
The results for their average are 

B e = (17.821 ±0.052)% , (9) 
Bf, = (17.332 ±0.049)% , (10) 

with a relative precision of 0.3%. The values for B e and B^ agree with the universality assumption of equal 
leptonic couplings, g e — — g T . Their ratio 

= 0.9726 ±0.0041 (11) 

is consistent with the predicted value equal to f{m 2 ^/m 2 )/f{m e \/m? r ) = 0.972565 ± 0.000009, using lepton 
masses from (Particle Data Group, 2004). The uncertainty arises essentially from the r mass determination, 
m T = (1777.03t£i$) MeV, dominated by the BES result (BES Coll., 2002).' 
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FIG. 2 Measurements of the branching fraction for r — » e v^v T (left) and r — * p, v )1 v T (right). References for 
experiments are (ALEPH Coll., 2005; ARGUS Coll., 1993; CLEO Coll., 1997b; DELPHI Coll., 1999; L3 Coll., 2001; 
OPAL Coll., 1999a, 2003). The world averages are indicated by the shaded vertical bands. 



The result for B e can also be compared with the value obtained from the r and /i masses and lifetimes t TiAI , 
under the assumption of r—fi universality, according to 

^ _ r T /m T \ 5 f(m 2 Jm 2 T ) S^S^ ^ 

Tr 

(1 631.9 ± 1.4) fs ' 

Using the r lifetime (Particle Data Group, 2004), t t = (290.6 ± 1.1) fs, as well as the fi lifetime, one obtains 
B e = (17.807 ± 0.067)%, in excellent agreement with the direct measurements. 

The direct value for B e and the two derived values from B^ and t t are in good agreement with each other, 
providing a consistent and precise combined result for the electronic branching fraction, 

Bl ni = (17.818 ± 0.032)% . (13) 

Hence we find for the branching fraction of r — * v T hadrons, Shad, assuming lepton universality and unitarity 

B had = 1-^-0^ = 1- 1.97257 B e ini = (64.853 ± 0.063)% . (14) 



C. Decay modes with kaons 

Measurements of r decay modes with kaons in the final state are important on at least two counts; (i) 
to determine the Cabibbo-suppressed component of the hadronic decay width, leading to interesting physics 
applications (see Section VIII), and (ii) to perform the necessary corrections to the measured topological 
branching fractions, based on charged particles and ir° multiplicities. This latter step is made necessary 
by the difficulty separating charged pions and kaons, and identifying K®'s at the level of an event. The 
subtraction of strange final states is then performed on a statistical basis. 

The most complete measurements of r branching fractions involving kaons are from ALEPH (ALEPH Coll., 
1997a, 1998a, b, 1999a, b). From these results a comprehensive study of the strange sector has been performed, 
as well as some aspects of the dynamics in the final states with a KK pair. The experimental difficulty 
resides in the fact that both charged and neutral kaons have to be identified into multi-particle final states. 
The analyses involve the detection of charged kaons through the measurement of specific ionization losses and 
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TABLE I World average branching fractions (10 3 ) for r decays into strange final states and v T from (ALEPH Coll., 
1999b; CLEO Coll., 2003, 2005; OPAL Coll., 2004; Particle Data Group, 2004). The value for the (K3n)~ branching 
fraction takes into account measurements from ALEPH and CLEO, as well as estimates for unseen final states using 
isospin relations. The estimate for the very small (K4n)~ fraction, in addition to the measured K*~r\ rate (CLEO 
Coll., 1999a), is obtained from an empirical observation of pionic modes and Cabibbo suppression. 



Mode 


Z3(10" 3 ) 


K~ 


6.81 ±0.23 


K-7T° 


4.54 ± 0.30 




8.78 ±0.38 


K-TT°ir° 


0.58 ±0.24 


K°7v-7r° 


3.60 ±0.40 




3.30 ±0.28 


K~f] 


0.27 ±0.06 


(K3ir)~ (estimated) 


0.74 ± 0.30 


Ki(1270)- -> K~u) 


0.67 ±0.21 


(K4tk)~ (estimated) and K*~r\ 


0.40 ±0.12 


Sum 


29.69 ±0.86 



neutral kaons through cither their decay into 7r + 7r~ (Kg) or their interactions in a calorimeter (K °). Assuming 
CP invariance the K® and K® branching fractions can be averaged. The higher efficiency achieved for K® 
detection leads to smaller statistical uncertainties. Although the strange final states are Cabibbo-suppressed, 
the precision achieved in these measurements is sufficient for various significant tests, ranging from r— fi lepton 
universality to hadron dynamics, resonance production, QCD analyses, and the measurement of the Cabibbo 
angle. In addition, the consistency of the results is checked through the verification of isospin relations between 
the rates observed for different charge configurations of the same final states. 

New experimental results have been released by CLEO (CLEO Coll., 1999a,b, 2003, 2005) and OPAL (OPAL 
Coll., 2004) achieving significant progress in the Kixix modes. The new measurements tend to favor larger 
branching fractions. In particular, the branching fraction for the K~ir + Tr~ mode, (2.14 ± 0.47) x 10~ 3 for 
ALEPH, becomes in the average (3.30 ± 0.28) x 10~ 3 with a x 2 of 9.6 for 2 DF (about 1% probability). Also, 
the branching fraction for K°n-n°, (3.27 ± 0.51) x 10~ 3 for ALEPH, becomes (3.60 ± 0.40) x 10" 3 in the 
average. A compilation of the world average S = 1 branching fractions is given in Table I. 

The branching fraction for the decay t~ — > K~v, Bk, 

B k = (0.681 ±0.023)% , (15) 

can be compared to the theoretical prediction assuming r-/x universality in the charged weak current, Bj^ 1 , 
using the values (Particle Data Group, 2004) for the r and K lifetimes and masses, the K~ — ► branching 
fraction, and a small radiative correction (Decker and Finkemeier, 1993), 5 t /k = 1.0090 ± 0.0022, 

T T B K ^nu rn\ (I - m 2 K /rri 



t k 2mKmf 1 \l-rnf l /m z K ) ' 



(40 636 ± 162) fs 



(0.715 ±0.003)% . 



D. Nonstrange hadronic branching fractions 

Nonstrange hadronic branching fractions represent the largest part of r decays and the physics goals re- 
quire the achievement of small systematic uncertainties. The main experimental difficulty for this task is 
the photon/7r° reconstruction. The high collimation of t decays at LEP energies makes this task quite often 
difficult, since these photons are close to one another or close to the showers generated by charged hadrons. 



10 



Of particular relevance is the rejection of fake photons, which are due to hadronic interactions, fluctuations of 
electromagnetic showers, or an overlap of several showers. These problems reach a tolerable level in ALEPH 
thanks to the fine granularity of its electromagnetic calorimeter, comprising about 70,000 cells covering the 
detector angular acceptance, each cell being segmented threefold along the photon direction in order to discrim- 
inate between hadronic and electromagnetic showers. Even then, proper and reliable reconstruction methods 
had to be developed to correctly identify photon candidates. In the ALEPH analysis, a likelihood method 
is used for discriminating between genuine and fake photons, taking advantage of the detailed calorimeter 
information for each detected cluster of energy not associated to a charged track. Photons converted in the 
detector material are also included. Then 7r°'s are reconstructed: resolved 7r°'s from two-photon pairing, 
unresolved high-energy 7r°'s from merged clusters, and residual 7r°'s from the remaining single photons after 
removing radiative, bremsstrahlung and fake photons with a likelihood method. 

The 13 t decay classes defined in Section III. A are not fully exclusive as they contain contributions from 
final states involving kaons. They are coming from charged kaons included in the class definition using hadrons 
(h) and from neutral kaons not taken into account in the classification. These contributions are subtracted 
on a statistical basis using the branching fractions measured separately (Section III.C). 

The r decays involving r\ or u> mesons also require special attention because of their electromagnetic decay 
modes. Indeed the final state classification relies in part on the tt° multiplicity, thereby assuming that 
all photons — except those specifically identified as bremsstrahlung or radiative — originate from 7r° decays. 
Therefore the non-7r° photons from rj and ui decays are treated as ir° candidates in the ALEPH analysis 
and the systematic bias introduced by this effect must be evaluated. The corrections are based on specific 
measurements by ALEPH of r decay modes containing those mesons (ALEPH Coll., 1997a). Hence the 
final results correspond to exclusive branching fractions obtained from the values measured in the topological 
classification, corrected by the contributions from K, rj and lo modes measured separately taking into account 
through the simulation their specific selection and reconstruction efficiencies. This delicate bookkeeping takes 
into account all the major decay modes of the considered mesons (Particle Data Group, 2004), including the 
isospin-violating uj — ► 7r + 7r~ decay mode. The main decay modes considered are ir~u>, 7r~7r°a> and 7r~7r ?7 with 
branching fractions of (2.26 ±0.18) x 10~ 2 , (4.3 ±0.5) x 10" 3 , and (1.80 ±0.45) x 10" 3 (ALEPH Coll., 1997a), 
respectively. The first two values are derived from the branching fractions for the 27r~7r + 7r and 2n~ ir + 2ir° 
modes obtained in the global analysis and the measured u> fractions of 0.431 ± 0.033 from (ALEPH Coll., 
1997a) and the average value, 0.78 ±0.06, from (ALEPH Coll., 1997a) and (CLEO Coll., 1993b), respectively. 

Small contributions with r\ and ui have been identified and measured in (CLEO Coll., 1997a; Weinstein, 
2001) with the decay modes t~ — > v T i]ir~ n + tt~ (2.4 ± 0.5), r~ — > v T rjir~27r a (1.5 ± 0.5), t~ — ► ^ T a>7r~7r + 7r~ 
(1.2 ± 0.2), and t~ — ► ^ T a>7r _ 27r° (1.5 ± 0.5), where the branching fractions are given in units of 10~ 4 . For 
the sake of completeness these modes have been included. Another small correction has been applied to take 
into account the a\ radiative decay into ttj with a branching fraction of (2.1 ± 0.8) x 10~ 3 (Ziclinski, 1984). 

The final ALEPH results on the hadronic branching fractions for r decay modes not including kaons 
are given in Table II. The branching fractions obtained for the different channels are correlated with each 
other. On one hand the statistical fluctuations in the data and the Monte Carlo samples are driven by 
the multinomial distribution of the corresponding events, producing well-understood correlations. On the 
other hand the systematic effects also induce significant correlations between the different channels. All the 
systematic studies were done keeping track of the correlated variations in the final branching fraction results, 
thus allowing a proper propagation of errors. The full covariance matrices from statistical and systematic 
origins are given in (ALEPH Coll., 2005). 

The branching fraction B n for the tt^v decay mode, 



agrees well with a theoretical prediction only depending on the assumption of r-/x universality, S" m , using the 
values (Particle Data Group, 2004) for the r and n lifetimes and masses, the iv~ — > branching fraction, 

and a small radiative correction (Decker and Finkemeier, 1993), 6 T / V — 1.0016 ± 0.0014, 



= (10.83 ±0.11)% , 



(17) 




(18) 



(2 663.7 ± 4.0) fs 



= (10.910 ±0.044)% , 



where the error in the denominator of the prediction is dominated by the uncertainty on the structure- 
dependent radiative correction. 
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TABLE II Results (ALEPH Coll., 2005) on exclusive hadronic branching fractions for modes without kaons. The 
contributions from channels with r\ and u) are given separately, the latter one only for the electromagnetic to decays. 
The measured branching ratio for 27r~7r + 37r° is consistent with zero and the value listed in the table has been estimated 
with the help of CLEO results (CLEO Coll., 1993c). 



Mode B ± a stat ± a syst (in %) 



TV 


10.828 


± 


0.070 


± 


0.078 






25.471 




0.097 


-j- 


0.085 




7T~2tT° 


9.239 


± 


0.086 


± 


0.090 




-n 
TT oTT 


0.977 


± 


0.069 


± 


0.058 




7r _ 47T° 


0.112 


± 


0.037 


± 


0.035 




tt~tt~tt + 


9.041 


± 


0.060 


± 


0.076 




- - + 

TT TT TT TT 


4.590 


± 


0.057 


± 


0.064 




2tt-tt + 2tt" 


0.392 


± 


0.030 


± 


0.035 




2tt-tt + 3tt° 


0.013 


± 


0.000 


± 


0.010 


estimate 


3tt-2tt+ 


0.072 


± 


0.009 


± 


0.012 




3tt-2tt + tt" 


0.014 


± 


0.007 


± 


0.006 




TT~Tr"ri 


0.180 


± 


0.040 


± 


0.020 


(ALEPH Coll., 1997a) 


TT~2TT a ri 


0.015 


± 


0.004 


± 


0.003 


(CLEO Coll., 1997a; Weinstein, 2001) 


TT~ TT~ TT + 7] 


0.024 


± 


0.003 


± 


0.004 


(CLEO Coll., 1997a; Weinstein, 2001) 




0.040 


± 


0.000 


± 


0.020 


estimate 


TT~U>(^ 7T°7, TT + TT~) 


0.253 


± 


0.005 


± 


0.017 


(ALEPH Coll, 1997a) 


7r~7r () a;(^ 7r°7, tt + tt~) 


0.048 


± 


0.006 


± 


0.007 


(ALEPH Coll., 1997a; CLEO Coll., 1993b) 


tt~2tt°lu(^ 7r°7, tt + tt~) 


0.002 


± 


0.001 


± 


0.001 


(CLEO Coll., 1997a; Weinstein, 2001) 


TT~TT~TT + Uj(^ Tt"j, TT + TT~ ) 


0.001 


± 


0.001 


± 


0.001 


(CLEO Coll., 1997a; Weinstein, 2001) 



Concerning the {Ztt)~ v decay mode, which is dominated by the a\ resonance, it is interesting to compare the 
rates in the 27r~7T + and 7t~27t° channels. Dominance of the pu intermediate state in the a\ decay predicts equal 
rates, but (small) isospin-breaking effects are expected from different charged and neutral tt (and possibly 
p) masses, slightly favoring the 7r _ 27r° channel. However, a partial- wave analysis of the 7r"27r° final state 
performed by CLEO (CLEO Coll, 2000a) revealed that this decay is in fact rather complicated involving 
many intermediate states, in particular isoscalars, which amount to about 20% of the total rate and produce 
strong interference effects. A good description of the a\ decays was achieved in this study, which can be applied 
to the (37r) - final state. Including known isospin-breaking effects from the pion masses, it predicts a ratio of 
the 27r~7T + to the 7r~27r° rate of 0.985. This value is in agreement with the ALEPH measurement (ALEPH 
Coll, 2005) 

B J^^ + = 0.979 ±0.018 , (19) 

supporting the expected trend. 

A comparison of recent measurements of the dominant hadronic branching fractions is provided in Figs. 3. 
To compare with other experiments, results are given there without charged hadron identification, i.e., the 
branching fractions correspond to the sum of all branching fractions of the same charged particle topology, 
whether pions or kaons. 



E. Separation of vector and axial-vector contributions 

From the complete analysis of the t branching fractions presented in (ALEPH Coll, 2005) it is possible to 
determine the nonstrange vector (V) and axial-vector (A) contributions to the total r hadronic width. It is 
conveniently expressed in terms of their ratios to the leptonic width, denoted R T y and Rt,a, respectively, while 
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CLEO 97 






1 1.52±0.05±0.12 






OPAL 91 -95 

1 1 .98+0.1 3+0.16 






ALEPH 91-95 






1 1 .524±0.07±0.078 






average 






1 1 .58410.076 


11 


.5 


12 12.5 







CLEO 94 






25.87±0.12±0.42 






OPAL 91-95 






1 25.89±0.17±0.29 




- 


ALEPH 91-95 

25.924±0.097±0.085 






average 

25.91610.116 



B(hv) (%) 



25 26 27 28 29 
B(h7t°v) (%) 







Ill 

CLEO 






8.96±0.16±0.44 






ALEPH 91-95 






~" 9.295±0. 08410.088 






average 




- 


9.274±0.118 



8.5 9 9.5 10 10.5 11 
B(h2A) (%) 







CLEO 






9.51±0.07±0.2 






OPAL 91-95 






* 1 9.87±0.1±0.24 






ALEPH 91-95 






9.469±0.062±0.073 




'' — • — 1 


average 

9.516±0.083 



9.5 10 10.5 

B(3hv) (%) 



11 







CLEO 






4.23±0.06±0.22 






ALEPH 91-95 




1— • 


4.726±0.059±0.049 




- 


average 

4.676±0.073 



4.5 



5.5 







HRS 

0.051 ±0.02±0 






ARGUS 






0.064±0.023±0.01 






CLEO 






0.077±0.005±0.009 






OPAL 91-95 






* ' 0.091±0.014±0.006 






ALEPH 91-95 






0.072±0.009±0.012 




l-» 


BABAR 

H 0.0856±0. 000510.0042 




— •— 


average 

0.083±0.004 



B(3h7t°v) (%) 



0.025 0.05 0.075 0.1 0.125 0.15 

B(5hv) (%) 



FIG. 3 Comparison of results on the main branching fractions for r hadronic modes (ALEPH Coll., 2005; ARGUS 
Coll., 1988; BABAR Coll., 2005b; CLEO Coll., 1993a, 1994a,b, 1995a, 1997b; HRS Colla., 1987; OPAL Coll., 1998, 
1999c). The world averages are indicated by the shaded vertical bands. 
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the inclusive strange counterpart is denoted R T ,s- Only limited attempts to separate vector and axial- vector 
IAS 1 ! = 1 currents have been performed by the experiments so far (ALEPH Coll., 1999b). 

The ratio R T for the total hadronic width is calculated from the leptonic branching fractions alone, and 
eventually from the (massless) electronic branching fraction only 1 , 

R T = r(r ' - "ThadroBQ = 1 - B - g„ = 1 = 

T(r- -» v T e-V e ) B e B™ V ' 

for ALEPH, assuming universality in the leptonic weak current. All R values given below are rescaled 
so that the sum of the hadronic and the leptonic branching fractions, the latter one computed using the 
"universal" value B™ 1 (ALEPH), add up to one. The deviations introduced in this way are small, less than 
10% of the experimental error, but this procedure guarantees the consistency of all values. Using the ALEPH 
measurement of the strange branching fractions (ALEPH Coll., 1999b), supplemented by a (small) contribution 
from the K*~ij channel measured by CLEO (CLEO Coll., 1999a), one finds (see Footnote 1) 

B s = (2.85 ±0.11)% , 
R t S = 0.1603 ± 0.0064 , (21) 



and consequently 



B v +A = (62.01 ±0.14)% , 
Rt,v+a = 3.482 ± 0.014 , (22) 



for the nonstrange component. 

Separation of V and A components in hadronic final states with only pions is straightforward. Even number 
of pions has G = 1 corresponding to vector states, while an odd number of pions has G = — 1, which tags 
axial-vector states. This property places a strong requirement on the knowledge of the tt° reconstruction 
efficiency in the detector. 

Modes with a KK pair are not in general eigenstates of G-parity and contribute to both V and A channels. 
While the decay to K~K° is pure vector, the situation is a priori not clear in the KKtt channel, observed in 
three charged modes: K~ K + tt~ , K~K°tt° and K°K°tt~ . Three sources of information exist on the possible 
V / A content in these decays: 

1. In the ALEPH analysis of r decay modes with kaons (ALEPH Coll., 1999b), an estimate of the vector 
contribution was obtained using the available e + e~ annihilation data in the KKtt channel, extracted in 
the 1=1 state. This contribution was found to be small, yielding the branching fraction Z?cvc( T ~ ~~ * 



v t (KKtt) v ) = (0.26 ± 0.39) x 10 3 , which corresponds to an axial fraction of /a,cvc = 0.94^q 



.06 
0.08- 



2. The CVC result is corroborated by a partial- wave and lineshape analysis of the a\ resonance from r 
decays in the v t tt~2tt° mode in (CLEO Coll., 2000a). The observation through unitarity of the opening 
of the K*K decay mode of the cti is postulated and the branching fraction B(a\ — > K*K) = (3.3 ±0.5)% 
is derived. With the known decay rate of t~ — > v T a^ , it is easy to see that such a result saturates, and 
even exceeds, the total rate for the KKtt channel. The corresponding axial fraction is /a,3tt = 1.30±0.24. 

3. A new piece of information, also contributed by CLEO, but conflicting with the two previous results, 
was recently published (CLEO Coll., 2004). It is based on a partial-wave analysis in the K~K + tt~ 
channel using two-body resonance production and including many possibly contributing channels. A 
much smaller axial contribution is found here: f A = 0.56 ± 0.10. 

Since the three determinations are inconsistent, a conservative value of Ja — 0.75 ± 0.25 is assumed. 
It encompasses the range allowed by the previous results and still represents some progress over earlier 
analyses (ALEPH Coll., 1997b, 1998c), where a value of 0.5 ± 0.5 was used. For the decays into KKtttt no 
information is available in this respect, and a fraction of 0.5 ± 0.5 is used. 



1 The branching ratios B" ni (ALEPH) = (17.810 ± 0.039)% and (ALEPH) = (2.85 ± 0.11)% used in this section are from the 
ALEPH analysis only (ALEPH Coll., 2005) and differ slightly from the world average values given in Sections HLB and III.C. 
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The total nonstrange vector and axial- vector contributions obtained in (ALEPH Coll., 2005) are 



B v = (31.82 ± 0.18 ± 0.12)% , (23) 

R t V = 1.787 ± 0.011 ± 0.007 , (24) 

B A = (30.19 ± 0.18 ± 0.12)% , (25) 

R t A = 1.695 ± 0.011 ± 0.007 , (26) 



where the second errors reflect the uncertainties in the V/A separation in the channels with KK pairs. Taking 
care of the (anti)correlations between the respective uncertainties, one obtains the difference between the 
vector and axial- vector components, which is physically related to the amount of nonperturbative contributions 
in the nonstrange hadronic r decay width (see Section VII) 

B v -A = (1.63 ± 0.34 ± 0.24)% , (27) 
Rr.y-A = 0.092 ± 0.018 ± 0.014 , (28) 

where the second error has the same meaning as in Eqs. (23) and (25). The ratio 

R ^ y - A = 0.026 ± 0.007 (29) 

tir,V+A 

is a measure of the relative importance of nonperturbative contributions. Its smallness, in spite of our 
subtraction of observables with very different spectral contributions from hadronic resonances, represents a 
proof of the remarkable property that is referred to as quark-hadron duality 2 (Shifman et al., 1979). It is the 
basis of the QCD analyses of hadronic r decays presented in this review. 

The inclusive vector and axial-vector branching fractions provide the normalization of the corresponding 
spectral functions, which will be discussed below and which are the ingredients for the vacuum polarization 
calculations and the QCD analysis of hadronic r decays. 



IV. TAU HADRONIC SPECTRAL FUNCTIONS 



The t is the only lepton of the three-generation SM that is heavy enough to decay into hadrons. It is 
therefore an ideal laboratory for studying the charged weak hadronic currents and QCD. In analogy to the 
purely leptonic r decays, the invariant amplitude for the hadronic decays can be written in form of a factorized 
current-current interaction 



M(t- - hadrons" u T ) = ^L\V CKM \£^ , (30) 

v2 



with the corresponding nonstrange or strange CKM matrix element |Vckm|, an d where € M describes the 
leptonic V — A current 

^=V t ^(1-j 5 )t (31) 

of the weak interaction. The hadronic transition current is the piece of interest here. It probes the matrix 
element of the left-handed charged current between the vacuum and the hadronic final state. Restricting to 
•&V — A structure, one can write 

h„ = (hadrons 1^(0) - A M (0)|0> . (32) 
The differential r hadronic width can be expressed as follows 

dT(r- -» hadrons" u T ) = ^\V C km\ 2 L^H^ dPS , (33) 



2 Reference (Blok et al., 1998) gives the following definition of quark-hadron duality. "In a nutshell, a truncated OPE is 
analytically continued, term by term, from the Euclidean to the Minkowski domain. A smooth quark curve obtained in this 
way is supposed to coincide at high energies (energy releases) with the actual hadronic cross section." 
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with the leptonic (hadronic) tensor L^ u (H^ v ) and the Lorentz invariant phase space element dPS. The 
hadronic tensor obeys a Lorentz invariant decomposition, the dynamical structure being embedded in so- 
called structure functions (Kiihn and Mirkes, 1992). 

For the simplest hadronic decay modes t~ — > h~u T (h = tt,K), the structure functions reduce to 5- 
distributions, whereas resonances appear in r decays into multiple hadrons. The hadronic structure is de- 
scribed by spectral functions. The quantum number corresponding to the (conserved) vector and axial-vector 
currents is the isotopic G-parity. It is a generalized multiplicative symmetry of multi-pion systems under the 
successive operations G and R, where G conjugates the charge and R rotates the system around the second 
axis in the I — 1 isospace. The measurement of the r vector and axial-vector current spectral functions 
requires the selection and identification of t decay modes with a defined G-parity G = +1 and G = — 1, 
and hence hadronic channels with an even and odd number of neutral or charged pions, respectively. Since 
hadronic final states of different G-parity differ also in their J p quantum numbers, there is no interference 
between these two states. Hence the total hadronic width separates into r to t = Ty + Ty*. 



A. Definitions 



The spectral function v\ (ai, an), where the subscript refers to the spin J of the hadronic system, is defined 
for a nonstrange (|AS| = 0) or strange (|A5| = 1) vector (axial-vector) hadronic t decay channel V~v T 
(A~v T ). The spectral function is obtained by dividing the normalized invariant mass-squared distribution 
(l/N V / A )(dN V / A /ds) for a given hadronic mass y/s by the appropriate kinematic factor 



v 1 (s)/a 1 (s) = 



6 |Vckm| 2 Sew B(t 



dN, 



V/A 



N v/A ds 



B{t~ - 


+ V-/A-v T ) 




B{t~ 


— > e~ v e y T ) 




s 




_1 


m 2 T/ 


1 



(34) 



a (s) 



G\V C 



KM 



'EW 



B(t- -> n-(K-)v T ) dN A 
B(t~ — > e~ v e v T ) Na ds 



1 



mi 



(35) 



where Sew accounts for electroweak radiative corrections (Marciano and Sirlin, 1988) (see Section V.A for 
the relevant values of Sew)- Since CVC is a very good approximation for the nonstrange sector, the J = 
contribution to the nonstrange vector spectral function is put to zero, while the main contributions to ao are 
from the pion or kaon poles, with dN A /NA ds = S(s — m^ K ). They are connected through partial conservation 
of the axial-vector current (PCAC) to the corresponding decay constants, f n .K- The spectral functions are 
normalized by the ratio of the vector/axial- vector branching fraction B(t~ — > V~ / A~ v T ) to the branching 
fraction of the massless leptonic, i.e., electron, channel (discussed in Section III.B). 

The values for the CKM matrix elements V u d and V us have been recently much debated, because the 
determinations from nuclear (3 decays and semileptonic kaon decays did not satisfy the unitarity relation in 
the first row of the CKM matrix, that is \V u d\ 2 + \V US \ 2 + \V u t,\ 2 = 1 (for this comparison, the V u b element 
has a negligible contribution). The situation has improved with the availability of new results from Kg 3 and 
K^2 decays with revisited corrections from chiral perturbation theory on one hand, and more precise radiative 
corrections for the nuclear [3 decays (Marciano, 2005) on the other hand. The current results (Isidori, 2005), 
which are still preliminary 3 , 



\V ud \ = 0.9739 ±0.0003 
\V US \ = 0.2248 ±0.0016 



nuclear (3 decays , 



(36) 
(37) 



are consistent with CKM unitarity and provide a new level of precision in this sector. The \V u d\ value gives 
IK* | = 0.2269 ± 0.0013 if CKM unitarity is used. 



3 We thank V. Cirigliano and G. Isidori for their advice on the current best values for \V n J\ an d \V U s\- 
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Using unitarity and analyticity, the spectral functions are connected to the imaginary part of the two-point 
correlation (or hadronic vacuum polarization) functions 

n^fa) = ^ J d 4 xe^(0|T(C/(;.(x)^.(0)t)|0) 

= (-.^V + ^uiq 2 ) + g'V U%(q 2 ) (38) 

of vector (C/^- = Vfi = qjj^qi) or axial-vector (U£ — = 9^7^75^) color-singlet quark currents, and for 
time-like momenta-squared q 2 > 0. Lorentz decomposition is used to separate the correlation function into 
its J = 1 and J = parts. The polarization function II^ I/ c/ (s) have a branch cut along the real axis in 
the complex s = q 2 plane. Their imaginary parts give the spectral functions defined in (34), for nonstrange 
(strange) quark currents 

Imn S s) ,y/>) = ^Vai(s) , lmn^ (s) A ( S ) = i-ao(s) , (39) 

which provide the basis for comparing short-distance theory with hadronic data. 
The analytic vacuum polarization function ll[j\j(q 2 ) obeys the dispersion relation 

1 °? Inin^ (s) 

n, ( , J) r r(<7 2 ) = - / ds iLEll + subtractions , (40) 

J ' u 7r J s — q l — ie 



where the unknown but in general irrelevant subtraction constants can be removed by taking the derivative 
o{Ilij i u(q 2 ). The dispersion relation allows one to connect the experimentally accessible spectral functions to 
the correlation functions n\jjj(q 2 ), which can be derived from theory (QCD). 



B. The physical mass spectra and their unfolding 

The measurement of the r spectral functions defined in Eq. (34) requires the determination of the invariant 
mass-squared distributions, obtained from the experimental distributions after unfolding from the effects of 
measurement distortion. The unfolding procedure used by the ALEPH collaboration (ALEPH Coll., 1997b, 
1998c) follows a method based on the regularized inversion of the simulated detector response matrix using 
the Singular Value Decomposition (SVD) technique (Hocker and Kartvelishvili, 1996). The regularization 
function applied 4 minimizes the average curvature of the distribution. The optimal choice of the regularization 
strength is found by means of the Monte Carlo simulation where the true distribution is known and chosen 
to be close to the expected physical distribution. The OPAL collaboration (OPAL Coll., 1999b) employs a 
similar regularized unfolding technique (Blobel, 1984), which is however based on cubic splines instead of the 
binned histograms used by ALEPH. Due to the superior mass resolution, the CLEO collaboration considers 
it to be sufficient to only apply a bin-to-bin migration correction (CLEO Coll., 2000d). 

To measure exclusive spectral functions, individual unfolding procedures with specific detector response 
matrices and regularization parameters are applied for each r decay channel considered. An iterative procedure 
is followed to correct the Monte Carlo spectral functions used to subtract the cross-feed among the modes. 

Before unfolding the mass distributions, the r and non-r backgrounds are subtracted. In the case of r 
feed-through the Monte Carlo distributions normalized to the measured branching fractions are used. When 
the spectral functions are measured for the nonstrange exclusive final states, the contributions from strange 
modes classified in the same topology are subtracted using their Monte Carlo spectral functions normalized 
by the measured branching fractions. 

The systematic uncertainties affecting the decay classification of the exclusive modes are contained in the 
systematic errors on the branching fractions. Additional systematic uncertainties related to the shape of the 
unfolded mass-squared distributions, and not its normalization, are included. They are dominated by the 
photon and n° reconstruction. 



4 Regularization is necessary since in certain cases the unfolding produces results with unphysical behavior: statistically insignif- 
icant components in the detector response matrix can lead to large oscillations in the unfolded distribution. The regularization 
is achieved by applying a smooth damping function to the unfolded distribution. 
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FIG. 4 Left hand plots: the inclusive vector (upper) and axial-vector (lower) spectral functions as measured in (ALEPH 
Coll., 2005). The shaded areas indicate the contributing exclusive r decay channels. The curves show the predictions 
from the parton model (dotted) and from massless perturbative QCD using a s (Ml) = 0.120 (solid). Right hand plots: 



comparison of the inclusive vector (upper) and axial-vector (lower) spectral functions obtained by ALEPH and 
OPAL (OPAL Coll., 1999b). 



C. Inclusive nonstrange spectral functions 

1. Vector and axial-vector spectral functions 



The inclusive r vector and axial- vector spectral functions are shown in the upper and lower plots of Fig. 4, 
respectively. The left hand plots give the ALEPH result (ALEPH Coll., 1997b, 1998c, 2005) together with its 
most important exclusive contributions, and the right hand plots compare ALEPH with OPAL (OPAL Coll., 
1999b). The agreement between the experiments is satisfying with the exception of the 7r~7r° (vector) mode, 
where some discrepancies occur. We note that the branching fraction used for this mode to obtain the inclusive 
spectral function is not the same in both experiments. The branching fraction used by ALEPH (ALEPH Coll., 
2005) exceeds the one used by OPAL (OPAL Coll., 1999b) by 0.9% (relative). 

The curves in the left hand plots of Fig. 4 represent the parton model (leading-order QCD) prediction 
(dotted) and the improved (massless) perturbative QCD prediction (solid), assuming the relevant physics to 
be governed by short distances (Section VII). The difference between the two curves is due to higher order 
terms in the strong coupling (a s (s)/7r)™ with n = 1, 2, 3. At high energies the spectral functions are assumed 
to be dominated by continuum production, which locally agrees with perturbative QCD. This asymptotic 
region is not yet reached at s = rn^ for the vector and axial-vector spectral functions. 
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FIG. 5 Inclusive vector plus axial-vector (left) and vector minus axial-vector spectral function (right) as measured 
in (ALEPH Coll., 2005) (dots with errors bars) and (OPAL Coll., 1999b) (shaded one standard deviation errors). The 
lines show the predictions from the parton model (dotted) and from massless perturbative QCD using a s (Mf ) = 0.120 
(solid). They cancel to all orders in the difference. 



To build the vector spectral function ALEPH has measured the two- and four-pion final states exclusively, 
while only parts of the six-pion state. The total six-pion branching fraction has been determined using 
isospin symmetry (ALEPH Coll., 1997b). However, one has to account for the fact that the six-pion channel 
is contaminated by isospin-violating t~ — > ^{Stt^i/t decays, as reported by CLEO (CLEO Coll., 1997a; 
Weinstein, 2001). A small fraction of the lott~v t decay channel that is not reconstructed in the four-pion final 
state is added using the simulation. Similarly, one corrects for the ?77r"7r ^ r decay mode where the r\ decays 
into pions. For the r\ — > 27 mode, the r decay is classified in the h^3Tr v T final state, since the two-photon 
mass is inconsistent with the ir° mass and consequently each photon is reconstructed as a ir°. The K~ K°v T 
mass distribution is taken entirely from the simulation (recall that this mode is pure vector). The spectral 
functions for the final states KKir and KKtttt (see Section III.E for the discussion on the G-parity of these 
modes) are obtained from the Monte Carlo simulation, where large systematic errors are applied to cover this 
approximation (ALEPH Coll., 1997b). 

In analogy to the vector spectral function, the inclusive axial-vector spectral function is obtained by sum- 
ming up the exclusive axial- vector spectral functions with the addition of small unmeasured modes taken from 
the Monte Carlo simulation. The small fraction of the un~ir Q v T decay channel that is not accounted for in 
the 27r - 7r + 2'7r c V T final state is added from the simulation. Also considered are the axial-vector ?7(37r) _ i/ r final 
states (CLEO Coll., 1997a; Weinstein, 2001). CLEO observed that the dominant part of this mode issues 
from the r~ — * /i(1285)7r~^ r intermediate state, with B{t~ — > f\ir~v T ) = (0.068 ± 0.030)%, measured in 
the /1 — ► t]tt + tt^ and fi — > t]2tt° decay modes (CLEO Coll., 1997a; Weinstein, 2001). Since the /1 meson is 
isoscalar, the branching fractions relate as B{t~ — > ry7r _ 7r + 7r _ i/ r ) = 2 x B{t~ — > r\Tx^2-K a v T ). The distributions 
are taken from the ordinary six-pion phase space simulation accompanied by large systematic errors. The 
KKn final states contribute mostly to the inclusive axial-vector spectral function, with errors that are fully 
anticorrelated with the inclusive vector spectral function. The invariant mass distributions for these channels 
are taken from the simulation. 



2. Inclusive V ± A spectral functions 

For the total v\ + a\ hadronic spectral function it is not necessary to distinguish in the experiment whether 
a given event belongs to one or the other current. The one, two and three-pion final states dominate and 
their exclusive measurements are added with proper accounting for the (anti)correlations. The remaining 
contributing topologies are treated inclusively, i.e., without separation of the vector and axial- vector decay 
modes. This reduces the statistical uncertainty. The effect of the feed-through between r final states on the 
invariant mass spectrum is described by the Monte Carlo simulation and resolution effects are corrected by 
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the data unfolding. In this procedure the simulated mass distributions are iteratively corrected using the 
exclusive vector/axial- vector unfolded mass spectra. Also, one does not have to separate the vector/axial- 
vector currents of the KKir and KKtttt modes. The v\ + 04 spectral functions for ALEPH and OPAL are 
plotted in the left hand plot of Fig. 5. The improvement in precision when comparing to a sum of the two 
parts in Fig. 4 is significant at higher mass-squared values. 

One nicely identifies the oscillating behavior of the spectral function and it is interesting to observe that, 
unlike the vector/axial-vector spectral functions, it does approximately reach the asymptotic limit predicted 
by perturbative QCD at s — > m^. Also, the V + A spectral function, including the pion pole, exhibits the 
features expected from global quark-hadron duality: despite the huge oscillations due to the prominent ir, 
p(770), a\ and p(1450) resonances, the spectral function qualitatively averages out to the quark contribution 
from perturbative QCD. 

In the case of the v\ — a\ spectral function, uncertainties on the V/A separation are reinforced due to their 
anticorrelation. In addition, anticorrelations in the branching fractions between r final states with adjacent 
numbers of pions increase the errors. The V\ — a x spectral functions for ALEPH and OPAL are shown in the 
right hand plot of Fig. 5. The oscillating behavior of the respective v\ and a\ spectral functions is emphasized 
and the asymptotic regime is not reached at s = rrq.. However again, the strong oscillation generated by the 
hadron resonances mostly averages out to zero, as predicted by perturbative QCD. 



D. Strange spectral functions 

The total |A5| = 1 spectral functions from the ALEPH and OPAL analyses are shown in Fig. 6. For 
ALEPH, the spectra for the Kir and Kirir modes are obtained from the corrected mass spectra, taking into 
account acceptance and bin migration corrections. The spectrum for the Kir mode is taken from the K~ir° 
and K°ir~ final states because of their superior mass resolution as compared to the higher statistics K®ir state 
where the K ° detection suffers from poor energy resolution. For the other modes a phase-space generator is 
used to simulate the mass distributions. The contributions from the respective channels are normalized to 
the corresponding branching fractions. 

Since OPAL normalizes each channel using the corresponding world average branching fraction, which is 
dominated by the ALEPH measurement, the added information is on the shape, however limited by the 
relatively low statistics. As a consequence the OPAL and ALEPH spectral functions are correlated. 

The K*(892)~ resonance stands out clearly on the low mass side, while within a large uncertainty (domi- 
nated by background subtraction in the (Kirir) ~ modes) the higher energy part is consistent with the parton 
model (or perturbative QCD) expectation. 



V. COMPARING TAU SPECTRAL FUNCTIONS WITH ELECTRON-POSITRON ANNIHILATION DATA VIA 
THE CONSERVED VECTOR CURRENT 

The vacuum polarization functions (38) for the various types of quark currents ij carry the dynamical 
information for the theoretical evaluation of total cross sections and decay widths. In the limit of isospin 
invariance the vector current is conserved (CVC), so that the spectral function of a vector r decay mode 
X~v T in a given isospin state for the hadronic system is related to the e + e~ annihilation cross section of the 
corresponding isovector final state X Q 

vlTt-^xois) = ^fv hX -(s), (41) 
where a is the electromagnetic fine structure constant. 



A. Isospin breaking and radiative corrections 

Since the breaking of isospin symmetry is expected at some level, in particular from electromagnetic effects, 
it is useful to carefully write down all the factors involved in the comparison of r and e + e~ spectral functions 
in order to make explicit the possible sources causing the breakdown of CVC. For the dominant irir spectral 
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FIG. 6 Inclusive hadronic vector plus axial-vector spectral function from r decays into S = ±1 final states with its 
different contributions as indicated from (ALEPH Coll., 1999b) and (OPAL Coll., 2004). The kaon pole is not included 
in the ALEPH plot. The parton model prediction is given by the dashed line. 



functions, we have on the e + e side 



< + - + -\ 4na2 i \ 
a(e e — > it n ) = vq(s) , 

s 

where the spectral function vq(s) is related to the pion form factor F®(s) by 



vo{s) 



12 



(42) 



(43) 



and where /3q(s) is the threshold kinematic factor. Similarly, on the r side the relation between spectral 
function (34) and the charged pion form factor reads 



v -( s ) 



12 



Isospin symmetry implies v~(s) = vq(s). The threshold functions /?o,- are defined by 
where 



f3(s,mi,m 2 ) = 



1 - 



(mi + m 2 )' 



1 - 



(mi - m 2 f 



1/2 



(44) 
(45) 

(46) 
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The experiments measure r decays inclusively with respect to radiative photons, i.e., for t~ — > ^ T 7r~7r°(7). 
The measured spectral function is hence w* (s) = V-(s) G(s), where G(s) is a radiative correction. 
Several levels of SU(2) breaking can be identified: 

• Electroweak radiative corrections to t decays are contained in the Sew factor (Braaten and Li, 1990; 
Marciano and Sirlin, 1988), which is dominated by short-distance effects. It is expected to be weakly 
dependent on the specific hadronic final state, as verified for the t~ — > (n~ , K~)v T decays (Decker 
and Finkemeier, 1993). Detailed calculations have been performed for the 7r~7r° channel (Cirigliano et 
al, 2001, 2002), which also confirm the relative smallness of the long-distance contributions. The total 
correction is Sew = S'ew^em/^em' where S^ w is the leading-log short-distance electroweak factor 
(which vanishes for leptons) and S^ 1 ' 1015 are the nonleading electromagnetic corrections. The latter 
corrections have been calculated at the quark level (Braaten and Li, 1990), at the hadron level for the 
-k~tt° decay mode (Cirigliano et at, 2001), and for leptons (Braaten and Li, 1990; Marciano and Sirlin, 
1988). The total correction amounts to (Davier et at, 2003a) S\g$ = 1.0198 ± 0.0006 for the inclusive 
hadron decay rate and S^w — (1.0232 ± 0.0006) • G E ^ (s) for the tttt decay mode, where G^m (s) is an 
s-dependent long-distance radiative correction (Cirigliano et al, 2001). 

• The pion mass splitting breaks isospin symmetry in the spectral functions (Alcmany et al, 1998; Czyz 
and Kiihn, 2001) since /3_(s) + (3 (s). 

• Isospin symmetry is also broken in the pion form factor due to the n mass splitting (Alemany et al, 
1998; Cirigliano et al, 2001). 

• A similar effect is expected from the p mass splitting. The theoretical expectation (Bijnens and Gos- 
dzinsky, 1996) gives a limit (< 0.7 MeV), but this is only a rough estimate. Hence the question must be 
investigated experimentally, the best approach being the explicit comparison of t and e+e~ 2-7T spectral 
functions, after correction for the other isospin-brcaking effects. No correction for p mass splitting is 
applied initially. 

• Explicit electromagnetic decays such as nj, 777, £ + l~ and irirj introduce small differences between the 
widths of the charged and neutral p's. 

• Isospin violation in the strong amplitude through the mass difference between u and d quarks is expected 
to be negligible. 

• When comparing r with e + e~ data, an obvious and locally large correction must be applied to the 
t spectral function to introduce the effect of p-uj mixing, only present in the neutral channel. This 
correction is computed using the parameters determined by the e + e~ experiments in their form factor 
fits to the 7r + 7r~ lincshape modeling p-oj interference (CMD-2 Coll., 2004). 

B. Tau spectral functions and e + e data 

Data from r decays into two- and four-pion final states t~ — > v T 7r~7r°, t~ — > v T 7r~37r° and t~ — > v T 
2tt-tt+tt°, are available from ALEPH (ALEPH Coll., 2005), CLEO (CLEO Coll., 2000c,d) and OPAL (OPAL 
Coll., 1999b). For the two-pion final states, they are compared in Fig. 7 in the region around the p resonance. 
For this comparison, each data set is normalized to the world average branching fraction and plotted with 
respect to their weighted average. The two most precise results from ALEPH and CLEO are in agreement. 
The statistics are comparable in the two cases, however due to a flat acceptance in ALEPH and a strongly 
increasing one (with the mass-squared) in CLEO, ALEPH data are more precise below the p peak, while 
CLEO has the better precision above. In the following we will use the r~ — ► ir~-K°v T spectral function 
averaged over the three experiments. 

The exclusive low-energy e + e~ cross sections have been mainly measured by experiments running at e + e~ 
colliders in Novosibirsk and Orsay. The most precise e + e~ — ► 7r + 7r~ measurements come from CMD-2, which 
arc now available in their final form (CMD-2 Coll., 2004), after a significant revision where problems related 
to radiative corrections have been removed. The overall systematic error of the final data is quoted to be 
0.6% and is dominated by the uncertainties in the radiative corrections (0.4%). 

Recently, new data on the 7r + 7r~ spectral function in the mass region between 0.60 and 0.97 GeV 
were presented by the KLOE collaboration (KLOE Coll., 2005), using the — for the purpose of precision 
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FIG. 7 The left figure shows a relative comparison of the -k + tt~ spectral functions extracted from r data from the 
different experiments, expressed as a ratio to the average r spectral function. The right figure emphasizes the p region. 
For CLEO only the statistical errors are shown. 

measurements — innovative technique of the radiative return (Binncr et al, 1999; Czyz et ai, 2002). The 
statistical precision of these data by far outperforms the Novosibirsk sample, but the systematic errors are 
about twice as large as those assigned by CMD-2. 

C. Comparing r with e + e _ data 

1. The 2-7T spectral function 

The 2-7T spectral functions extracted from r decays and e + e~ annihilation data are compared in Fig. 8. The 
e+e- data are taken from (CMD-2 Coll., 2004; DM1 Coll., 1978; DM2 Coll., 1989; OLYA Coll., 1979; OLYA, 
CMD Colls., 1985; TOF Coll., 1981). All error bars shown contain statistical and systematic errors. Visually, 
the agreement appears to be satisfactory, however the large dynamical range involved does not permit an 
accurate test. To do so, the e + e _ data are plotted as a point-by-point ratio to the r spectral function in Fig. 9. 
In these latter plots we also show more recent e + e~-annihilation results from the SND collaboration (SND-2 
Coll., 2005), and from the radiative return analysis performed by the KLOE collaboration (KLOE Coll., 2005). 
Several observations can be made: 

• A significant discrepancy, mainly above the p peak is found between t and the e + e~ data from CMD-2 
as well as older data from OLYA. 

• Overall, the KLOE data seem to confirm the trend exhibited by the other (older) e + e~ data. 

• Some disagreement between KLOE and CMD-2 occurs on the low mass side (KLOE data are large), on 
the p peak (KLOE below CMD-2) as well as on the high mass side (KLOE data are low). 

• A significant discrepancy between the (most recent) SND data and KLOE is observed in the energy 
domain above the p(770) peak. The SND data largely dissolve the discrepancy observed with the r 
data. 

At this stage, the r spectral function has not been corrected for a possible p~-p° mass and width split- 
ting (ALEPH Coll., 2005; Davier, 2003; Ghozzi and Jegerlehner, 2004). In contrast to earlier experimen- 
tal (ALEPH Coll., 1997b) and theoretical results (Bijnens and Gosdzinsky, 1996), a combined pion form 
factor fit to the new precise data on r spectral functions and e + e~ leads to m p - — m p a — (2.4 ± 0.8) MeV 
(Section V.D.2), while no significant width splitting is observed within the fit error of 1.0 MeV. 

Considering the mass splitting in the isospin-breaking correction of the r spectral function tends to locally 
improve though not restore the agreement between r and CMD-2 data, leaving an overall normalization 
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FIG. 8 Comparison of the -k + tt~ spectral functions from isospin-breaking corrected r data (world average) expressed 
as e + e~ cross sections, and e + e~ annihilation. The band indicates the combined r and e + e~ result within 1<t errors. 
It is given for illustration purpose only. A compendium of references for the e + e~ data is given in (Davier et al, 
2003a,b). 



discrepancy. Increasing the T p - — r p o width splitting by +3 MeV improves the agreement between r and 
KLOE data in the peak region, while it cannot correct the discrepancies in the tails. Note that a correction 
of the mass splitting alone would increase the discrepancy between the r and e + e~-based results for a h ad > LO 
(Section VI.F). 

In light of the new SND data (SND-2 Coll., 2005), it seems appropriate to consider the possibility of a bias 
in the KLOE results, and one may also question the very small systematic uncertainties claimed by CMD-2. 
We also point out that the application of the r long-distance correction G^m ( s ) (Cirigliano et al., 2001) 
worsens the agreement with the SND data in the energy domain above the p(770) peak. 



2. The 4-7T spectral functions 



The spectral function measurements of the t vector-current final states it 3ir° and 2ir 7r + 7r° are compared 
to the cross sections of the corresponding e + e~ annihilation into the isovector states 27r~27r + and 7r~7r + 27r°. 
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FIG. 9 Relative comparison of the isospin-breaking corrected r data (world average) and 7r + 7r~ spectral functions 
from e + e~ annihilation, expressed as a ratio to the r spectral function. The shaded band gives the uncertainty on the 
r spectral function. The right hand plot zooms into the p(770) peak region where the data density is high. The p-ui 
mixing has been phenomenologically corrected using a fit to the CMD-2 data. The larger mixing observed in the SND 
data, expressed as a 40% increase in the value for the u> — » tt + ty~ branching fraction (SND-2 Coll., 2005), is responsible 
for the residual effect after correction seen in the right hand plot. 
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FIG. 10 Comparison of the 27r + 27r (left) and it + tt 2ty° (right) spectral functions from e + e and isospin-breaking 
corrected r data (ALEPH), expressed as e + e~ cross sections. References for the e + e~ data are given in (Davier et al, 
2003a,b). The BABAR results are taken from (BABAR Coll., 2005a). 



Using Eq. (35) and isospin invariance, the following relations hold 

I=1 Ana , 

1=1 47ra 2 r -, , . 

Oe+e-^ + Tr-^OTrO — — 2tt—k+tv° v t — v l^- 3tt° v t \ • ( 4 <5j 

The comparison of these cross sections is given in Fig. 10 for 27r + 27r~ (left hand plot) and 7T + 7r _ 27r° (right hand 
plot). The latter mode suffers from discrepancies between the results from the various e + e~ experiments. The 
r data, combining two measured spectral functions according to Eq. (48), and (coarsely) corrected for isospin 
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TABLE III Branching fractions of r vector decays into 2 and 4 pions in the final state (ALEPH Coll., 2005). Second 
column: r decay measurements. Third column: results inferred from e + e~ spectral functions using the isospin 
relations (41,44,47) and correcting for isospin breaking (Davier et ai, 2003a,b). The 2 pion prediction does not 
include the recent KLOE and SND data. Experimental errors, including uncertainties on the integration procedure, 
and theoretical errors (missing radiative corrections for e + e~, and isospin-breaking corrections and \V u d\ for r) are 
shown separately. Last column: differences between the direct measurements in r decays and the CVC evaluations, 
where the various errors have been added in quadrature. 



Mode 



T 



Branching fractions (in %) 

e+e-viaCVC A(t - e+e") 



v t -k~-k° 25.47 ±0.13 24.52 ±0.26 CX p±0.11 rad ±0.12 su(2) +0.95 ± 0.33 



0.31 

t~ ->■ v^-i-K 0.98 ± 0.09 1.09 ± 0.06 cxp ± 0.02 rad ± 0.05 SU (2) -0.11 ±0.12 



0.08 

t~ -> u T 2-K-Tx+Tx a 4.59 ± 0.09 3.63 ± 0.19 cxp ± 0.04 rad ± 0.09 SU (2) ±0.96 ±0.23 



breaking originating from the pion mass splitting (Czyz and Kiihn, 2001), lie somewhat in between, with large 
uncertainties above 1.4 GeV because of the lack of statistics and a large feed-through background in the mode 
with three neutral pions. In spite of these difficulties the tt~3tt° spectral function is in agreement with e + e~ 
data as can be seen in Fig. 10. Due to the inconsistencies among the e + e~ experiments a quantitative test of 
CVC in the 7r + 7r~27r° channel is premature. 



3. Branching fractions in r decays and CVC 

It is instructive to compare the r and e + e~ spectral functions in a more quantitative way by calculating 
weighted integrals over the mass range of interest up to the r mass. One convenient choice is provided by the r 
branching fractions, which for the spin one part involve as a weight the kinematic factor (1 — s/m^) 2 (l±2s/m^) 
coming from the V — A charged current in r decay. It is then possible to directly compare the measured r 
branching fractions to their prediction through isospin invariance (CVC), with the e+e~ isovector spectral 
functions as input. 

Using the universality- improved branching fraction (13), one finds the results given for the dominant chan- 
nels in Table III. The errors quoted for the CVC values are split into uncertainties from (z) the experimental 
input (the e+e~ annihilation cross sections) and the numerical integration procedure, (ii) additional radiative 
corrections applied to some of the e + e~ data (Davier et ai, 2003a), and (Hi) the isospin-breaking corrections 
when relating r and e + e~ spectral functions. 

As expected from the preceding discussion, a discrepancy is observed for the t~ — > tt~tt°v t branching 
fraction, with a difference of (0.95 ± 0.13 r ± 0.26 cc ± 0.11 ra d ± 0.12 SU (2))%, where the uncertainties are 
from the r branching fraction, e + e~ cross section, e + e~ missing radiative corrections and isospin-breaking 
corrections (also including the uncertainty on |Vud|), respectively. Adding all errors in quadrature gives a 2.9a 
effect. Including the new SND measurements (SND-2 Coll., 2005) into the 7r + 7r~ data sample would decrease 
this discrepancy to approximately 2.5c. The improvement would have been stronger, if we had discarded 
older data sets. In effect, using only SND data in the region where these are available [0.39-0.97 GeV], and all 
other data elsewhere, one finds for the predicted branching fraction: Bcvc( T ~ — * v T ^^ Q ) = (25.12 ± 0.36)%, 
which is in agreement with the r result 5 . More information on this comparison is displayed in Fig. 11, where 
we also give the results for using only CMD-2 and only KLOE data, respectively, in the available energy 



5 However, we point out that the scarce SND data points in the tails of the p(770) resonance, where the line shape is concave, 
may lead to an overestimate of the integral computed by the trapezoidal rule as it is done here. 
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FIG. 11 The measured branching fractions for r" -» 7r _ 7rV (ALEPH Coll., 2005; CLEO Coll., 1994a; OPAL Coll., 
1998) compared to the predictions from the e + e~ — > 7r + 7r~ spectral functions, applying the isospin-breaking correction 
factors discussed in Section V.A. For the e + e~ results we have used only the data from the indicated experiments 
where they are available (0.61-0.96 GeV for CMD-2, 0.60-0.97 GeV for KLOE, and 0.39-0.97 GeV for SND), and all 
data combined in the remaining energy domains below m T . Since the compatibility between the e + e~ data is marginal 
(c/. Fig. 9), we have refrained from taking their average. 



regions (and the average of all other data elsewhere). The discrepancy between the e + e _ -based prediction 
using the KLOE data and the r result is at the 3.8a level. 

The situation in the Att channels is different. While there is agreement for the ir~3ir° mode within a relative 
accuracy of 12%, the comparison is not satisfactory for 2tt~tt + tt°. In the latter case, the relative difference is 
very large, (23 ± 6)%, compared to a reasonable level of isospin symmetry breaking. As such, it rather points 
to experimental problems that require further investigation. These are emphasized by the scatter observed 
among the different e + e~ results. 



D. Fits to the nn spectral function 

Phenomenological fits to the pion form factor have been performed by the ALEPH and CLEO collabora- 
tions (ALEPH Coll., 1997b, 2005; CLEO Coll., 2000d). Appropriate parameterizations, e.g., the one proposed 
in (Kiihn and Santamaria, 1990), employ relativistic P-wave Breit-Wigner propagators, and coherent sum- 
mation of the resonance amplitudes using the isobar model. The 7T7t spectral function is dominated by the 
broad p resonance, parameterized by both collaborations following Gounaris-Sakurai (Gounaris and Sakurai, 
1968) (GS). The GS parameterization takes into account analyticity and unitarity properties. We will in the 
following discuss the fits performed in (ALEPH Coll., 2005). 

Assuming vector dominance, the pion form factor is given by interfering amplitudes from the known isovector 
meson resonances p(770), p(1450) and /o(1700) with relative strengths 1, /3, and 7. Although one could expect 
from the quark model that (3 and 7 are real and respectively negative and positive, the phase of (3, 4>p, is left 
free in the fits, and only the much smaller 7 is assumed to be real for lack of precise experimental information 
at large masses. Taking into account p-u mixing, one writes for the 7r + 7r~ form factor in e + e~ annihilation 



„/ =ll o, , _ BW p(770) ( S ) gg^gg + /3BW p(1450) ( 5 ) + 7BW p(1700) ( S ) 
K [S) ~ 1 + (3 + 7 



with the Breit-Wigner propagators 

/p(m " )V " ; " ml - s + f(s) - i^T p (s) 



„ m„( 1 + d ■ r„/m„l 
BW«s (s) = — p \ p J- -, , (50) 
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where 



/(*) 



fc 3 (m p ) 



k 2 (s){h(s)-h(m 2 p )) + (m 2 p -s)k 2 (m 2 p ) 



dh 



ds 



and T p = T p (m 2 ) . The P-wave energy-dependent width is given by 



= r ^ f fc(g) 

p UK) 



(51) 



(52) 



where k(s) — \ \fs (3 n (s) and k(m 2 ) are pion momenta in the p rest frame. The function h(s) is defined as 



Ma)= 2*w In yi + 2*w 



7T y^S 



2m7r 



(53) 



with dh/ds\ m 2 = h(m 2 ) [(8k 2 (m 2 )) 1 — (2m p ) x ] + (2nm 2 ) 1 . Since interference with the isospin-violating 
electromagnetic w — » 7r + 7r~ decay occurs only in e + e~ annihilation <5 is fixed to zero when fitting r data. The 
normalization BW^j(O) = 1 fixes the parameter d = f(0)/(T p m p ), which is found to be (Gounaris and 
Sakurai, 1968) 



3 m 
d = — 



2 m p + 2k(m ) 
In — + 



m„ 



m\m p 



tt k 2 (m 2 p ) 2m n 



2-Kk(m 2 ) nk 3 (m 2 ) 



(54) 



1. Fit to r data 

The fit of the GS model to the isovector r data establishes the need for the p(1450) contribution to the 
weak pion form factor (ALEPH Coll., 1997b, 2005). Only weak evidence is found for a p(1700) contribution 
lying close to the r end-point. Most of the fit parameters exhibit large correlations, which have to be taken 
into account when interpreting the results. The p mass uncertainty is found to be dominated by systematic 
effects, the largest being the knowledge of the n° energy scale (calibration). 



2. Combined fit to r and e + e data 

The ALEPH collaboration has performed a combined fit using r (ALEPH and CLEO) and e + e~ data (not 
yet including KLOE and SND) in order to better constrain the lesser known parameters in the phenomeno- 
logical form factor (ALEPH Coll., 1997b, 2005). For this study, the r spectral function is duly corrected 
for the isospin-breaking effects identified in Section V.A. Also photon vacuum polarization contributions are 
removed in the e + e~ spectral function since they are absent in the r data. In this way, the mass and width 
of the dominant p(770) resonance in the two isospin states can be determined. For the subleading amplitudes 
from the higher vector mesons, isospin symmetry is assumed so that common masses and widths can be used 
in the fit. 

The result of the combined fit is given in Table IV. The differences between the masses and widths of the 
charged and neutral p(770)'s are found to be (ALEPH Coll., 2005) 

m p - - m p0 = (2.4 ± 0.8) MeV , (55) 
F p --r p0 = (0.2 ± 1.0) MeV . (56) 

The mass splitting is somewhat larger than the theoretical prediction (< 0.7 MeV) (Bijnens and Gosdzinsky, 

1996) , but only at the 2a level. The expected width splitting, from known isospin breaking, but not taking 
into account any p mass splitting, is (0.7 ±0.3) MeV (Cirigliano et at, 2001; Davier et al., 2003a). However, 
if the mass difference is taken as an experimental fact, a larger width difference would be expected. From 
the chiral model of the p resonance (Cirigliano et al., 2001; Gomez Dumm et at, 2000; Guerrero and Pich, 

1997) , one expects a total width difference of (2.1 ± 0.5) MeV, which is only marginally consistent with the 
observedvalue. 

As seen in the previous applications, the situation should improve with the inclusion of the new SND 
data (SND-2 Coll., 2005). Also the observed p mass difference (55) is expected to decrease in this case. 
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FIG. 12 Phenomenological fit to the ALEPH it ty° spectral lunction (left) and to e + e annihilation data (right) using 
the Gounaris-Sakurai parameterization. The fits are taken from (ALEPH Coll., 2005). 



TABLE IV Combined fit to the pion form factor-squared to ALEPH, CLEO r and all e + e~ data, where vacuum 
polarization and known sources of isospin violation have been excluded from the latter data. The parameterization of 
the p(770), /o(1450) and p(1700) line shapes follows the Gounaris-Sakurai formula. Separate masses and widths are fit 
for the p(770), while common values are kept for the higher vector mesons. All mass and width values are in units of 
MeV and the phases are in degrees. The value of the \ 2 estimator per degree of freedom (DF) is also quoted. 



Simultaneous fit to r and e + e spectral functions (x 2 /DF = 383/326) 



m p-(770) 
«V(770) 
TO p(1450) 
TO p(1700) 



775.5 ± 0.6 
773.1 ± 0.5 
1409 ± 12 
1740 ± 20 



-p-(770) 
"p°(770) 
^p(1450) 
^p(1700) 



148.2 ± 0.8 
148.0 ± 0.9 
501 ± 37 
= 235 



5 (2.03 ± 0.10) x 10" 3 
(j> 5 (13.0 ±2.3) 

(3 0.166 ± 0.005 

4>f3 177.8 ± 5.2 

7 0.071 ± 0.006 

aL =0 



VI. TAU DECAYS AND HADRONIC VACUUM POLARIZATION 

Hadronic vacuum polarization in the photon propagator plays an important role in precision tests of the SM. 
This is the case for the evaluation of the electromagnetic coupling at the Z mass scale, a(M|), which receives 
a contribution Aahad(Af§) of the order of 2.8 x 10~ 2 . This correction must be known to a relative accuracy 
of better than 1% so that it does not limit the accuracy on the indirect determination of the Higgs boson 
mass from the measurement of sin 2 #iy Another example is provided by the anomalous magnetic moment 
°p = (ffp, — 2)/2 of the muon, where the lowest-order hadronic vacuum polarization component a h ad > LO j s the 
leading contributor to the uncertainty of the theoretical prediction. 

Starting from (Bouchiat and Michel, 1961; Cabibbo and Gatto, 1960, 1961) there is a long history of 
calculating the contributions from hadronic vacuum polarization in these processes. As they cannot be 
obtained from first principles because of the low energy scale involved, the computation relies on analyticity 
and unitarity so that the relevant integrals can be expressed in terms of an experimentally determined spectral 
function, which is proportional to the cross section for e + e~ annihilation into hadrons. The accuracy of the 
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calculations has therefore followed the progress in the quality of the corresponding data (Eidelman and 
Jegerlehner, 1995). Because the data were not always suitable, it was deemed necessary to resort to other 
sources of information. One such possibility was the use of the vector spectral functions derived from the study 
of hadronic r decays for the energy range less than m T (Alemany et al, 1998). Also, it was demonstrated 
that essentially perturbative QCD could be applied to energy scales as low as 1-2 GcV (ALEPH Coll., 1998c; 
OPAL Coll., 1999b), thus offering a way to replace poor e + e~ data in some energy regions by a reliable 
and precise theoretical prescription (Davier and Hockcr, 1998a, b; Erler, 1999; Groote et al, 1998; Kiihn and 
Steinhauscr, 1998; Martin and Zeppcnfcld, 1995). Finally, without any further theoretical assumption, it 
was proposed to use QCD sum rules (Davier and Hocker, 1998b; Groote et al, 1998) in order to improve 
the evaluation in energy regions dominated by resonances where one has to rely on experimental data. The 
complete theoretical prediction includes in addition quantum electrodynamics (QED), weak and higher order 
hadronic contributions. 

Since in this review of hadronic r decays we are mostly dealing with the low energy region 6 , common 
to both r and e + e~ data, and because of the current interest in the muon magnetic moment prompted by 
the excellence of the recent experimental results (see below), the emphasis will be on a^ ad ' LO rather than 
Aahad(-Mf)- We note that the presently achieved accuracy on Aahad(-^f) i s meeting the goals for the 
LEP/SLC/FNAL global electroweak fit. However the situation will change in the long run when very precise 
determinations of sin 2 #w, as could be available from the beam polarization asymmetry at the International 
Linear Collider, require a significant increase in the accuracy of Aah a d(Af§) (Abe et al, 2001). 



A. Muon magnetic anomaly 

One of the great successes of the Dirac equation (Dirac, 1928) was its prediction that the magnetic dipolc 
moment, [I, of a spin |s| = 1/2 particle such as the electron (or muon) is given by 

fli=9t7, — £ = e,/j,... (57) 

with gyromagnetic ratio gi — 2, a value already implied by early atomic spectroscopy. Later it was realized 
that a relativistic quantum field theory such as QED can give rise via quantum fluctuations to a shift in gg 

at = — ^— , (58) 

called the magnetic anomaly. In a now classic QED calculation, Schwinger (Schwinger, 1948) found the 
leading (one loop) effect (Fig. 13, lower-left) a e = a/(2n) ~ 0.00116, with a = e 2 /(47r) ~ 1/137.036, which 
agreed beautifully with experiment (Nagle et al, 1947a, b), thereby providing strong confidence in the validity 
of perturbative QED. Today, the tradition of testing QED and its SU(3)c x SU(2)l x U(1)y SM extension 
is continued (which includes strong and electroweak interactions) by measuring ct^ xp for the electron and 
muon even more precisely and comparing with a| M expectations, calculated to much higher order in per- 
turbation theory. Such comparisons test the validity of the SM and probe for new physics effects, which if 
present in quantum loop fluctuations should cause disagreement at some level. An experiment underway at 
Harvard (Gabriclsc and Tan, 1994) aims to improve the best present measurement (Van Dyck et at, 1987) 
of a e by about a factor of 15. Combined with a much improved independent determination of a, it would 
significantly test the validity of perturbative QED. It should be noted, however, that a e is in general not 
very sensitive to new physics at a high mass scale A because its effect on a e is expected to be quadratic in 
1/A (Czarnecki and Marciano, 2001) 



Aa e (A) ~°{-&) ( 59 ) 

and, hence, highly suppressed by the smallncss of the electron mass. It would be much more sensitive if Aa e 
were linear in 1/A; but that is unlikely if chiral symmetry is present in the m e — > limit. 



It is the property of the dispersion relation describing the hadronic contribution to (g — 2) M that the nir spectral function 
provides the major part of the total hadronic vacuum polarization contribution. Hence the experimental effort focuses on this 
channel. 
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FIG. 13 Representative diagrams contributing to a M . First column: lowest-order diagram (upper) and first order QED 
correction (lower); second column: lowest-order hadronic contribution (upper) and hadronic light-by- light scattering 
(lower); third column: weak interaction diagrams; last column: possible contributions from lowest-order supersymme- 
try. 



The muon magnetic anomaly has recently been measured for positive and negative muons with a relative 
precision of 5 x 10~ 7 by the E821 collaboration at Brookhaven National Laboratory (Muon (g — 2) Coll., 
2004). Combined with the older, less precise results from CERN (Bailey et ai, 1977), and averaging over 
charges, gives 

a c * p = (11 659 208.0 ±5.8) x 10~ 10 . (60) 

Although the accuracy is 200 times worse than Og Xp , a M is about m^/ml ~ 40, 000 times more sensitive to new 
physics and hence a better place (by about a factor of 200) to search for a deviation from the SM expectation. 
Of course, strong and electroweak contributions to are also enhanced by mi/m^ relative to et e ; so, they 
must be evaluated much more precisely in any meaningful comparison of a^ M with Eq. (60). Fortunately, the 
recent experimental progress in a® xp has stimulated much theoretical improvement of a^ M , uncovering errors 
and inspiring new computational approaches along the way, among these the use of hadronic r decays. 

It is convenient to separate the SM prediction for the anomalous magnetic moment of the muon into its 
different contributions, 

SM QED | weak . had la-\\ 

a M = <V + « M + a M , (61) 

where a^ ED = (11 658 472.0 ± 0.2) x 10~ 10 is the pure electromagnetic contribution (see (Czarnecki and 
Marciano, 1999; Hughes and Kinoshita, 1999) and references therein), 7 a™ ak = (15.4 ± 0.1 ± 0.2) x 10~ 10 , 
with the first error being the hadronic uncertainty and the second due to the Higgs mass range, accounts for 
corrections due to exchange of the weakly interacting bosons up to two loops (Czarnecki et ai, 2003) (see 
third column in Fig. 13). The term a^ ad can be further decomposed into 

had had.LO had. HO had,LBL ff{ r )\ 

where aJj ad ' LO is the lowest-order contribution from hadronic vacuum polarization and fl h ad ' HO the correspond- 
ing higher-order part (Section VI. E). At the 3-loop level in a, the so-called hadronic light-by- light (LBL) 



7 An improved calculation including all mass-dependent a 4 QED contributions has been published recently with a slightly differ- 
ent result (Kinoshita and Nio, 2004): a^ ED = 116 584 719.58(0.02)(1.15)(0.85) X 10" 11 . Here, 0.02 and 1.15 are uncertainties 
in the a 4 and a 5 terms, and 0.85 is from the uncertainty in a measured by atom interferometry. 
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scattering contributions, a!j ' LBL , must be estimated in a model-dependent approach. Those estimates have 
been plagued by errors, which now seem to be sorted out (Section VI. E). Nevertheless, the remaining overall 
uncertainties from hadronic vacuum polarization and LBL scattering in ajj ad represent the main theoretical 



error m a 



SM 



Owing to unitarity and to the analyticity of the vacuum-polarization function 8 , the lowest-order hadronic 
vacuum polarization contribution to can be computed via the dispersion integral (Gourdin and de Rafael, 
1969) 



a 



had.LO 



oo 

" 5 



K(s) 



(63) 



where K(s) is the QED kernel (Brodsky and Rafael, 1968), 

„2 X / x 



K{s) 



ln(l + x) - x + — 



x 2 ln.x 



l + x 



(64) 



with x = (1 - (3^/(1 + and = (1 - ^/s) 1 / 2 . In Eq. (63), R^{s) denotes the ratio of the "bare" 
cross section for e+e - annihilation into hadrons to the pointlike muon-pair cross section at center-of-mass 
energy \fs. The bare cross section is defined as the measured cross section corrected for initial-state radiation, 
electron- vertex loop contributions and vacuum-polarization effects in the photon propagator. As the only QED 
effect, photon radiation in the final state is to be included. The reason for using the bare (i.e., lowest-order) 
cross section is that a full treatment of higher orders is anyhow needed at the level of a^, so that the use of 
the "dressed" cross section would entail the risk of double-counting some of the higher-order contributions. 

The function K(s) ~ 1/s in Eq. (63) gives a strong weight to the low-energy part of the integral. About 
91% of the total contribution to a h ad > LO j g accum ulated at center-of-mass energies y/s below 1.8 GeV and 73% 



of a 



had.LO 



is covered by the itit final state, which is dominated by the p(770) resonance. 



The photon vacuum polarization function is the spin-one component of the time-ordered product of two electromagnetic 
currents, n 7 (q 2 ) (cf. Eq. (38)). Due to Ward identities, all contributions from self energy and vertex correction graphs to 
n 7 (q 2 ) cancel, and only vacuum polarization modifies the charge of an elementary particle. One can therefore write for the 
running electron charge (charge screening) at energy scale s 

a(s) 



1 - Aa(s) 

where q(0) is the fine structure constant in the long- wavelength Thomson limit and Aa(s) = Acq op (s) + ActhadM = 
— 47raRe[n 7 (s) — LT 7 (0)]. The leptonic part Aoi ep (s) is calculable within QED and known up to three loops (Stcinhauser, 
1998). However, quark loops are modified by long-distance hadronic physics that cannot be calculated within QCD. Instead, 
the optical theorem 

„ , . <r(e+e~ — » hadrons) , % 
127rlmn 7 (s) = \ = R(s) , 

and, owing to the analyticity of IT 7 , the dispersion relation 



n 7 (s)-n 7 (o) = - / ds'^ 

IT J S'(s' 



Imn 7 (s') 



■ s) — is 



lead to 



a(0)s f , , R(s') 
Aa had (s) = ^-Re / ds' - 



3-7T / s'(s' — s) — ie ' 

o 

Its r.h.s. can be obtained with the use of experimental data for R(s). Similarly, one obtains Eq. (63) for the hadronic vacuum 
polarization contribution to a M . 
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B. Input data and their treatment 

The information used for the evaluation of the integral (63) comes mainly from direct measurements of 
the cross sections in e + e~ annihilation and via CVC from r spectral functions. In addition to the data with 
the two- and four-pion final states shown in Section V.B, other final states and at different energies are also 
used (References for experimental data can be found in (Davier et al., 2003a)). However, due to the higher 
hadron multiplicity at energies above ~ 2.5 GeV, the exclusive measurement of the many hadronic final states 
is not practicable. Consequently, the experiments at the high-energy colliders ADONE, SPEAR, DORIS, 
PETRA, PEP, VEPP-4, CESR and BEPC have measured the total inclusive cross section ratio R. In some 
cases measurements are incomplete, as for example e + e~ — ► KKnn or e + e~ — ► 6tt, and one has to rely on 
isospin symmetry to estimate or bound the unmeasured cross sections (Davier et al., 2003a). To obtain the 
corresponding e + e~ annihilation cross sections using the r data, the isospin rotations (41, 47) are applied, 
and all identified sources of isospin breaking are corrected. 

To exploit the maximum information from available data, weighted measurements of different experiments 
at a given energy should be combined instead of calculating the integrals for every experiment and finally 
averaging them. Different averaging methods arc used in the literature varying from classical \ 2 minimization 
methods (Alemany et al., 1998; Davier et al., 2003a; Menke, 2001), to the so-called clustering method (Hagi- 
wara et al., 2004) in which a cluster of arbitrary size is predefined and precise data are used to rescale imprecise 
measurements within the (correlated) systematic errors of the imprecise data set. 

To overcome the lack of precise data at threshold energies, a phenomenological expansion in powers of s is 
used: 



This interpolating function satisfies the known form factor constraints at s = 0. Above the 2-7T threshold, where 
a branch cut develops in the form factor, such a parameterization has been derived within the framework of 
Chiral Perturbation Theory (Colangelo et al., 1996). Exploiting precise results from space- like data (NA7 
Coll., 1986), the pion charge radius-squared is constrained to 



and the two parameters c\.2 are fit to the data in the range [2m w , 0.6 GeV]. Good agreement is observed 
in the low energy region and this parameterization is used for evaluating the integrals in the range up to 
0.5 GeV (Davier et al, 2003a). 

C. QCD for the high energy contributions 

Since problems still remain regarding the spectral functions from low-energy data, the prediction from 
essentially perturbative QCD is used above an energy of 1.8 and 5 GeV, respectively. The details of the 
calculation can be found in (Davier and Hocker, 1998a, b) and in the references therein. 

The perturbative QCD prediction uses a next-to- next-to-leading order 0(a|) expansion of the Adler D- 
function (Gorishnii et al., 1991; Surguladze and Samuel, 1991), with second-order quark mass corrections 
included (Chetyrkin et al, 1996). R(s) is obtained by evaluating numerically a contour integral in the 
complex s plane (c/. the detailed discussion of this technique in Section VII). Nonpcrturbative effects are 
considered through the Operator Product Expansion, giving small power corrections controlled by gluon and 
quark condensates. The value a s (M§) = 0.1193 ± 0.0026, used for the evaluation of the perturbative part, 
is taken as the average of the results from the analyses of r decays (ALEPH Coll., 1998c) (see Section VII) 
and of the Z width in the global electroweak fit (LEP EWWG, 2002). Uncertainties on the QCD prediction 
are taken to be equal to the common error on a s (M§), to half of the quark mass corrections and to the full 
nonperturbative contributions. A local test of the QCD prediction can be performed in the energy range 
between 1.8 and 3.7 GeV. The contribution to a h ^ d < LO in this region is computed to be (33.87±0.46) x lO" 10 
using QCD, to be compared with the result, (34.9 ± 1.8) x 10~ 10 from the data. The two values agree within 
the 5% accuracy of the measurements. 

The evaluation of a h ad > LO was shown to be improved by applying QCD sum rules (Davier and Hocker, 
1998b; Groote et al., 1998). The latest full reanalyses (Davier et al., 2003a, b) however did not consider this 
possibility since the main problem at energies below 2 GeV was the inconsistency between the e + e~ and r 




| = l + -(r 2 ).s + Cl .s 2 + c 2 ,s 3 . 
6 



(65) 



(r 2 ) = (0.439 ± 0.008) fm 2 



(66) 



33 




FIG. 14 Fit of the pion form factor from 4mJ to 0.35 GeV 2 using a third-order expansion with the constraints at s = 
and the measured pion r.m.s. charge radius from space-like data (NA7 Coll., 1986). The result of the fit is integrated 
only up to 0.25 GeV 2 . The e+e - data do not yet include the recent SND measurement (SND-2 Coll., 2005). 



input data, which had to be resolved with priority. Also the improvement provided by the use of QCD sum 
rules results from a balance between the experimental accuracy of the data and the theoretical uncertainties. 
With an increase in the experimental precision, the gain from the theoretical constraint would be smaller than 
before. 



D. Results for the lowest order hadronic vacuum polarization 

Many calculations of the hadronic vacuum polarization contribution have been carried out in the past, 
taking advantage of the e + e~ data available at that time and benefiting from a more complete treatment of 
isospin-breaking corrections (Cirigliano et al., 2001, 2002). In this review, we will follow the latest published 
analysis by (Davier et al., 2003b) (and the update (Hocker, 2004) including data from KLOE, see also the 
result in (Troconiz and Yndurain, 2004)), which considers input from both r and e + e™ data. 

Figure 15 gives a panoramic view of the e + e~ data in the relevant energy range. The shaded band below 
2 GeV represents the sum of the exclusive channels. The QCD prediction is indicated by the cross-hatched 
band. Note that the QCD band is plotted taking into account the thresholds for open flavor B states, in order 
to facilitate the comparison with the data in the continuum. However, for the evaluation of the integral, the 
bb threshold is taken at twice the pole mass of the b quark, so that the contribution includes the narrow T 
resonances, according to global quark-hadron duality. 

Compiling all contributions, the e + e~-based result (including the KLOE data, but not yet SND) for the 
lowest-order hadronic contribution is (Hocker, 2004) 

a had,LO = ( 693 4 ± 5 3 ± 3 5rad ) x 1Q -10 ^ (g 7 ) 

where the second error is due to the treatment of (potentially) missing radiative corrections in the older 
data (Davier et al., 2003a). The new KLOE data decreased the contribution by about 3 x 10~ 10 . For a 
comparison we also give the results without KLOE for the e + e~ and r-based analyses (Davier et al., 2003b): 
fl had,LO[ e+e -_ based ] = (- 696 3 ± g 2exp ± 3 . 6rad ) x 10 -io and a had,LO [ T _based] = (711.0 ± 5.0 oxp ± 0.8 rad ± 

2.8 SU(2) ) x 10- 10 . 
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FIG. 15 Compilation of the data contributing to aJj ad,L °. Shown is the total hadronic over muonic cross section ratio 
R. The shaded band below 2 GeV represents the sum of the exclusively measured channels, with the exception of the 
contributions from the narrow resonances that are given as dashed lines. All data points shown correspond to inclusive 
measurements. The cross-hatched band gives the prediction from (essentially) perturbative QCD (see text). 



E. Hadronic three-loop effects 

The three-loop hadronic contributions to a^ M involve one hadronic vacuum polarization insertion with 
an additional loop (either photonic or another leptonic or hadronic vacuum polarization). They can be 
evaluated (Krause, 1997) with the use of the same e + e~ — > hadrons data sets used for a^ ad,L . Denoting that 
subset of C(a/7r) 3 hadronic contributions a h ad . NLO j we quote here the result of a recent analysis (Hagiwara 
et al., 2004), 

a had,NLO = ( _ 9 _g ± Q1) x 1Q -10 ^ (g8) 

which is consistent with earlier studies (Alemany et al., 1998; Krause, 1997). It would change by about 
—0.3 x 10~ 10 if the r data were also used. 

More controversial arc the hadronic light-by-light scattering contributions illustrated in the lower diagram of 
the second column in Fig. 13. Since it invokes a four-point correlation function, a dispersion relation approach 
using data is not possible and a first-principles calculation (e.g., lattice gauge theory (Blum, 2003)) has so far 
not been carried out. Instead, calculations involving pole insertions, short distance quark loops (Bijnens et 
al., 1996; Hayakawa et al., 1996) and charged-pion loops have been individually performed in a large Nc QCD 
approach. The pseudoscalar poles (7r°,ry and r]') dominate such a calculation. Unfortunately, in early studies 
the sign of the contribution was incorrect. Its correction (Bijnens et al, 2002; Blokland et al, 2002; Hayakawa 
and Kinoshita, 1998; Knecht and Nyffeler, 2002; Knecht et al, 2002) led to a large shift in the o„ M prediction. 
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A representative estimate (Davier et al, 2003a) of the LBL contribution, which includes 7r, 77 and 77' poles as 
well as other resonances, charged pion loops and quark loops, currently gives a^ad.LBL ^ ^g_g _|_ 3 ^ x j^q-io 
(representative) . 

A new analysis of LBL scattering (Melnikov and Vainshtein, 2004) takes into account the proper matching 
of the asymptotic short distance behavior of pseudoscalar and axial-vector contributions with the free quark 
loop behavior. It yields the somewhat larger result a h ad < LBL — (13.6 ± 2.5) x 10~ 10 . 

As pointed out in (Davier and Marciano, 2004), several small but (likely) negative contributions such as 
charged pion loops and scalar resonances were not included in (Melnikov and Vainshtein, 2004) and could 
reduce somewhat the magnitude of the result. The authors (Melnikov and Vainshtein, 2004) provide a con- 
sistency check on their result much in the spirit of the EW hadronic triangle diagram study (Czarnecki et al., 
2003; Vainshtein, 2003) discussed in (Davier and Marciano, 2004). Using constituent quark masses in the LBL 
diagram combined with a pion pole contribution that properly accounts for the chiral properties of massless 
QCD and avoiding short distance double counting, they find a ^ ad > LBL ~ 12.0 x 10~ 10 with about half of the 
contribution coming from quark diagrams and the other half from the pion pole. Using that result along with 
an inflated error assigned so that the representative result and the one obtained in (Melnikov and Vainshtein, 
2004) overlap within their errors, (Davier and Marciano, 2004) suggest the value 

a hadiBL _ (12 ± 3 5) x 1Q - 10 (6Q) 

Further resolution of the LBL scattering contribution is very important. 
Combining Eqs. (68) and (69), one finds 

a had,3-loop = 2 ± 3 5 ) x 1Q -10 ^ ( 7Q ) 

which we identify with a|f d < HO in Eq. (62). 



F. Comparing a M between theory and experiment 

Summing the results from the previous sections on a^ EU , a^ w , aJf d ' LO , a ^ ad < NLO , and a h ad > LBL 7 one obtains 
the SM prediction for a^. The newest e + e~-based result reads (Hocker, 2004) 

al M = (11 659 182.8 ±6.3 had XO+NLO ± 3.5had,LBL ±0.3 Q ed+ew) x HT 10 . (71) 

This value can be compared to the present measurement (60); adding all errors in quadrature, the difference 
between experiment and theory is 

< P - al M = (25.2 ± 9.2) x 10~ 10 , (72) 

which corresponds to 2.7 "standard deviations" (to be interpreted with care due to the dominance of exper- 
imental and theoretical systematic errors in the SM prediction). A graphical comparison of the result (71) 
with previous evaluations (also those containing t data) and the experimental value is given in Fig. 16. 

Whereas the evaluation based on the e + e~ data only disagrees with the measurement, the evaluation 
including the tau data is consistent with it. The dominant contribution to the discrepancy between the two 
evaluations stems from the 7T7T channel with a difference of (— 11.9±6.4 oxp ±2.4 ra d±2.6 S u(2) (±7.3 to tai)) x 10~ 10 , 
and a more significant energy-dependent deviation 9 . As a consequence, during the previous evaluations of 
fl had,LO^ -(- ne resu its using respectively the r and e + e~ data were quoted individually, but on the same footing 
since the e + e _ -based evaluation was dominated by the data from a single experiment (CMD-2). 

The seeming confirmation of the e + e _ data by KLOE could lead to the conclusion that the r-based result 
be discredited for the use in the dispersion integral (Hocker, 2004). However, the newest SND data (SND-2 
Coll., 2005) alter this picture in favor of the r data, along with prompting doubts on the validity of the 
KLOE results (see discussion in Section V.C). Comparing the SND and CMD-2 data in the overlapping 



The systematic problem between r and e+e~ data is more noticeable when comparing the t~ — > tt~it v t branching frac- 
tion with the prediction obtained from integrating the corresponding isospin-brcaking-corrected e+e - spectral function (c/. 
Section V). 
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FIG. 16 Comparison of the result (71) (Hocker, 2004) labelled DEHZ 04 with the BNL measurement (Muon (g - 2) 
Coll., 2004). Also given are the previous estimate (Davier et al, 2003b), where the triangle with the dotted error 
bar indicates the r-based result, as well as the estimates from (Hagiwara et al, 2004; Jegerlehner, 2003; Troconiz and 
Yndurain, 2004), not yet including the KLOE data. 

energy region between 0.61 GeV and 0.96 GeV, the SND-based evaluation of a h ad > NLO j s f ounc l to be larger by 
(9.1 ± 6.3) x 10~ 10 . However, once these two experiments are averaged using the trapezoidal rule, the increase 
over the CMD-2 only evaluation is reduced to 2 x 10~ 10 due to the much smaller systematic errors and the 
higher density of the CMD-2 measurements. 

The proper response to the problems in the LO hadronic contribution is to further improve the experimental 
precision by a renewed experimental effort, and to continue to improve theory. In this respect, new e + e~ — > 
hadrons data and further study of LBL scattering could potentially reduce the overall theoretical uncertainty. 
The excitement that was caused by the deviation (72) stems from the expectation that new physics could 
cause an effect of similar magnitude, some examples of which are discussed in (Davier and Marciano, 2004; 
Marciano, 1996). 

G. Conclusions and perspectives 

After considerable effort, the experiment E821 at Brookhaven has improved the determination of by 
about a factor of 14 relative to the classic CERN results of the 1970's (Bailey et al., 1977; Farley and Picasso, 
1990). The result is still statistics limited and could be improved by another factor of 2 or so (to a precision 
of 3.0 x 10~ 10 ) before systematic effects become a limitation. Pushing the experiment to that level seems to 
be an obvious goal for the near term. In the longer term, a new experiment with improved muon acceptance 
and magnetic fields could potentially reach 0.6 x 10~ 10 (Hertzog, 2005; Roberts, 1999). 

In spite of the new and precise data on the two-pion spectral function from the KLOE and SND collabora- 
tions, the lowest-order hadronic vacuum-polarization contribution remains the most critical component in the 
SM prediction of a M . The publication of the KLOE results led to an outward confirmation of the discrepancy 
between the r data and e + e _ annihilation observed in the it + tt~ channel (Davier et al, 2003b). Including 
the KLOE data in the hadronic evaluation, the e + e~-based SM prediction of a M was found to differ from the 
experimental value by 2.7 "standard deviations". 

An opposite point view is advocated in (Maltman, 2005). A QCD sum-rule analysis of low-energy e + e~ 
data results in the value a s (M|) = 0.1150±0.0024, which is somewhat lower than the world average (Particle 
Data Group, 2004), 0.1187± 0.0020, and even less consistent with the average 0.1200±0.0020, which excludes 
r data (for the sake of comparison) and heavy quarkonium results (the accuracy of which has been questioned 
in (Brambilla et al., 2004)). The same analysis using the r vector spectral function provides agreement, 
with a value of a s (M§) = 0.1201 ± 0.0021. Although the disagreement between the e + e~ data and other 
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measurements is not yet significant, the observed trend is worth more investigations. It is supported by the 
newest measurement from SND that tend to amend the compatibility with the r data in the 2 pion channel. 

The further improvement of the e + e~ data is the better route and one is looking forward to the forthcoming 
results on the low- and high-energy two-pion spectral function from the CMD-2 collaboration. These data will 
help to significantly reduce the systematic uncertainty due to the corrective treatment of radiative effects, often 
omitted by part by the previous experiments. The initial-state-radiation program of the BABAR collaboration 
has already proved its performance by publishing the spectral functions for ir + ir~ir° (BABAR Coll., 2004) and 
2ir + 2ir~ (BABAR Coll., 2005a), while results for the two-pion final state are expected. All sources considered, 
and assuming all discrepancies to be resolved, it appears however difficult to reach an uncertainty much better 
than 3 x 10~ 10 in this sector 10 . 

The 30% uncertainty in the hadronic LBL contribution (3.5 x 10~ 10 ) is largely model-dependent. Here 
one could conceive that further work following the approach of (Melnikov and Vainshtein, 2004) or lattice 
calculations could help to reduce the error. Doing much better will be difficult, but an improvement in the 
precision by about a factor of 2 would well match short-term experimental capabilities. 

If the above experimental and theoretical improvements do occur, they will lead to a total uncertainty in 
Aa M = a° xp — a^ M of about 5.0 x 10~ 10 , i.e., a factor of 2 reduction with respect to today. This would provide 
an important step towards a better test of the SM. The result could be used in conjunction with future collider 
discoveries to sort out the properties of new physics (e.g., the size of tan/3 in SUSY) or to further constrain 
possible extensions of the SM. 

VII. HADRONIC TAU DECAYS AND QCD 

Tests of Quantum Chromodynamics and the precise measurement of the strong coupling constant a s at 
the t mass scale (Braaten, 1989; Braaten et ai, 1992; Le Diberder and Pich, 1992b; Narison and Pich, 1988), 
carried out for the first time by the ALEPH (ALEPH Coll., 1993) and CLEO (CLEO Coll., 1995b) collab- 
orations have triggered many theoretical developments. They concern primarily the perturbative expansion 
for which different optimized rules have been suggested. Among these are contour-improved (resummed) 
fixed-order perturbation theory (Le Diberder and Pich, 1992a; Pivovarov, 1991, 1992), effective charge and 
minimal sensitivity schemes (Dhar, 1983; Dhar and Gupta, 1984; Grunberg, 1980, 1984; Stevenson, 1981), the 
large-/?o expansion (Altarclli et ai, 1995; Ball et al, 1995; Neubert, 1996), as well as combinations of these 
approaches. They mainly distinguish themselves in how they deal with the fact that the perturbative series 
is truncated at an order where the missing part is not expected to be small. Also, in the discussion of the 
so-called contour improved approach we will point out the importance to fully retain the complete information 
from the renormalization group when computing integrals of perturbatively expanded quantities. With the 
publication of the full vector and axial-vector spectral functions by ALEPH (ALEPH Coll., 1997b, 1998c) 
and OPAL (OPAL Coll., 1999b) it became possible to directly study the nonperturbative properties of QCD 
through vector minus axial-vector sum rules. These analyses are described in Section IX. 

One could wonder how r decays may at all allow us to learn something about perturbative QCD. The 
hadronic decay of the r is dominated by resonant single particle final states. The corresponding QCD interac- 
tions that bind the quarks and gluons into these hadrons necessarily involve small momentum transfers, which 
are outside the domain of perturbation theory. On the other hand, perturbative QCD describes interactions 
of quarks and gluons with large momentum transfer. Indeed, it is the inclusive character of the sum of all 
hadronic r decays that allows us to probe fundamental short distance physics. Inclusive observables like the 
total hadronic r decay rate R T (20) can be accurately predicted as function of a s (m^.) using perturbative 
QCD, and including small nonperturbative contributions within the framework of the Operator Product Ex- 
pansion (OPE) (Shifman et ai, 1979). In effect, R T is a doubly inclusive observable since it is the result of 
an integration over all hadronic final states at a given invariant mass and further over all masses between 
and m T . Hadronic r decays are more inclusive than e + e~ annihilation data, since vector and axial-vector 
contributions are of approximately equal size, while the isospin-zero component in e + e _ data is three times 
smaller than the isovector part, which is (about) equal to the vector contribution in r decays. The scale m T 



We also point out that lattice gauge theories with dynamical fermions can in principle provide a determination of a M a ' (Blum, 
2003). 
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lies in a compromise region where a s (m 2 ) is large enough so that R T is sensitive to its value, yet still small 
enough so that the perturbative expansion converges safely and nonperturbative power terms are small. 

If strong and electroweak radiative corrections are neglected, the theoretical parton level prediction for 
SUc(Nc), N c = 3 reads 

Rr = N c (\V ud \ 2 + \V us \ 2 )=3 , (73) 

so that we can estimate a perturbative correction to this value of approximately 21% to obtain (20), assuming 
other sources to be small. One realizes the increase in sensitivity to a s compared to the Z hadronic width, 
where due to the three times smaller a s (M§ ) the perturbative QCD correction reaches only about 4%. 

The nonstrange inclusive observable R T can be theoretically separated into contributions from specific quark 
currents, namely vector (V) and axial-vector (A) ud and us quark currents. It is therefore appropriate to 
decompose 

Rr = Rr.V + Rr.A + Rr.S , (74) 

where for the strange hadronic width i? T .s vector and axial- vector contributions are not separated so far due 
to the lack of the corresponding experimental information for the Cabibbo-suppressed modes. Parton- level 
and perturbative terms do not distinguish vector and axial- vector currents (for massless partons). Thus the 
corresponding predictions become R t .v/a = {Nc/2)\V u d\ 2 and R T: s = Nc\V us \ 2 , which add up to Eq. (73). 

A crucial issue of the QCD analysis at the r mass scale concerns the reliability of the theoretical description, 
i.e., the use of the OPE to organize the perturbative and nonperturbative expansions, and the control of 
unknown higher-order terms in these series. A reasonable stability test is to continuously vary m T to lower 
values y/s~o < m T for both theoretical prediction and measurement, which is possible since the shape of the full 
r spectral function is available. The kinematic factor that takes into account the r phase space suppression 
at masses near to m T is correspondingly modified so that ^/sq represents the new mass of the r. 



A. Renormalization group equations 

Like in QED, the subtraction of divergences in QCD is equivalent to the renormalization of the coupling 
(a s = g 2 /4:Tr), the quark masses (m q ), etc, and the fields in the bare (superscript B) Lagrangian such as 
= s £ Z as a s , m q — s £ Z m m q , etc. Here s is the renormalization scale, e the dimensional regularization 
parameter, and Z denotes a series of renormalization constants obtained from the generating functional of 
the bare Green's function. The renormalization procedure introduces an energy scale, s, which represents the 
point at which the subtraction to remove the divergences is actually performed. This leads to the differential 
Renormalization Group Equations (RGE) 



da s 
dins 

1 dm q 
m q dins 



(3(a s ) = -a 2 s J2Pn< , (75) 

n 

7(a s ) = -a s J2lna" , (76) 



with a s — a s /n. Expressed in the modified minimal subtraction renormalization scheme (MS) (Bardeen et al, 
1978; t'Hooft et at, 1972) and for three active quark flavors at the r mass scale, the perturbative coefficients, 
known to four loops (van Ritbergen et al., 1997a, b), are 11 

Po = 2.25 , Pi = 4 , p 2 « 10.05989 , fa w 47.22804, (77) 



The full expressions for an arbitrary number of quark flavors (n/) are (van Ritbergen et al., 1997a,b) 
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70 = 1 , li= 3.79166 , 72 « 12.42018 , 73 « 44.26278 . (78) 

The solutions for the evolutions (75) and (76) are obtained from integration 12 (Chetyrkin et al, 1997b; 
Stevenson, 1981) 

a s (s) oc 

A 2 J s \P(a' s ) /3 (a' s ) 2 +/3i(4) 3 / J A)«) 2 + /3i(a' s ) 3 

o»(s) 
o«(«) 

w 9 (s) _ /" , 7(4) 



In^H = / da',-£-% , (80) 

Os(so) 

where the MS asymptotic scale parameter /l and the RG-invariant quark mass rh q (which appears after 
evaluating the integral (80)) are fixed by referring to known values a s (s ) and m q (s ), respectively, at scale 
s Q . 

For most practical purposes, it is however unnecessary to explicitly perform the integrations (79, 80), since 
the evolution of the differential RGEs can be solved numerically by means of single-step integration using a 
modified Euler method. It consists of Taylor-developing, say, Eq. (75) around the reference scale sq in powers 
of r} = ln(s/so), and reordering the perturbative series in powers of a s = a s (sn), 

a s (s) = a s - poVa 2 s + {-piV + Pffi) a 3 s + (-fori + \fofcr? - flfo 3 ) ^ 

+ (~fov + \pW + 3/?o/V ~ §P 2 oPiV 3 + /?oV) a 5 s (81) 

+ (-far) + \fafor? + \fofor? ~ ~ Wofc 7 ! 3 + ^/?o/W " /?oV) a* + 0(a 7 s ) . 

For later use, we have expressed the expansion up to order a^, which involves the unknown five- loop coefficient 
/?4. An efficient integration of the RGEs can also be obtained using the CERNLIB routine "RKSTP", which 
performs a 4 th order Rungc-Kutta single step approximation with excellent accuracy. 

B. Theoretical prediction of R T 

According to Eq. (39) the absorptive parts of the vector and axial-vector two-point correlation functions 
ni"2 v/a( s )' Wlt ^ * ne s P m J 01 th° hadronic system, are proportional to the r hadronic spectral functions 
with corresponding quantum numbers. The nonstrange ratio R T y+A can be written as an integral of these 
spectral functions over the invariant mass-squared s of the final state hadrons (Braaten et al, 1992) 



so 2 

RtV+a(s ) = 12itS ew f— (l- — ] (l + 2—)lmn^(s + ie) +lmU^(s + ie) , (82) 
J s \ s J [\ s J 



where n( J ) can be decomposed as n( J ) = | V^.^ | 2 [J^JJy + n^^j . The lower integration limit is zero because 
the pion pole is at zero mass in the chiral limit. 
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12 Note that the integration is only denned for a non- vanishing derivative /3(a 3 ) jt- 
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Vanishing f3(a s ) may indeed occur for 
example in effective charge schemes with large K4 or K$ coefficients (see later in this section). 
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FIG. 17 Integration contour for the r.h.s. in Eq. (83). 



The correlation function n^ J ^ is analytic in the complex s plane everywhere except on the positive real 
axis where singularities exist. Hence by Cauchy's theorem, the imaginary part of n( J ) is proportional to the 
discontinuity across the positive real axis 



so 



Sp 

dsw(s)ImU(s) = -— I dsw(s)U(s) , (83) 
2m J 

\s\=s 



where w(s) is an arbitrary analytic function, and the contour integral runs counter-clockwise around the circle 
from s = so + ie to s = s — ie as indicated in Fig. 17. 

The energy scale sq = m 2 is large enough that contributions from nonperturbative effects are expected 
to be subdominant and the use of the OPE is appropriate. The kinematic factor (1 — s/s ) 2 suppresses the 
contribution from the region near the positive real axis where (s) has a branch cut and the OPE validity 
is restricted due to large possible quark-hadron duality violations (Braaten, 1988; Poggio et al, 1976). 

The theoretical prediction of the vector and axial-vector ratio R t ,v/a can hence be written as 

Rry,A = \\V u tfS^[l + 5^+5^ + C m v %+ Y.C]v/ A ) > (84) 

\ D=4,6,... / 

with the massless perturbative contribution 5^°\ the residual non-logarithmic electroweak correction <5^ w = 
0.0010 (Braaten and Li, 1990) (c/. the discussion on radiative corrections in Section V.A), and the dimension 
D = 2 perturbative contribution 5^ji™ q ) A from quark masses, which is lower than 0.1% for u, d quarks. The 
term denotes the OPE contributions of mass dimension D 

S ud'v/A= 2^ C udy/A {s,n) , (85) 



where the scale parameter // separates the long-distance nonperturbative effects, absorbed into the vacuum 
expectation elements (O u d(n)}, from the short-distance effects that are included in the Wilson coefficients 
C ud , v/A (s,») (Wilson, 1969). Note that 6$ V+A = (5 { u % + 6$ A )/2. 



1. The Perturbative Prediction 



The perturbative prediction used by the experiments follows the work of (Le Diberder and Pich, 1992a). 
The perturbative contribution is given in the chiral limit. Effects from quark masses have been calculated 
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in (Chetyrkin and Kwiatkowski, 1993) and are found to be well below 1% for the light quarks. As a conse- 
quence, the contributions from vector and axial-vector currents coincide to any given order of perturbation 
theory and the results are flavor independent. 

For the evaluation of the perturbative series, it is convenient to introduce the analytic Adler function (Adler, 
1974) 

OO 



where the second identity is the dispersion relation (c/. Section VI. A). The derivative of the correlator avoids 
extra subtractions (renormalization) on the r.h.s. of Eq. (86), which are unrelated to QCD dynamics. The 
function D(s) calculated in perturbative QCD within the MS renormalization scheme depends on a non- 
physical parameter /i occurring as ln(^ 2 /s). Furthermore it is a function of a s . On the other hand, since 
D(s) is connected to a physical quantity, the spectral function ImiT(s), it cannot depend on the subjective 
choice of fi. This can be achieved if a s becomes a function of /i providing independence of D(s) of the choice 
of /i. Nevertheless, in the realistic case of a truncated series, some [i dependence remains and represents an 
irreducible systematic uncertainty. D(s) is then a function of fi 2 /s, a s and can be understood as an effective 
charge (see below) that obeys the RGE. Choosing = s, the perturbative series reads 

oo 

£>(a s ( S ))=5> n <(«). (87) 

n=0 

with renormalization scheme dependent coefficients r n . 

To introduce the Adler function in Eqs. (82) and (83) one uses partial integration 



j> dsg(s)U(s) = - j> ^ (G(s) - G(s )) s 



<m(s) 



ds 

\s\ = s \s\=so 



(88) 



with the kernel g(s) = (1 — s/s ) 2 (l + 2s/s )/s and G(s) = ds' g(s'). This gives (Le Diberder and Pich, 
1992a) 



1 + *(°> = -27TZ J> - 



s 



i - 2 - +2 (-T-(- 4 

SO \SoJ \So 



D(s) . (89) 



The perturbative expansion of the Adler function can be inferred from the third-order calculation of the 
e + e~ inclusive cross section ratio R e + e -(s) — a(e + e~ — > hadrons (7))/cr(e + e~ — > (7)) (Celmaster and 

Gonsalves, 1980; Chetyrkin et al, 1979; Dine and Sapirstein, 1979; Gorishnii et al, 1991; Le Diberder and 
Pich, 1992a; Surguladze and Samuel, 1991) 

00 

n=0 

with the K n (0 functions up to order n — 5 (Le Diberder and Pich, 1992a) 

K (0 - K , ^i(0 = K X , £ 2 (0 =K 2 + L, 

K z {£,)=K z + {b 1 + 2K 2 )L + L 2 

K 4 (0 = K 4 + (b 2 + 2hK 2 + 3 K 3 ) L + (^h+SK^j L 2 + L 3 , (91) 
K 5 (0 = K 5 + (63 + 2 b 2 K 2 + 3 hK 3 + AK 4 )L+ (h 2 1+ 3b 2 + 7 b x K 2 + 6K 3 )L 
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where bi — Pi/Po and L — /3oln£. The factor £ = sq/s in Eq. (90) represents the arbitrary renormalization 
scale ambiguity of which the Adler function (86) is independent. The £ dependence of the Ki(£) coefficients 
is determined by inserting the expansion (81) into (90), and equating the two series obtained for £ = 1 and 
arbitrary £ at each order in a s . 13 The coefficients K n are known up to third order in a§. For n > 2 they 
depend on the renormalization scheme used, while the first two coefficients are universal 



K n = 1 



K x = 1 



K 2 (MS) = F 3 (MS) , 



K 3 (MS) = F 4 (MS) + -tt 2 /? 2 , 



(92) 



where F 3 (MS) = 1.9857-0.1153/1/ and F 4 (MS) = -6.6368 - 1.2001 n f - 0.0052 n 2 f are the coefficients of the 
perturbative series of R e + e -. For n/ = 3 one has Ki = 1.640 and K$ = 6.371. With the series (90), inserted 
in the r.h.s of Eq. (89), one obtains the perturbative expansion 



(93) 



n=l 



with the functions 



1 f ds 



A (n) (a s ) = — & — 



2iri J -; 
|«|=»o 



1-2— + 2 ( — 

so \ so 



^- J dip [1 + 2e iv - 2e i3v - a^(£s a e llp ) , 



(94) 



where s = — s e 4V has been substituted in the second line and the integral path proceeds according to Fig. 17. 



2. Fixed-order perturbation theory (FOPT) 

Inserting the series (81) in Eq. (93), evaluating the contour integral, and collecting the terms with equal 
powers in a s leads to the familiar expression (Le Diberder and Pich, 1992a) 

* (0) =£ [*n(fl+ft.(0] (^f^T , (95) 

where the g n are functions of K m<n and /3 m < n -i 5 and of elementary integrals with logarithms of power 
m < n in the integrand. Setting £ = 1 and replacing all known and Ki coefficients by their numerical 
values, Eq. (95) simplifies to 

= a s (s ) + (1.6398 + 3.5625) a 2 s {s ) + (6.371 + 19.995) a 3 (s ) 
+ (ICi + 78.003) a^so) 
+ {K 5 + 14.250 Ki - 391.54) a 5 s {s ) 

+ (K 6 + 17.813 if 5 + 45.112 K 4 + 1.58333 /3 4 - 8062.1) a^(s ) , (96) 

where for the purpose of later studies we have kept terms up to sixth order. When only two numbers are 
given in the parentheses, the first number corresponds to K n , and the second to g n . 



13 Note that this procedure corresponds to a fixed-order perturbation theory rule since known higher order pieces of the expan- 
sion (81) are neglected when solving 

N N N 

n— n — n — 

for a truncated series with maximum known order N. A resummation of the subleading logarithms has been proposed in (Dams 
and Kleiss, 2004). 
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FIG. 18 Comparison between the iteratively evolved asim^-e^) (exact solution) and the fourth-order Taylor expan- 
sion (81) on the circle ip — —ty . . . ir (Menke, 1998). This example uses a s (mj) = 0.35. 

The FOPT series is truncated at given order despite the fact that parts of the higher coefficients g n >4(0 
are known to all orders and could be resummed. These known parts are the higher (up to infinite) order terms 
of the expansion (81) that are functions of (3 n <3 and K n < 3 only. In effect, beyond the use of the perturbative 
expansion of the Adler function (90), two approximations have been used to obtain the FOPT series (96): (i) 
the RGE (75) has been Taylor-expanded and terms higher than the given FOPT order have been truncated, 
and (ii) this Taylor expansion is used to predict a s (—s) on the entire \s\ = so contour. 

3. Contour-improved fixed-order perturbation theory (CIPT) 

A more promising approach to the solution of the contour integral (94) is to perform a direct numerical 
evaluation by means of single-step integration and using the solution of the RGE to four loops as input for the 
running a s (— £s) at each integration step (Le Diberder and Pich, 1992a; Pivovarov, 1991, 1992). It implicitly 
provides a partial resummation of the (known) higher order logarithmic integrals and improves the convergence 
of the perturbative series. While for instance the third order term in the expansion (96) contributes with 17% 
to the total (truncated) perturbative prediction, the corresponding term of the numerical solution amounts 
only to 6% (assuming a s {m^.) = 0.35). This numerical solution of Eq. (93) will be referred to as contour- 
improved fixed-order perturbation theory (CIPT) in the following. 

Single-step integration also avoids the Taylor approximation of the RGE on the entire contour, since a s is 
iteratively computed from the previous step using the full known RGE. 14 Figure 18 shows the comparison 
between the iteratively evolved a s (m^.e tv ) (referred to as the exact solution), and the fourth-order Taylor 
expansion (81) on the circle ip — —it . . .it (Menke, 1998). The differences appear to be significant far away 
from the reference value at (tp = 0). However, as shown below, the numerical effect of this discrepancy on the 
contour integral is small compared to the main difference between FOPT and CIPT, which is the truncation 
of the perturbative series after integration. 



It has been tested that the single-step integration, i.e., as^soe 1 ^^^^) is computed from a s (gsoe 11 *'), using a modified Euler 
method (that is, a Taylor expansion of a s in each integration step), or using explicit 4 th order Runge-Kutta integration, leads 
to the same result, as far as the step size is chosen to be small enough (1000 integration steps on the contour are found to be 
largely sufficient). 
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4. Effective charge perturbation theory (ECPT) 

Inspired by the pioneering work in (Dhar, 1983; Dhar and Gupta, 1984; Grunberg, 1980, 1984; Stevenson, 
1981) the use of effective charges to approach the perturbative prediction of R T has become quite popular in 
the recent past (Korner et al, 2001a; Maxwell and Tonge, 1996, 1998; Raczka, 1998). The advocated advantage 
of this technique is that the perturbative prediction of the effective charge is renormalization scheme (RS) 
and scale invariant since it is a physical observable. The effective t charge is defined by 

a T = <5 (0) , (97) 

and obeys the RGE 

, CO 

= /3 T (a r ) = -«?E^« a " ■ ( 98 ) 

n=0 

The first two terms in the expansion of [3 T are universal {i.e., RS-indcpcndent), while the higher order terms 
must be computed using perturbation theory. The prescription for this computation is obtained from the 
relation between the renormalization scheme and scale used to obtain the FOPT coefficients in the series (95), 
and the requirement of RS and scale invariance. Since a T is scheme and scale invariant, its evolution coefficients 
(3 T ,i have the same property. "Another way of looking at this is to say that the method of effective charges 
involves a specific choice of renormalization scheme and scale (...), so that the energy evolution of the 
observable is identical to the beta-function evolution of the coupling" (Dinsdale and Maxwell, 2004). 
The function (3 T is obtained from the equation 

da T d[3 T (a T ) 
da s d/3(a s ) 
which, setting £ — 1, can be rewritten in the form 



(99) 



, . da T da T da s 
T T dins da. dins 



da T 

= -i—PM 



(100) 

=a s (a T ) 



Inserting Eq. (95) for a T on the l.h.s. and the r.h.s. of Eq. (100), and equating the coefficients for each order 
in a s leads to the coefficients 15 

0r,2 =02- PlC2 + Po(-C 2 2 + Cs) , (101) 

0t,3 = /33-2/3 2C2 +/3 lC 2 + 2/3o(2c3-3c 2 c 3 + C 4) , (102) 

0r,4 = (34-3f3 3 C2 + f32(^4-c 3 )+f3 1 (~2c 2 c 3 + c 4 ) 

+ A)(-14 c\ + 28 c\c 3 - 5 c§- 12 c 2 c 4 + 3 c 5 ) , (103) 

where Cj = Ki + gi(l) (cf. Eq. 95) and [3 Tt o = (3q, (3 t .\ = [3\ because K = K\ = 1. Inserting the known fa 
and a coefficients, one finds (Korner et at, 2001a) 

/3 T ,2 = -12.320, 

p r 3 = -182.72 + 4.5 K A , 

(3 T i = -236.2 + f3i- 40.27 K 4 + 6.75 K b , (104) 

/3 T)5 = 16427 - 20.80 /3 4 + /%- 521.6 K 4 -65.79 K 5 + 9 K % , 

with a noticeable sign alteration for /3 Tj2 (and also for /3 T:3 if K 4 < 40). Neglecting nonperturbative contribu- 
tions, one has 

a T (s ) = i? T ,y +j4 (s )/(3 \V ud \ 2 ) - 1 , (105) 
which using Eq. (22) and the \V u d\ value given in Section IV yields a r (m^) = 0.2219 ± 0.0055. 



For completeness we also give the sixth-order coefficient, which is later needed for extrapolation tests 

Pt,5 = 05 ~ 4 /?4C 2 + 2 /3 3 (4 4 - c 3 ) + 4 /3 2 (-2 4 + c 2 c 3 ) + 0i (-6 4 + 16 4c 3 - 3 c§ - 8 c 2 c 4 + 2 c 5 ) 
+ 4/3 (12c| - 30 4 C 3 + 12 c 2 c| + 14 c|c 4 -4c 3 c 4 - 5c 2 c 5 + c 6 ) . 
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5. Estimating unknown higher order perturbative coefficients 

Compared with the MS RGE (75) the (3 T (a T ) series is significantly less convergent, taking into account that 
a T (m 2 ) ~ 1.8 x a s (m 2 ). This property has been used in the past to estimate an effective K% s coefficient, 
which was found to be Kf 1 <~ 27 (Kataev and Starshcnko, 1995). It approximates the true value of K4 in the 
limit of vanishing higher order contributions. In effect, the derivative of a T (m^) with respect to K4 is large 16 



da T (m%) 



~ 0.0020 , (106) 

a T (m2)=0.22, if 4 =25 



so that the precise experimental error on a T (m 2 ) allows us to predict K4 with an accuracy of a few units, if 
we assume that the subleading contributions can be neglected. However, this assumption is unfounded since, 
taking naively the unknown higher order coefficients to grow like a geometric series, i.e., z n <~ z^_ 1 /z„_2, 
with z = (3,K, and using the estimate K4 = 25 so that K 5 ~ 98, /3 T) 4 <~ —360 (with P4 ~ 222), one finds 



da T (m^.) 



~ 0.0004 , (107) 

a T (m*),K i ,K B ,0 i ... 



which is not small with respect to (106), in particular considering the coarseness of the assumptions that 
went into the estimate (107). It is found that even the sixth-order coefficient is not negligible. Assuming the 
same geometric growth, one would expect K§ ~ 395 with a derivative of da r (m~) / dKg ~ 5 x 10~ 5 , which 
appears small. However, considering that the uncertainty on Kq is at least of the same size as its estimated 
value (which is probably still optimistic) generates an uncertainty in a T (m^) that is equivalent to a{Ki) ~ 10. 
Taking correspondingly the fifth-order uncertainty of the order of the estimated value of K 5 , we conclude 
that the systematic uncertainty on K4 obtained with the effective charge scheme is at least of the order of 
a{K A ) ~ 23. 

Another method to estimate if| ff has been suggested in (Le Diberder, 1995). It consists of increasing the 
sensitivity of the CIPT prediction of 5^ on K4 by reducing the renormalization scale £ in Eq. (93). The result 
27±5 is consistent with the one obtained from ECPT. However, these methods are not independent as can be 
seen in the following. Using the integral (79), and taking s — > 00 so that a T , a s — > 0, we can write (Stevenson, 
1981) 

ln T2 - ln i = 21ni T = - + °(°«) = C2 + °(°») ' ( 108 ) 

A 2 A 2 . A l3 a s (s) Paa T (s) 

A T ~ Ae C2/f3 ° w 3.2/1 . (109) 

Since the asymptotic scale parameters are reciprocal to the scale, transforming FOPT into ECPT is equivalent 
to reducing the renormalization scale, at lowest order. 

In terms of the measurement of a s (m 2 ) the three perturbative methods FOPT, CIPT and ECPT use similar 
information. However their dependence on the unknown higher order coefficients K n >4 can be dramatically 
different. Figure 19 shows the evolution of a r (s ) for FOPT, CIPT with renormalization scales £ = 1 and 
£ = 0.4, and ECPT, all normalized to a T (m 2 ) = 0.22. The solid lines correspond to fifth order perturbation 
theory with assumed K4 = 25, K$ = 98, and for variations of these (dashed and dotted lines). One observes 
relatively small dependence for FOPT and CIPT with £ = 1, while large and very large effects are found with 
CIPT (£ = 0.4) and ECPT. The latter effect provides the sensitivity to estimate the unknown coefficients. To 
better analyze this sensitivity, we have compared the AK4 variation in fifth-order ECPT to the fourth-order 
effect (dashed-dotted line in the lower right hand plot of Fig. 19). Including the bold estimate for the next 
order significantly reduces the dependence on the K4 term. This cancellation becomes obvious when looking 
at the numerical expansion (104): the K4 coefficient enters /3 Tj 3 and /3 T .4 with different signs, and with a larger 
coefficient at fifth order to compensate for the a T suppression. A similar behavior is observed for K 5 when 
turning on the sixth-order term. We conclude that as for the K4 term, the strong K 5 dependence observed 



One can obtain the solution for a T (m^) as a function of the /3 T coefficients without explicitly integrating the RGE (98), by 
using as fixed point of the differential equation the one-loop solution at s — » 00, A 2 = se~ 1 /(' 3 o a -r(s)) > so that the four-loop 
a T (00 — * m 2 ) evolution reproduces the experimental value. This gives A = 0.41 GeV. 
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FIG. 19 Evolution ol a T (so) for the different perturbative methods, and for variations of the unknown higher order per- 
turbative coefficients K4 (dashed) and K5 (dotted). All methods are normalized to a T {rrfc) = 0.22. The inserted frames 
give the as(m?) values derived from this a T value in the respective methods. The error bars show the uncertainties, 
added in quadrature, from the K4 and K5 variations quoted on the plots. 



in ECPT is (by part) an artifact of the truncation of the perturbative series, and it is therefore dangerous to 
use it for an estimate of the unknown coefficients. 

The embedded frames in Fig. 19 give the a s (m^.) values derived from the given a T value in the respective 
methods. The error bars indicate the uncertainty, added in quadrature, from the K4 and K§ variations 
quoted on the figures. One observes that the difference in a s between FOPT and CIPT is not covered by 
the uncertainties on K± and K§. This is because the main contribution to the difference in these methods is 
due to the truncation of the known higher order terms of the Taylor expansion after the contour integration, 
rather than due to the unknown coefficients. Due to the good convergence of these series, the dependence on 
the K§ coefficient is marginal. 

It should be mentioned that a large computational effort is underway with the goal to calculate the K4 
coefficient. The large number of five-loop diagrams needed to calculate the two-point current correlator at 
this order is somewhat discouraging, however the results on two gauge invariant subsets are already available. 
The subset of order O(afn^) was evaluated long ago through the summation of renormalon chains (Beneke, 
1993), while the much harder subset 0(a^nj) was recently calculated (Baikov et ai, 2002). The known parts 
amount to about half of the expected K4 coefficient assuming a geometric growth. Using the results of (Kataev 
and Starshenko, 1995) for different n/ obtained with the effective charge method, it is possible to compare 
the nj estimated term at the a\ order of the Adler function expansion (3.64ny), with the newly computed 
one (1.875n^). The difference provides a conservative estimate for the uncertainty of the ECPT procedure to 
predict the full K4 term, as one could expect the sum of all a\n\ terms to be closer to the true value than the 
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term-by-term comparison. A successful test of the procedure has been performed for K3, where all the a^n^ 
terms are exactly known. The two unknown remaining terms 0{ct^nj) and 0(afn°), the direct calculation 
of which is prevented by the present computing power available to this project, have been estimated (Baikov 
et al, 2003) to be 

Ki = -0.01009n^ + 1.875 n) -31.8 n f + 105.7 , (110) 

where the first two terms are known exactly. The calculation for n/ = 3 yields K4 ~ 27 with an uncertainty 
estimated to ±16 in (Baikov et al, 2003). 

6. Infrared behavior of the effective charge 

In the effective charge scheme, the physical coupling should have a finite infrared behavior and hence 
must not suffer from the Landau-pole that lets MS (but also momentum subtraction schemes) break down 
below a scale \fs <~ 1 GeV. The universal two-loop function (3 T (a T ) does however suffer from a singularity at 
y/s ~ 1.1 GeV (for a T (m 2 ) = 0.22). Only adding the negative three-loop coefficient regularizes fi T {a T ) and 
leads to a freezing a T (s — > 0) ~ 0.62. Adding higher orders does not qualitatively alter this behavior, however 
the freezing value is very unstable against variations of the unknown coefficients (which again exhibits the 
mediocre convergence of the ECPT series). The lower right hand plot of Fig. 19 shows the evolution of a T (s) 
for different orders and assumptions upon the unknown higher order coefficients (cf. the discussion in (Brodsky 
et al, 2002)). 

7. Renormalons and the large-/3o expansion 
Perturbative series in quantum field theory 17 

N 

R N = ]T r n a n+1 , (111) 

n=0 

be they abelian or not, are divergent at N — > 00 for any value of a. Because the number of Feynman graphs 
increases at each order n, one expects a factorial growth of the perturbative coefficients of the corresponding 
S'-matrix series 

r n ~Ka n nW , (112) 

with constant K, a, 7. The divergent behavior of perturbative expansions often has physical significance: it 
indicates the nonperturbative structure of the vacuum and its excitations. To be a useful approximation of 
the true value R, the series (111) must be asymptotic to R for given a, i.e., it exists a number Km so that 18 
\R — Rn\ < K N+ ia N+2 . If r n <~ K n , the best approximation of R is achieved when the series is truncated 
at its minimal (absolute) term and the truncation error can be roughly estimated by the size of the minimal 
term. The accuracy in the point of minimal sensitivity depends on the strength of the perturbative expansion 
parameter (the coupling). The larger it is, the lower the order when minimal sensitivity is reached, and the 
worse is the approximation. Note that as long as the computed order is below the order of minimal sensitivity, 
there is no actual difference between a divergent or a convergent series. 

A convenient way to deal with this divergence is to consider the Borel transform 

00 

B ^](*) = E^* n . ( 113 ) 

n=0 



The suppression of the constant term in the definition of Eq. (Ill) is convenient to avoid 1/a terms in the definition of the 
Borel integral (115). It is without consequence for the following discussion. 

Our discussion of renormalons has greatly benefited from the comprehensive review given in (Beneke, 1999). 
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FIG. 20 Multi-fermion loop insertion (renormalons) into a fermion anti-fermion vacuum polarization diagram. 



which is believed to have a finite radius of convergence in the i-plane (Mueller, 1993). The n th fixed order 
perturbation coefficient is then generated by the n th derivative 



d n B[R](t) 



dt n 



(114) 

t=o 



The explicit factor n! in the denominator of B[R](t) makes the Borel-transformed series much better behaved. 
Summing up all orders leads to the integral representation 19 

oo 

R(a) = J dte- t/a B[R](t) . (115) 
o 

If the integral exists, it gives the Borel sum of the original divergent series. Existence requires that B[R](t) has 
no singularities in the integration range. However, for a factorial behavior akin to Eq. (112) (non-alternating 
series), the Borel transform is given by B[R] (t) oc (1 — at) ~ 1-7 . Hence, if 7 is positive there exists a singularity 
t = l/a so that the integral (115) does not exist. Indeed, "the divergent behavior of the original series is 
encoded in the singularities of its Borel transform!' (Beneke, 1999). One notices that the larger a, i.e., the 
faster diverging the series is, the closer to zero moves the singularity. 

Let us consider a renormalizcd fermion-loop insertion into a quark-antiquark vacuum polarization diagram 
(c/. one bubble in Fig. 20). Assuming dominance of the (/?o a s(~ s )) n term (denoted as the large-f3o expansion), 
one can resum this chain with a large number of n bubbles. If one neglcts higher order terms of the (3- function, 
it is thus possible to derive estimates for the FOPT coefficients of a given perturbative series at all orders. This 
assumption is supported empirically by the observation that the [3q term dominates second order radiative 
corrections for many observables in the MS scheme (Ball et ai, 1995). (We shall however anticipate here 
that the accuracy of these estimates is insufficient for the purpose of a precise perturbative prediction of R T , 
and difficult to control.) The procedure provides a naive non-abelianization of the theory since lowest order 
radiative corrections do not include gluon self-coupling. 

Following this line, the Adler function can be expanded as (Beneke and Braun, 1995) 

N 

D(s) = l + a s Y / < (dnffi + 5 n ) , (116) 

71=0 

where the coefficients d n are computed in terms of fermion bubble diagrams 20 (see Fig. 20). The MS Borel 
transform of Eq. (116) is given by (Beneke, 1993; Broadhurst, 1993) 

„ r „ n/ . 327T 6,. u ^ (-l) fe fc 

B[D]{U) = — 63 u i-(i- U ) 2 E (fc2 _ (1 _ u)2)2 - ( 117 ) 



The Borel integral is easily derived by using the integral representation n! = J dv cxp(— v)v n , and exchanging summation 
and integration in Eq. (Ill) (see, e.g., (Lovett-Turner and Maxwell, 1994)). Note however that this interchange is formally 
not justified in presence of the Landau pole for the coupling constant. 

One identifies the coefficients d n in Eq. (116) with their leading-n/ pieces in the expression d n = d^n'j + . . . + d^ . It 
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FIG. 21 Ultraviolet (UV) and infrared (IR) renormalons in the Borel plane of the Adler function. 



where u ~ — /?ot, and where we have chosen the renormalization scale to be equal to the momentum scale. 
Neglecting the corrections S n , the series (116) leads to the large-/3o expansion of D(s). The first elements of 
the series are 

d = 1 , di/3 = 0.496 , d 2 Pl = 15.71 , 

d^Pl = 24.83 , d 4 /3o = 787.8 , d 5 ffi = -1991 , (118) 

which compare only coarsely with the exact terms (92) where these are known. The Borel transform B[D](u) 
has singularities on the real axis, which can be distinguished in infrared (IR renormalons) for small virtuality, 
and ultraviolet singularities (UV renormalons) for high virtuality of the exchanged gluon. One finds an infinite 
sequence of IR (UV) renormalons at positive (negative) integer u, with the exception of u = 1 (see Fig. 21). 

The integral (115) may still be solvable in presence of renormalons by acquiring an imaginary part from 
moving the integration contour above or below the real axis. The associated ambiguity is exponentially small 
in the expansion parameter a, and is therefore of nonperturbative origin. Correspondingly, the IR poles give 
rise to ambiguities when one uses the generators (114) to recompute D(s) from its Borel transform. Using in 
first order 

as{s) = am!m 2 ) ' (119) 

one obtains for the associated ambiguity for u = 2, 3, . . . 

/ A2\ u 

AD(s) ~ e - u /(/3oa s ) = _ . ( 12 o) 

These IR renormalons proportional to (A 2 / s) u are reabsorbed into the nonperturbative terms of the OPE. The 
absence of a u = 1 IR rcnormalon is thereby related to the impossibility to build a gauge invariant operator 
of dimension D = 2, which is an important property that is exploited for the determination of a s from R T . 
The lowest IR renormalon is hence of the order (A 2 /s) 2 . On the contrary to the UV renormalons, this term 
is scale independent and cannot be decreased (Beneke and Zakharov, 1992). It should therefore have physical 
meaning and can be identified with the gluon condensate (0|G Mi ,G A " y |0) as the leading infrared contribution to 
the Adler function. However, while IR renormalons lead to nonperturbative operators for power corrections, 
it does not necessarily contain all of them. For example, at D = 4 one misses the quark condensate (0|<7<7|0), 
which is generated by chiral symmetry breaking that does not occur in perturbation theory. 



k yln+2-k] 
n+2-k 

k\(n-k)\ ' 



has been shown in (Broadhurst, 1993) that the djj 1 ' in MS are given by (s = p?) 
with the generating function 

r l r>2 — m / J yn-2 



and P (x) ■ 



32 >r-^ (-l) k k 



3(1 + 2!) (k 2 -x 2 ) 2 
k=2 



50 



Due to asymptotic freedom, governed by the negative sign in front of the RGE /3-function, the UV renor- 
malons arise on the negative real axis 21 (on the contrary to QED where they occur on the positive side) so 
that they are outside the integration range of (115) and insofar harmless, that is, Borel-summable (Mueller, 
1993). Let us return for a moment to the series (111) in the large-/3o expansion. For large n, the factorial 
growth of the perturbation series is dominated by the contribution of the leading UV renormalon k = — 1 with 
alternating coefficients r n ~ n!n 7 (— /?o)™- A guess at which order iV m j n minimal sensitivity is achieved can be 
obtained from the argument that the series is convergent, i.e., reliable at order n+ 1 if a(r n+ i/r n < 1). Con- 
sidering leading IR and UV renormalons only, one finds from the asymptotic behavior r n+ i/r n ~ npo/2 (IR) 
and r n+ i/r n ~ — n(3 (UV). Hence the breakdown of convergence is first caused by the UV renormalon and 
-V m in ~ l/(/3o fl s) is of order one. A series truncated at finite order brings an intrinsic limitation of accuracy 
along with it. The associated error is reasonably estimated with the magnitude of the order N = 7V m ; n term 
of the perturbative series. That gives using for the leading UV renormalon 

A 2 

|rjv|af ~ N\N^(-/3 a s ) N x logarithms , (121) 

s 

where the r.h.s. is obtained with the use of the Stirling formula, and where we have fixed the renormalization 
scale at [i 2 = s. The truncation of the perturbative series is accompanied by an uncertainty which scales 
like 1/s, i.e., with apparent dimension D = 2. This contribution is however not to be confounded with 
nonperturbative IR renormalon ambiguities 22 . Similarly to Eq. (116) the R T FOPT series (95) can be written 

as 

N 

5<°> (s) = l + a s J2< KM + Sn) , (122) 

n=0 

where d^/3g + = K n+ \ + g n+ i(l). Inserting Eq. (117) into (89), and taking advantage of the factorized s 
dependence in the large-/3o expansion, the Borel transform of (122) is (Beneke, 1993) 

*(")■ (; + + ro) • c») 

The first elements of the series are 

do = 1 , d[P = 5.119 , d T 2 Pl = 28.78 , 

dlffi = 156.6 , dl$ = 900.9 , d T 5 f3% = 4867 . (124) 

They compare reasonably well with the exact terms (96) where these are known. In this approach, the 
unknown fixed-order fourth-order coefficient is estimated to be K4 ~ 79, which is significantly larger than 
the estimate based on ECPT, and also larger than the one obtained from the large-/3 expansion of the Adlcr 
function (118). 

Apart from the conceptual interest in the study of renormalons, the main question regarding the concrete 
use of the large-/3o expansion certainly is: how large are the uncertainties? From Eq. (124) we see that the 
series is positive and that minimal sensitivity is reached at n ~ 5. On the other hand, the resummation of 
the known g n (£) coefficients (96) shows that large negative coefficients occur at higher order FOPT so that 
one would need a far larger than a geometric growth of the K n coefficients to keep the series positive. A 
direct comparison between the large-/?o expansion coefficients for the Adler function (118) with the exact 
calculation (92) reveals an oscillating behavior of the large-/3o series, i.e., small d n coefficients for n-even and 
large d n for n-odd. For example, the large-/3o dg, coefficient is substantially overestimated. One also has to 
worry about important contributions from the ds term, which is found to be 788. As a consequence, our 
conclusion is that the uncertainties associated with the large-/?o expansion are not under sufficient control 
to improve the perturbative prediction of R T . It is also an unreliable estimator of the uncertainty of the 
perturbative series. 



In QCD, the presence of IR renormalons implies that nonperturbative corrections should be added to define the theory 
unambiguously. The same is true for UV renormalons in QED. This corresponds to the observation of large coupling constants 
in the respective energy regimes. 

It has been shown in (Bcnckc and Zakharov, 1992) that the truncation uncertainty at ./V m j n scales actually as A 2 s/fj, 4 once the 
renormalization scheme dependence of A is taken into account. 
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TABLE V Massless perturbative contribution to R T {m1) for the various methods considered, and at orders n > 1 
with a s (m 2 T ) — 0.35. The value of K4 is set to 25, while all unknown higher order K n> 4 and f3 n >'j. coefficients are 
assumed to follow a geometric growth. 



6° 



Pert. Method 


n = 1 


n = 2 


n = 3 


(n = 4) 


(n = 5) 


(n = 6) 


T 4 




FOPT (£ = 1) 


0.1114 


0.0646 


0.0365 


0.0159 


0.0010 


-0.0086 


0.2283 


0.2208 


CIPT (Taylor RGE, £ = 1) 


0.1573 


0.0317 


0.0126 


0.0042 


0.0011 


0.0001 


0.2058 


0.2070 


CIPT (full RGE, £ = 1) 


0.1524 


0.0311 


0.0129 


0.0046 


0.0013 


0.0002 


0.2009 


0.2025 


CIPT (full RGE, £ = 0.4) 


0.2166 


-0.0133 


0.0006 


-0.0007 


0.0010 


-0.0007 


0.2032 


0.2048 


ECPT 


0.1442 


0.2187 


-0.1195 


-0.0344 


-0.0160 


-0.0120 


0.2090 


0.1810 


Large-/3o expansion 


0.1114 


0.0635 


0.0398 


0.0241 


0.0155 


0.0093 


0.2388 


0.2636 



8. Comparison of the perturbative methods 

Some comparisons between the fixed-order PT methods have already been given for the discussion of the 
estimate of the unknown higher order coefficients in Section VII. B. 5. The slow convergence of ECPT (and 
CIPT at small renormalization scales, ( < 1) has been pointed out in this respect. To further study the 
convergence of the perturbative series, we give in Table V the contributions of the different orders in PT to 
r)(°) for the various approaches and using a s (m^.) — 0.35. We assumed here K t = 25 and a geometric growth 
of all other unknown PT and RGE coefficients. In the case of CIPT the results are given for the various 
techniques used to evolve a s (soe tlf ): the truncated Taylor expansion (81), Runge-Kutta integration of the 
RGE with £ = 1, and Runge-Kutta integration with £ = 0.4. 

As advocated in (Le Diberder and Pich, 1992a), faster convergence is observed for CIPT compared to FOPT 
yielding a significantly smaller error associated with the renormalization scale ambiguity (while, somewhat 
counter intuitive, the uncertainty due to the unknown K4, and higher order coefficients is similar in both 
approaches). Our coarse extrapolation of the higher order coefficients could indicate that minimal sensitivity 
is reached at n ~ 5 for FOPT, while the series further converges for CIPT. Although the Taylor expansion in 
the CIPT integral exhibits significant deviations from the exact solution on the integration circle (c/. Fig. 18), 
the actual numerical effect from this on r)(°) is small (cf. second and third column in Table V). As discussed in 
Section VII. B. 4, the convergence of the ECPT series is much worse than for FOPT and CIPT. Consequently, 
the difference between truncation at n = 4 and n — 6 is expected to be significant. A similar behavior is 
observed for the large-/?o expansion. 

The CIPT series is found to be better behaved than FOPT (as well as ECPT) and is therefore to be preferred 
for the numerical analysis of the r hadronic width. As a matter of fact, the difference in the result observed 
when using a Taylor expansion and truncating the perturbative series after integrating along the contour 
(FOPT) with the exact result at given order (CIPT) is exhibiting the incompleteness of the perturbative 
series. However, it is even worse than that since large known coefficients are neglected in FOPT so that the 
difference between CIPT and FOPT is actually overstating the perturbative truncation uncertainty. This can 
be verified by studying the behavior of this difference for the various orders in perturbation theory given in 
Table V. The CIPT- vs. -FOPT discrepancy increases with the addition of each order, up to order four where 
a maximum is reached. Adding the fifth order does not reduce the effect, and only beyond fifth order the two 
evaluations become asymptotic to each other. As a consequence varying the unknown higher order coefficients 
and using the difference between FOPT and CIPT as indicator of the theoretical uncertainties overemphasizes 
the truncation effect. 
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9. Quark-mass and nonperturbative contributions 



Following SVZ (Shifman et ai, 1979), the first contribution to R T beyond the D = perturbative expansion 
is the non-dynamical quark mass correction of dimension 23 D = 2, i.e., corrections in powers of 1/sn- The 
leading D = 2 corrections induced by the light-quark masses are computed using the running quark masses 
evaluated at the two-loop level. Evaluation of the contour integral in FOPT leads to (Braaten et ai, 1992) 



c(2 — mass) 
°ud,V/A 



■Mi + y^ )) m " (so) + m ^ o) 



±4 1 



25 A m„(sn)m d (sn) 
3 I s 



(125) 



where mj(so) are the running quark masses evaluated at the scale s using the RGE 7-function (76). The 
following values arc used by ALEPH for the renormalization group invariant quark mass parameters rhi defined 
in (80) 



m u = (8.7 ± 1.5) MeV , m d = (15.4 ± 1.5) MeV , m s = (270 ± 30) MeV . 



(126) 



The dimension D — 4 operators have dynamical contributions from the gluon condensate (a s GG) and light 
u,d quark condensates {unqiqj), which are the matrix elements of the gluon field strength-squared and the 
scalar quark densities, respectively. Remaining D = 4 operators arc running quark masses to the fourth 
power. Inserting the Wilson coefficients of these operators (Becchi et ai, 1981; Broadhurst, 1981; Chetyrkin 
et al., 1985; Generalis, 1989) in the contour integral one obtains (Braaten et ai, 1992) 



J ud,V/A 



(a s GG) 



+ 16tt 2 



1 + |o2(*o) 



(K T m d ) (uu =F dd)) 



18^ 2 a 2 (s ) 
48 1 



(m u uu + niddd) 4 



(m k q k q k ) 



k—u,d.. 



22 

7 a s (s ) 7 



(m„(s ) T m d (so)) (^(s ) T ^( s o)) 



± 6 m M (s )m d (so) (m M (s ) T m d (s )) 2 + 36 to|(so)to|(so) 



(127) 



Where two signs are given the upper (lower) one is for V (A). The gluon condensate vanishes in first order 
a s (so). However, there appear second order terms in the Wilson coefficients due to the s dependence of a s , 
which after integration becomes a 2 s . 

The contributions from dimension D — 6 operators are complex. The most important operators arise from 
four-quark dynamical effects of the form q i Tiqjq k T2qi- Other operators, such as the triple gluon condensate 
whose Wilson coefficient vanishes at order a s , or those which are suppressed by powers of quark masses, are 
neglected in the evaluation of the contour integrals performed in (Braaten et ai, 1992). The large number of 
independent operators of the four-quark type occurring in the D = 6 term can be reduced by means of the 
factorization (or vacuum saturation) assumption (Shifman et ai, 1979) to leading order a s . The operators 
are then expressed as products of scale dependent two-quark condensates a s (fj,)(q i qi(fi))(qjqj(fj 1 )). To take 
into account possible deviations from the vacuum saturation assumption, one can introduce an effective scale 
independent operator pa s (qq) 2 that replaces the above product. The effective D = 6 term obtained in this 
way is (Braaten et ai, 1992) 



x(6) 

°ud,V/A 



7 

-11 



2567T 4 pa s (qq)^ 
sT~ 



(128) 



Corrections of dimension D = 2 refer to the mass dependence of the perturbative OPE coefficients of the unit operator. 
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predicting a different sign and hence a partial cancellation between vector and axial-vector contributions 2 

The dimension D = 8 contribution has a structure of nontrivial quark-quark, quark-gluon and four-gluon 
condensates whose explicit form is given in (Broadhurst, 1985). For the theoretical prediction of R T it is 
custom to absorb the complete long and short distance part into the scale invariant phcnomcnological D = 8 
operator (0$), which is fit simultaneously with a s and the other unknown nonperturbative operators to date. 

Higher order contributions from D > 10 operators are expected to be small as, equivalent to the gluon 
condensate, constant terms and terms in leading order a s vanish in Eq. (83) after integrating over the contour. 



C. Results 

It was shown in (Le Diberder and Pich, 1992b) that one can exploit the shape of the spectral functions to 
obtain additional constraints on a s (so) and — more importantly — on the nonperturbative effective operators. 
The t spectral moments at s = m 2 arc defined by 



n T.y/A 



ds 1 - -=■ 



mi 



mi 



dR r 



V/A 



ds 



(129) 



where R^ v / A = R T ,v/A- The factor (1 — s/m 2 ) fc suppresses the integrand at the crossing of the positive real 
axis where the validity of the OPE less certain (Braaten, 1988) and the experimental accuracy is statistically 
limited. Its counterpart (s/to 2 )^ projects upon higher energies. The spectral information is used to fit 
simultaneously a s (m 2 ) and the effective operators (a s GG), pa s (qq) 2 and (Od) f° r dimension D = 4, 6 and 
8, respectively. Due to the intrinsic experimental correlations only five moments are used as input to the fit. 

In analogy to R T , the contributions to the moments originating from perturbative and nonperturbative 
QCD are decomposed through the OPE 



^ M = |iK d | 2 5 EW (i+^) + ^ + e>7r H) + E 



kl) 

V/A 



(130) 



D=4,6,. 



The prediction of the perturbative contribution takes the form 

3 



(131) 



71=1 



with the functions 



2r(3 + fc )[ J?i^L + 2; r ^ 



T{A + k + i) r(5 + fc + ^) 



-/(— ,l+^3 + fc)-2/(— ,2 + £,3 + k 

\so ) \s 



(132) 



which make use of the elementary integrals 1(7, a, b) = J 7 t a_1 (l — t) b ~ 1 dt. The contour integrals are numer- 
ically solved for the running a s (— £s) (CIPT). 

In the chiral limit 25 and neglecting the small logarithmic s dependence of the Wilson coefficients, the 



It has been pointed out in (Chetyrkin and Kwiatkowski, 1993) that Eq. (128) assumes that the same p parameter can be used 
for both vector and axial-vector contributions. 

In the chiral limit, vector and axial-vector currents arc conserved so that slly^ = sll^ 1 = 0. 
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dimension D nonperturbative contributions in Expression (130) read 
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dimO=Z3 



(133) 



where the matrix is defined by the choice of the coefficients for the moments k = 1, £ = 0,1,2,3 and the 
corresponding dimension D. For completeness, we also give the coefficients for dimension D = 10, which is 
not considered by the experiments since they do not contribute to R®° v / A at leading order. With increasing 
weight I the contributions from low dimensional operators vanish. For example, the only nonperturbative 
contribution to the moment R^y/A stems from dimension D = 8 and beyond. Hence, in a fit using spectral 
moments, D = 8 will be determined by R]. 3 V / A - This observation is in some sense a "paradox", as higher 
moments project upon higher masses on the contrary to R T and the spirit of the OPE, where the higher 
dimension terms are enhanced at small so. 

For practical purpose it is more convenient to define moments that are normalized to the corresponding 
Rt.v/a in order to decouple the normalization from the shape of the r spectral functions 



n ki 

u t.v/a 



R 



kl 



,V/A 



R 



t,V/A 



(134) 



The two sets of experimentally almost uncorrelated observables — R t .v/a on one hand and the spectral mo- 
ments on the other hand — yield independent constraints on a s (m^.) and thus provide an important test of 
consistency. The correlation between these observables is completely negligible in the V + A case where 
Rt.v+a is calculated from the difference R T — R T ,s, which has no correlation with the hadronic invariant 
mass spectrum. One experimentally obtains the D^v/A by integrating the weighed normalized invariant 
mass-squared spectrum. The corresponding theoretical predictions are easily modified. 



1. The ALEPH determination of a s (m T ) and nonperturbative contributions 

Combined fits to experimental spectral moments and the extraction of a s (m^) together with the leading 
nonperturbative operators have been performed by ALEPH, CLEO and OPAL (ALEPH Coll., 1993, 1998c, 
2005; CLEO Coll., 1995b; OPAL Coll., 1999b) using similar strategies and inputs. Let us follow in this report 
the most recent analysis performed by the ALEPH collaboration (ALEPH Coll., 2005). 

Since the determination of a s should be model- independent, the experiments proceed by fitting simultane- 
ously the nonperturbative operators, which is possible since the correlations between these and a s turn out 
to be small enough. The theoretical framework and its intrinsic uncertainties were the subject of the previous 
section. Minor additional contributions originate from the CKM matrix element IK^I, the electroweak radia- 
tive correction factor S'ew? the light quark masses m u ^ and the quark condensates (Section VII. B. 9). The 
largest contributions have their origin in the truncation of the perturbative expansion. Although it introduces 
some double-counting for the systematic error, the procedure used in (ALEPH Coll., 2005) considers separate 
variations of the unknown higher order coefficient K4 and the renormalization scale. The renormalization 
scale is varied around m T from 1.1 to 2.5 GeV with the variation over half of the range taken as systematic 
uncertainty. Taking advantage of the new theoretical developments discussed above, the value K4 = 25 ± 25 
is used in (ALEPH Coll., 2005). 

The fit minimizes the \ 2 of the differences between measured and adjusted quantities contracted with the 
inverse of the sum of the experimental and theoretical covariance matrices. The results of (ALEPH Coll., 
2005) are listed in Table VI. The precision of a s (m 2 ) obtained with the two perturbative methods employed 
is comparable, however their central values differ by 0.01-0.02, which is larger than the theoretical errors 
assigned. The 8^ term is theoretical only with quark masses varying within their allowed ranges (126). The 
quark condensates in the 5^ term are obtained from PCAC, while the gluon condensate is determined by the 
fit. The main contributions to the theoretical uncertainty on a s (m 2 ) are listed in Table VII. The fact that 
the x 2 value of the axial-vector fit is better with FOPT than with CIPT should not — in our view — affect 
the conclusion in favor of the latter method, reached in Sections VII. B. 8 and VII. B. 3 on the basis of a better 
perturbative behavior, but it may require further studies. 
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TABLE VI Fit results (ALEPH Coll., 2005) for tts(mj) and the nonperturbative contributions for vector, axial- vector 
and V + A combined fits using the corresponding experimental spectral moments as input parameters. Where two 
errors are given the first is experimental and the second theoretical. The 8^ term is theoretical only with quark masses 
varying within their allowed ranges (see text). The quark condensates in the term are obtained from PCAC, while 
the gluon condensate is determined by the fit. The total nonperturbative contribution is the sum <5np = 
Full results are listed only for the CIPT perturbative prescriptions, except for a s (m 2 ) where both CIPT and FOPT 
results are given. 



Parameter 


Vector (V) 


Axial- Vector (A) 


V ± A 


a s (m 2 T ) (CIPT) 
a s K) (FOPT) 


0.355 ± 0.008 ± 0.009 
0.331 ± 0.006 ±0.012 


0.333 ± 0.009 ± 0.009 
0.327 ±0.007 ±0.012 


0.350 ± 0.005 ± 0.009 
0.331 ±0.004 ±0.012 


5 m (CIPT) 


(-3.3 ±3.0) x 10~ 4 


(-5.1 ±3.0) x 10" 4 


(-4.4 ±2.0) x 10~ 4 


(a a GG) (GeV 4 ) (CIPT) 


(0.4 ±0.3) x 10~ 2 


(-1.3 ±0.4) x 10~ 2 


(-0.5 ±0.3) x 10~ 2 


$W (CIPT) 


(4.1 ± 1.2) x 10~ 4 


(-5.7 ±0.1) x 10~ 3 


(-2.7 ±0.1) x 10~ 3 


S (6) (CIPT) 


(2.85 ±0.22) x 10~ 2 


(-3.23 ±0.26) x 10~ 2 


(-2.1 ± 2.2) x 10~ 3 


5 (8) (CIPT) 


(-9.0 ±0.5) x 10~ 3 


(8.9 ±0.6) x 10~ 3 


(-0.3 ±4.8) x 10~ 4 


Total 5 NP (CIPT) 


(1.99 ±0.27) x 10~ 2 


(-2.91 ±0.20) x 10~ 2 


(-4.8 ± 1.7) x 10~ 3 



X 2 /DF (CIPT) 0.52 4.97 3.66 



TABLE VII Sources of theoretical uncertainties and their impact on R t ,v+a and a s (m 2 ) (ALEPH Coll., 2005). 



Error source 


Value ± A th 


CIPT 


A th R T ,v+A 


FOPT 


CIPT 


A th a s (m 2 T ) 


FOPT 


Sew 


1.0198 ±0.0006 




0.002 






0.001 




\v ud \ 


0.9745 ± 0.0004 




0.002 






0.001 




K 4 


25 ± 25 


0.015 




0.009 


0.007 




0.003 


R-scale \x 


m 2 ± 2 GeV 2 


0.012 




0.035 


0.006 




0.011 




Total errors 


0.019 




0.036 


0.009 




0.012 



The final result on a s (m^.) given by (ALEPH Coll., 2005) is taken as the arithmetic average of the CIPT and 
FOPT values given in Table VI, with half of their difference added as additional theoretical error a s (m^) = 
0.340 ± 0.005 oxp ± 0.014th- The first error accounts for the experimental uncertainty and the second error 
gives the uncertainty related to the theoretical prediction of the spectral moments including the ambiguity 
between CIPT and FOPT. 

There is a remarkable agreement within statistical errors between the a s (m^) determinations using the 
vector and axial- vector data. This provides an important consistency check of the results, since the two corre- 
sponding spectral functions are experimentally independent and manifest a quite different resonant behavior. 

The advantage of separating the vector and axial-vector channels and comparing to the inclusive V ± A 
fit becomes obvious in the adjustment of the leading nonperturbative contributions of dimension D = 6 and 
D = 8, which approximately cancel in the inclusive sum. This cancellation of the nonperturbative terms 
increases the confidence in the a s (m^) determination from the inclusive V + A observables. The gluon 
condensate is determined by the first k = 1, I = 0, 1 moments, which receive lowest order contributions. It is 
observed that the values obtained in the V and A fits are not very consistent, which could indicate problems in 
the validity of the OPE approach used once the nonperturbative terms become significant. Taking the value 
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obtained in the V + A fit, where nonperturbative effects are small, and adding as systematic uncertainties 
half of the difference between the vector and axial-vector fits as well as between the CIPT and FOPT results, 
ALEPH measures the gluon condensate to be 

(a s GG) = (0.001 ± 0.012) GcV 4 . (135) 

This result does not provide evidence for a nonzero gluon condensate, but it is consistent with and has 
comparable accuracy to the independent value obtained using charmonium sum rules and e + e _ data in the 
charm region, (0.011 ± 0.009) GcV 4 in a combined determination with the c quark mass (Ioffe and Zyablyuk, 
2003). 

The D — 6, 8 nonperturbative contributions are obtained after averaging the FOPT and CIPT values: 
5 { y ] = (2.8 ±0.3) x 10- 2 , 5 { y ] = (-8.8 ±0.6) x lO" 3 , 

<^ 6) = (-3.1 ±0.3) x 10- 2 , <^ 8) = (8.7 ±0.6) x lO" 3 , (136) 
tfy+A = (-1.8 ±2.4) x 10- 3 , 5 { y ] +A = (0.5 ±5.1) x 10~ 4 . 

The remarkable feature is the approximate cancellation of these contributions in the V+A case, both for D = 6 
and D — 8. This property was predicted (Braaten et al., 1992) for D = 6 using the simplifying assumption 

of vacuum saturation for the matrix elements of four-quark operators, yielding 5^/6^ = —7/11 = —0.64, 
in fair agreement with the above results which reads — 0.90 ± 0.18. The estimate (Braaten et at, 1992) for 

Sy ^ = (2.5 ±1.3) x 10~ 2 agrees well with the experimental result. 

The total nonperturbative V + A correction, ^Np.y+A = (—4.3 ± 1.9) x 10~ 3 , is an ordcr-of-magnitude 
smaller than the corresponding values in the V and A components, (5np.v = (2.0 ± 0.3) x 10~ 2 and <5np,t1 = 
(-2.8 ±0.3) x 10~ 2 . 

2. The OPAL results 

The OPAL analysis (OPAL Coll., 1999b) makes use of the spectral functions shown in Fig. 4 and provides 
also a complete description in terms of QCD of their vector and axial-vector data. It proceeds along very 
similar lines as the ALEPH analysis, with the difference that OPAL has chosen to quote their results separately 
for three prescriptions for the perturbative expansion, respectively CIPT, FOPT and with renormalon chains, 
i.e. 



a s (m 2 ) 



0.348 ± 0.009 oxp ± 0.019th CIPT , 

0.324 ±0.006 cxp ± 0.013th FOPT, (137) 
0.306 ± 0.005 exp ± 0.011th renormalon . 



The OPAL results using CIPT and FOPT agree with those from ALEPH. The quoted theoretical uncer- 
tainties are consistently treated between the two analyses , with the important difference that the more recent 
ALEPH analysis (ALEPH Coll., 2005) uses a reduced uncertainty range for the K4 parameter (K4 = 25 ±25) 
as compared to K 4 = 50 ± 50 in the 1998 analysis (ALEPH Coll., 1998c) and K 4 = 25 ± 50 in OPAL. The 
justification of the smaller K4 uncertainty was provided in Section VII.B.l. As for the value obtained with 
the renormalon chains, we refer to our discussion in Sections VII. B. 7 and VII. B. 8 where we concluded that 
however interesting this approach is, its performance for a reliable and precise determination of a s (m 2 ) is not 
adequate. OPAL also derives values for the gluon condensate: averaging the results of the CIPT and FOPT 
analyses, one obtains {a s GG) = (0.005 ± 0.017) GeV 4 , in agreement with ALEPH. 

3. Running of a s (s) below 

Using the spectral functions, one can simulate the physics of a hypothetical r lepton with a mass t/so 
smaller than m T through Eq. (82). Assuming quark-hadron duality, the evolution of R T (so) provides a direct 
test of the running of a s (sn), governed by the RGE /^-function. On the other hand, it is also a test of the 
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FIG. 22 Left: The ratio R Tt v+A versus the square "r mass" so. The curves are plotted as error bands to emphasize 
their strong point-to-point correlations in so. Also shown is the theoretical prediction using CIPT and the results for 
R Tt v+A and the nonperturbative terms from Table VI. Right: The running of a s (so) obtained from the fit of the 
theoretical prediction to R T: v+a{so) using CIPT. The shaded band shows the data including only experimental errors. 
The curve gives the four-loop RGE evolution for three flavors. 



stability of the OPE approach in r decays. The studies performed in this section employ only CIPT. Results 
obtained with FOPT are similar 26 and differ mostly in the central a s (m 2 ) value. 

The functional dependence of R T y + A(so) is plotted in the left hand plot of Fig. 22 together with the 
theoretical prediction using the results of Table VI. The spreads due to uncertainties are shown as error 
bands. The correlations between two adjacent points in sg are large as the only new information is provided 
by the small mass difference between the two points and the slightly modified weight functions. Moreover 
they are reinforced by the original experimental and theoretical correlations. Below 1 GeV 2 the error of the 
theoretical prediction of R t ,v+a{so) starts to blow up due to the increasing uncertainty from the unknown 
Ki perturbative term; errors of the nonperturbative contributions are not contained in the theoretical error 
band. Figure 22 (right) shows the plot corresponding to Fig. 22 (left), translated into the running of a s (sn), 
i.e., the experimental value for a s (sn) has been individually determined at every sq from the comparison of 
data and theory. Also plotted is the four-loop RGE evolution using three quark flavors. 

It is remarkable that the theoretical prediction using the parameters determined at the r mass and 
Rt,v+a(so) extracted from the measured V + A spectral function agree down to s ~ 0.8 GeV 2 . The agree- 
ment is good to about 2% at lGeV 2 . This result, even more directly illustrated by the right hand plot of 
Fig. 22, demonstrates the validity of the perturbative approach down to masses around 1 GeV, well below the 
r mass scale. The agreement with the expected scale evolution between 1 and 1.8 GeV is an interesting result, 
considering the relatively low mass range, where a s is seen to decrease by a factor of 1.6 and reaches rather 
large values ~ 0.55 at the lowest masses. This behavior provides confidence that the a s (m 2 ) measurement is 
on solid phenomenological ground. 



4. Final assessment on the a s (m 2 ) determination 

Although the recent ALEPH evaluation of a s (w 2 ) represents the state-of-the art, several remarks can be 
made: 

• The analysis is based on the ALEPH spectral functions and branching fractions, ensuring a good con- 
sistency between all the observables, but not exploiting the full experimental information currently 



As already indicated by the smaller \ 2 value of the spectral moments fit (see Table VI), the agreement between the FOPT 
evolution computation and the data is better than it is when using CIPT. 
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available from other experiments. Since the result on a s (m 2 ) is limited by theoretical uncertainties, 
one should expect only a small improvement of the final error in this way, however it can influence the 
central value. 

• One example for this is the evaluation of the strange component. Some discrepancy is observed between 
the ALEPH measurement of the {K-k-k^v mode and the CLEO and OPAL results, as discussed in 
Section III.C. Although this could still be the result of a statistical fluctuation, their average provides a 
significant shift in the central value compared to using the ALEPH number alone. Another improvement 
is the substitution of the measured branching fraction for the K~v mode by the more precise value 
predicted from r-fi universality (see Eq. (16)). Both operations have the effect to increase the strange 
R T .s ratio from O.f 603 ± 0.0064, as obtained by ALEPH, to 0.1686 ± 0.0047 for the world average. 

• One can likewise substitute the world average value for the universality-consolidated value of the elec- 
tronic branching fraction given in Eq. (13), £>" m = (17.818 ± 0.032)%, to the corresponding ALEPH 
result, B e = (17.810 ± 0.039)%, with little effect on the central value, but some improvement in the 
precision. 

• Most of the theoretical uncertainty originates from the limited knowledge of the perturbative expansion, 
only predicted to third order. This problem was discussed in detail in Section VII. B.l and, among the 
standard FOPT and CIPT approaches, strong arguments have been presented in favor of the latter 
one. We therefore take the result from the CIPT expansion, not introducing any additional uncertainty 
spanning the difference between FOPT and CIPT results. The dominant theoretical errors are from the 
uncertainty in K\ and from the renormalization scale dependence, both covering the effect of truncating 
the series after the estimated fourth order. 

From this analysis, one finds the new value for the nonstrange ratio, 

Rt,V+A = Rt — Rt,S 

= (3.640 ±0.010) - (0.1686 ±0.0047) 

= 3.471 ±0.011 , (138) 

to be compared to the ALEPH value of 3.482±0.014. The result (138) translates into an updated determination 
of a s {m 2 ) from the inclusive V + A component using the CIPT approach 

a s (m 2 ) = 0.345 ± 0.004 cxp ± 0.009 th , (139) 

with improved experimental and theoretical precision over the ALEPH result. 

5. Evolution to M| 

It is customary to compare a s values, obtained at different renormalization scales, at the scale of the Z- 
boson mass. To do this, one evolves the a s result using the RGE (75). A difficulty arises from the quark 
thresholds: in the MS scheme, or any of its modifications such as MS, the beta function governing the 
running of the strong coupling constant is independent of quark masses. This simplifies the calculation of 
QCD corrections beyond the one-loop level. On the other hand, decoupling of heavy quarks (Appclquist and 
Carazzone, 1975) is not manifest in each order of perturbation theory (Bcrnreuther, 1983). A heavy quark 
decoupling is realized in momentum-space dependent renormalization schemes (MO). The RGE in an MO 
scheme involves the scaling- function /3mo = /?(o4, a') that depends on the coupling a' s (p?), the quark 
masses m-(/u 2 ) and a gauge parameter a'(/x 2 ), where the primes stand for the scheme dependence of the 
parameters. Quark-loop calculations in this scheme appear rather complicated, however the mass dependence 
of the /3-function provides a suppression of heavy-quark effects at scales much smaller than the masses of these 
quarks, i.e., it decouples heavy quarks from the light-particle Green's function to each order in perturbation 
theory. 

To obtain decoupling in MS schemes, one builds in the decoupling region, [i <C mj n/ ' (p 2 ), where m^ n f\^ 2 ) 
is the running MS mass of the heavy quark with flavor /, an effective field theory that behaves as if only 
light quarks up to flavor rif — 1 were present. Matching conditions connect the parameters of the low energy 
effective Lagrangian to the full theory. The coupling constant of the effective theory can then be developed 
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in a power series of the coupling constant of the full theory with coefficients that depend oni = ln(/i 2 /m^). 
Doing so one obtains for the matching of a s and the light quark masses between the "light" flavor n/ — 1 and 
the heavy-quark flavor rif (Bernreuther and Wetzel, f982; Chetyrkin et al, f997b, f998a; van Ritbergen et 
al, 1997a; Rodrigo et al, f998; Wetzel, 1982) 



a< n 'V) - ai" /_ V) 
m< n/) (M 2 ) - m< n, - 1) (p 2 ) 



C k (x) (a^V))' 



OO 

E 

k=i 

OO ^ 

^H k {x) (a< n '"V)) 
fe=i 



(140) 



There is no straightforward choice of the matching scale /x. However, to acquire good convergence of the 
perturbative expansion, one should satisfy fi/m q ~ 0(1). Since vacuum polarization of quark pairs modifies 
the coupling constant, we believe that the scale fi = 2m q is quite meaningful. The effect from this scale 
ambiguity is part of the systematic error assigned to the RGE evolution (see below). 

The functions Ck{x) and Hk(x) are derived by inserting the RGE-based Taylor expansions (Eq. (81) and 



equivalent for the quark masses) into Eq. (140), and solving for each order in a 
coupled differential equations (Rodrigo et al., 1998) 
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where the integration coefficients c^o, c?i,o arc computed in the MS scheme at the scale of the quark masses 
(i.e., m q — m q (m q )) (Bernreuther and Wetzel, 1982; Chetyrkin et al, 1997b; Rodrigo et al., 1998; Wetzel, 
1982) 
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The evolution of the a s (m^) measurement from the inclusive V + A observables given in Eq. (139), based 
on Runge-Kutta integration of the RGE (75) to N 3 LO (see Section VILA), and three-loop quark-flavor 
matching (140), gives 



a s (M|) = 0.1215 ±0.0004 cxp ± 0.0010th ± 0.0005 O voi , 
= 0.1215 ±0.0012 . 



(146) 



The first two errors originate from the a s (m 2 ) determination given in Eq. (139), and the last error stands for 
ambiguities in the evolution due to uncertainties in the matching scales of the quark thresholds (Rodrigo et 
al, 1998). This evolution error receives contributions from the uncertainties in the c-quark mass (0.00020, m c 
varied by ±0.1 GeV) and the 6-quark mass (0.00005, nib varied by ±0.1 GeV), the matching scale (0.00023, 
H varied between 0.7 m q and 3.0m g ), the three-loop truncation in the matching expansion (0.00026) and the 
four-loop truncation in the RGE equation (0.00031), where we used for the last two errors the size of the 
highest known perturbative term as systematic uncertainty. These errors have been added in quadrature. 
The result (146) is a determination of the strong coupling at the Z mass scale with a precision of 1%. 

The evolution path of a s (m^.) is shown in the upper plot of Fig. 23. The two discontinuities are due to 
the quark-flavor matching at [i = 2m q . One could prefer to avoid the discontinuities by choosing fi = m q so 
that the logarithms in Eq. (140) vanish, and the matching becomes (almost) smooth. However, in this case, 
one must first evolve from m T down to m c to match the c-quark flavor, before evolving to Tn^. The effect on 
<x s {M\ ) from this ambiguity is within the assigned systematic uncertainty for the evolution. 
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FIG. 23 Top: The evolution of a s {m 2 T ) (139) to higher scales fi using the four-loop RGE and the 3Toop matching 
conditions applied at the heavy quark-pair thresholds (hence the discontinuities at 2m c and 2m,b). The evolution is 
compared with other independent measurements (see text) covering scales varying over more than two orders magni- 
tude. The experimental values are briefly discussed in the text; they are mostly taken from the compilation (Bethke, 
2004). Bottom : The corresponding extrapolated a s values at Mz- The shaded band displays the r decay result within 
errors. 



6. Comparison with other determinations of as(M§ ) 

The evolution of a s (m^) in Fig. 23 is compared with other independent measurements compiled in (Bethke, 
2004). Listed below are some outstanding measurements made over a large energy range in a variety of 
processes with different experimental and theoretical precisions: 

DIS (Bj-SR): a s (1.58GeV) = 0-375±g;g|? (Ellis and Karliner, 1995) was obtained using the polarized struc- 
ture function data for protons gf(x) and neutrons g"(x) from the CERN SMC (Adams et al., 1994) 
and SLAC E142 (E143 Coll., 1995) experiments, based on the Bjorken sum rule (Bjorken, 1970) at 
N 2 LO (Larin et al., 1991; Larin and Vermaseren, 1991) 
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Here gv and gA are the vector and axial- vector coupling constants of the neutron decay, \gA/gv\ = 
— 1.2573 ± 0.0028 (determined assuming CVC), and the last two coefficients are calculated for three 
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active flavors. No explicit corrections for nonperturbative higher twist effects 27 were applied to derive 
this result. However an estimate of the size of the perturbative 0{oq ) term was taken into account. 



DIS (GLS-SR): a s (1.73 GeV) = 0.280 ± 0.061 cxp _o!o3olth is an u P date 01 thc analysis using structure 
function F$ data in the Q range from 1 to 15.5 GeV 2 performed by the CCFR collaboration (CCFR 
Coll., 1998). It uses the Gross-Llewellyn-Smith (GLS) sum rule (Gross and Llewellyn-Smith, 1969) at 
N 2 LO (Chyla and Kataev, 1992; Larin et a/., 1991; Larin and Vermaseren, 1991) 
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(148) 



given for three active flavors. The systematic uncertainty is dominated by the extrapolation of the GLS 
integral to the region x < 0.01 where no measurements exist, and to x > 0.5, which is substituted by 
the structure function F 2 from SLAC data (Withlow et at, 1992). The theoretical error is dominated 
by uncertainties in the higher twist corrections. 

Y sum rules: a s (mb)\m b =4.75GeV = 0.217 ± 0.021 was found together with the bottom quark pole mass 
rrib from sum rules for the T system in N 2 LO, resumming all O (a|, ct s v, v 2 ) terms with v being the 
velocity of the heavy quark (Penin and Pivovarov, 1998). The precision of this result on a s (nib) was not 
confirmed by a similar analysis (Hoang, 1999), where the theoretical parameter space has been scanned, 
which leads to a more conservative estimate of the theory uncertainty. 



e+e (<T had ): 



e+e 



a s (10.52GeV) = 0.20±0.06 was determined in N 2 LO by CLEO (CLEO Coll., 1998) using the 
hadronic cross section according to 
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where Qi are the electric charges of quark flavors i that are produced at center-of-mass energy just below 
the T(45) resonance in the e + e _ continuum. Another measurement a s (42.4GeV) = 0.175 ± 0.028 was 
obtained with the use of data from PETRA and TRISTAN in the center-of-mass energy range from 20 
to 65 GeV (Bethke, 2000; Haidt, 1995). 

" (jet & event-shape): Using data from PETRA, TRISTAN, SLC and LEP1.2 over a large energy 
range, several determinations of a s were obtained based on thc jet rate and event-shape variables 
that are sensitive to gluon radiation governed by the strong coupling constant. The observables arc 
calculated in resummed NLO (i.e., NLO matched with next-to-leading order logarithms). In all these 
determinations, the theoretical uncertainty from the renormalization scale uncertainty is the dominant 
error source. At low PETRA energy the hadronization uncertainty is also important. See references 
in (Bethke, 2000). 



pp — bbX: a s (20GeV) - 0.145 tHloLxp 



1-0.013 1 
-0.016 Ith 



was obtained by UA1 (UA1 Coll., 1996) from a measure- 
ment of the cross section of the process pp — ► bbX for which NLO QCD predictions exist. The theoretical 
error includes uncertainties due to different sets of structure functions, renormalization/factorization 
scale uncertainties and the 6-quark mass. 

pp,pp jX: a s (24.3GcV) = 0.135 ± 0.006 oxp taoosL was determined by UA6 (UA6 Coll., 1999) in NLO 
from a measurement of the cross section difference cr(pp — ► jX) — u{pp — » jX) such that the poorly 
known contributions of the sea quarks and gluon distributions in the proton cancel. The theoretical error 
includes uncertainties from the scale choice and from the variation of the parton distribution functions. 

e+e~ (Z width): a s (M z ) Zwidth = 0.1186 ± 0.0027 was determined from the hadronic width at the Z° 
resonance by a global fit to all electroweak data in N 2 LO (LEP EWWG, 2004). 



Interactions involving more than one parton, like secondary interactions of quarks with the target remnant, give rise to 
additional contributions to the structure functions. These "higher twist" terms are classified according to their Q 2 dependence 
1/Q 2n , where n = is leading twist (twist two), n = 1 is twist four, etc. 
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A comparison is best achieved by extrapolating these measurements to Mz as shown in the lower plot in 
Fig. 23. The most precise result at the Mz scale stems from r decays (146), which is a 1.1% determination 
limited in accuracy by theoretical uncertainties in the perturbative expansion. The significant improvement 
in precision compared to the previous ALEPH result (ALEPH Coll., 1998c), 0.1202 ± 0.0008 exp ± 0.0024 th ± 
O.OOlOovoi, is due to the higher statistics and the more detailed experimental analysis, but mostly because of 
the smaller theory uncertainty assigned to the perturbative expansion and our conclusion to lift the ambiguity 
between CIPT and FOPT. 

Another precise determination (DIS {e//j,;F 2 ) [1.9 - 15.2] GeV), a s (M§) = 0.1166 ± 0.0009 sta t ± 
0.0020 sys t (Bethkc, 2003; Santiago and Yndurain, 2001), obtained using the structure function F 2 from deep 
inelastic electron (e) and muon (/x) scattering (DIS) data in the Q 2 range between 3.5 GeV 2 and 230 GeV 2 is 
also shown. The analyses use N 2 LO calculations wherever available. The systematic error is mostly theoretical 
including higher twist effects and an estimate of the N 3 LO corrections, and is doubled to account for missing 
contributions associated with the renormalization scale and scheme uncertainties. 

In addition to the measurements shown in Fig. 23, there is now a rather precise determination from lattice 
QCD simulations (Davis et al, 2004): a s (M|) = 0.121 ± 0.003. This is the first lattice determination with 
realistic quark vacuum polarization and an 0(a 2 ) improved staggered-quark discretization of the light-quark 
action. 



7. A measure of asymptotic freedom between m 2 T and M| 

The r-decay and Z-width determinations have comparable accuracies, which are however very different in 
nature. The r value is dominated by theoretical uncertainties, whereas the determination at the Z resonance, 
benefiting from the much larger energy scale and the correspondingly small uncertainties from the truncated 
perturbative expansion, is limited by the experimental precision on the electroweak observables, essentially 
the ratio of leptonic to hadronic peak cross sections. The consistency between the two results provides the 
most powerful present test of the evolution of the strong interaction coupling, over a range of s spanning more 
than three orders of magnitude, as it is predicted by the nonabelian nature of the QCD gauge theory. The 
difference between the extrapolated r-decay value and the measurement at the Z is: 

a T s( M z) ~ ot z s {M 2 z ) = 0.0029 ± 0.0010 r ± 0.0027 z , (150) 

which agrees with zero with a relative precision of 2.4%. 

In fact, the comparison of these two values is valuable since they are among the most precise single mea- 
surements and they are widely spaced in energy scale. Thus it allows one to perform an accurate test of 
asymptotic freedom. Let us consider the following evolution estimator (DELPHI Coll., 2004) for the inverse 
of a s (s), 

r( Sl , S2) = 2.^M^M, (151) 
msi — ms2 

which reduces to the logarithmic derivative of 1 (s) when si — > ,s 2 , 

da' 1 2ir/3(s) 



2/3 , Pi a s 



(1 + ^— + ■■■), (152) 

with the notations of Eq. (75). At first order, the logarithmic derivative is driven by (3 . 
The r and Z experimental determinations of a s (s) yield the value 

r oxp (m 2 , M|) = 1.405 ± 0.053 , (153) 

which agrees with the QCD prediction using the RGE to N 3 LO, and three-loop quark-flavor matching (140), 
as discussed in Section VII. C. 5, 

r QC D(TO 2 ,M|) = 1.353 ±0.006 . (154) 

This is to our knowledge the most precise experimental test of the asymptotic freedom property of QCD at 
present. It can be compared to an independent determination (DELPHI Coll., 2004), using an analysis of 
event shape observables at LEP between the Z energy and 207 GeV, r(M|, (207 GeV) 2 ) = 1.11 ± 0.21, for a 
QCD expectation of 1.27. 
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VIII. SPECTRAL FUNCTIONS OF TAU FINAL STATES WITH STRANGENESS 

Complete experimental studies of the Cabibbo-suppressed decays of the r became available at LEP ( ALEPH 
Coll., 1999b), allowing to initiate systematic QCD studies (Chen, 1998; Chen et al., 1999; Davier, 1997; Davier 
et al., 2001; Hocker, 1997) of the effect induced by the strange quark mass in the r decay width (Braaten 
et al, 1992; Chctyrkin and Kwiatkowski, 1993; Chetyrkin and Kiihn, 1997; Chetyrkin et al, 1997a, 1998b; 
Kambor and Maltman, 2000, 2001; Korncr et al, 2001b; Maltman, 1998a; Pich and Prades, 1998, 1999, 2000; 
Prades and Pich, 1999; Prades, 1999). Indeed, from the separate measurement of the strangeness S = and 
S = — 1 r decay widths, it is possible to pin down the SU(3) breaking effects and to obtain information on the 
strange quark mass. Recent experimental information (CLEO Coll., 1999b; OPAL Coll., 2004) can now be 
included in addition to the ALEPH data. In this section the current situation on the physics of r decays with 
net strangeness is reviewed, with emphasis on the determination of the strange quark mass and the CKM 
element \V US \. We anticipate here that while the determination of \V US \ through QCD sum rules is found to 
be robust, a reliable and competitive determination of the strange quark mass m s is difficult because of poor 
convergence of the mass-dependent perturbative series, and lack of sufficient inclusiveness of the higher-order 
spectral moments that are most sensitive to m s . 



A. Strange spectral functions 

The experimental determination of nonstrange and strange (subscript S) spectral functions has been pre- 
sented and discussed in Section IV. As shown in Section VII.C, one uses the spectral functions, or equivalcntly 
the fully corrected mass-squared distributions, to construct the spectral moments 

^W^-^te)^* (155) 



where i?°° (s) = R T ,(S)- 

To reduce the theoretical uncertainties and emphasize the SU(3) breaking component originating from 
the larger s quark mass, it is convenient to consider the difference between nonstrange and strange spectral 
moments, properly normalized with their respective CKM matrix elements (ALEPH Coll., 1999b) 



6Rk/ " w^ Rk ^ s =° ~ r^F^- 1 ' (156) 



for which the massless perturbative contribution vanishes. The contributions from the various decay modes 
to 5R®° are visualized in Fig. 24. 

B. Theoretical analysis of R Ti s 

The inclusive r decay ratio into strange hadronic final states 

r(r~ — > hadronSo_ , v T ) 

Rr,s = 1 r( _ __ S ~' 1 ' , 157 

T(r — > e u e u T ) 

can be used through the framework of the OPE to determine m s (s) at the scale s = m\ (Chen, 1998; Chen 
et al., 1999; Davier, 1997; Hocker, 1997). However, the perturbative expansion computed for the mass term 
in (Chetyrkin and Kwiatkowski, 1993) was shown to be incorrect (Maltman, 1998a). After correction the 
series shows a problematic convergence behavior (Chetyrkin et al., 1998b; Maltman, 1998a; Pich and Prades, 
1998; Prades and Pich, 1999; Prades, 1999), which is discussed in detail below. 

Following (Braaten et al, 1992) we recall the theoretical prediction for the inclusive vector plus axial- vector 
hadronic decay rate, given by 



R T {mi) = 12ttS ew I % ( 1 - -Ay 
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(158) 
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FIG. 24 Integrand of Eq. (155) with (k = 0,1 = 0) for the difference (156) of the Cabibbo-corrected nonstrange 
and strange invariant mass spectra from (ALEPH Coll., 1999b). The contribution from massless perturbative QCD 
(pQCD) vanishes. To guide the eye, the solid line interpolates between the bins of constant 0.1 GeV 2 width. 



with the two-point correlation functions U.^ — \V u d\ 2 ^-udv+A + l^"s| 2 n^Py +/1 f° r the hadronic final state 
of spin J. The choice of the particular spin combination for the correlators taken in Eq. (158) is justified 
in (Braaten et al, 1992). Using the OPE (Shifman et al, 1979), Eq. (158) for the strange part can be 
decomposed as 

Rr.s = 3\v us \ 2 s EW ( i + <y<°> + s' EW + 4 2 < m « } + ]T si s D) 1 • ( 159 ) 

£>=4,6,... 

The 5^ term is the perturbative part of mass dimension D = 0, known to third order in the expansion 
with a s — a s (m 2 .)/7r, and discussed in detail in Section VII. B.l. The residual non-logarithmic electroweak 
correction S' EW ~ 0.0010 (Braaten and Li, 1990) is neglected in the following. Throughout this section, all 
QCD observables are expressed in the MS renormalization scheme. 

From the result obtained in Section VII. C. 4 for the world average branching fraction for inclusive r 
decays to strange hadronic final states (improved using lepton universality for the Kv mode), one ob- 
tains R T> s = 0.1686 ± 0.0047, using the value for B^ nl given in Section III.B. This result can be com- 
pared with the QCD prediction for a massless strange quark neglecting the nonperturbative contributions, 
R^ s = 3\V us \ 2 Sew (1 + <5 (0) ) = 0.1847 ± 0.0029, where the quoted error reflects the uncertainties on mostly 
\V US \ and to a lesser extent a s (m 2 .). The two values differ by 2.9ct, in the direction predicted for a non-zero 
m s value. 

The 4 2 ' m,) term in Eq. (159) is the mass contribution of dimension D — 2, which is in practice only relevant 
for the strange quark mass. Fixed-order perturbation theory (FOPT) gives (Chetyrkin and Kwiatkowski, 1993; 
Maltman, 1998a) 



(2,m s ) = 8 w s (m r ) 

S ™ 2 



1 + y a s + 46.00 a 2 + ( 283.6 + | 4 1+ ° } ) a l + 



(160) 



The third order coefficient x^ +0) occurring in the expansion of the massive J = 1 + correlator has been 
recently calculated (Baikov et al, 2004), while the J — correlator was already known up to third order. 
Previous analyses used an estimate for Sg 1+0 \ however it turns out to be of minor numerical importance for 
the result on m s . The true value (Baikov et al, 2004) is 4 1+0) = 202.3, while the naive estimate, assuming 
a geometric growth of the perturbative coefficients, gives a;^ ^ — x^ +0 \x^ + ^ / x^ +a ^ ) ~ 165, which was 
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taken with a 200% uncertainty. Setting a s {rn 2 r ) = 0.334, one finds that Eq. (160) docs not exhibit good 
convergence. 

A contour-improved analysis (CIPT) for the dimension D = 2 contribution has been used in (Chetyrkin 
et al, 1998b; Pich and Prades, 1998; Prades and Pich, 1999; Prades, 1999). It consists of a direct numerical 
evaluation of the contour integral derived from Eq. (158), using the solution of the renormalization group 
equation (RGE) to four loops as input to the running a s (s) and m s (s). This provides a resummation of 
all known higher order logarithmic integrals and has been found to improve the convergence of the masslcss 
pcrturbative series. This is also the case for Eq. (160), however not solving the convergence problem. Inde- 
pendently of whether FOPT or CIPT is used, the origin of the convergence problem is found in the J = 
component as defined in Eq. (158). 

The phenomenological analysis is performed using the OPE framework including the light quark masses 
and the D — 4 quark-mass corrections. Higher dimension terms are assumed to be numerically insignificant 
compared to the present experimental uncertainty so that for all practical purposes the strange quark mass 
can be obtained from the relation 



TO' 



;K) * 



kt 



5R k T 
24Sew 



2tt 



2 (* Q4K)) 

mi 



Qki(a s ) 



(161) 



where Ake{a s ) and Qu{a s ) are the pQCD series, defined in (Pich and Prades, 1999, 2000), associated with 
the D = 2 and D = 4 contributions to 5R* e , and 



(SOi{m 2 r )) = (0\m s ss -m d dd\0}(m 2 T ) 
~ -(1.5 ±0.4) x 10~ 3 GeV 4 



(162) 



As discussed above, the QCD scries for = A^ +0 ^ + A)^ are problematic, since they exhibit bad 

convergence originating from the longitudinal component (A^) (Maltman, 1998a; Pich and Prades, 1998, 
1999, 2000; Prades and Pich, 1999; Prades, 1999). The value of the strong coupling constant is taken to 
be a s (m^) = 0.334 ± 0.022, from the analyses of the nonstrange V + A moments R ke . Correlations between 
a s (m^) and the moments 5Rn are negligible since the uncertainty on the former value is dominated by theory 
and on the latter one by the measurement of the strange component. Using these inputs, the pQCD series 
Ake can be displayed up to third order using CIPT, separating the (1 + 0) and (0) series (a s = 0.106): 
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The exhibited behavior is that of asymptotic series close to their point of minimum sensitivity and a prescrip- 
tion is needed to evaluate the expansions and to make a reasonable estimate of their uncertainties. 



C. Several approaches 

The bad convergence properties of the perturbative QCD expansion can be alleviated in several ways. 
For instance, one could try to work with the full series up to the highest known order (third), observing 
that the terms in the overall expansion show a reasonable rate of decrease, at least in the CIPT approach. 
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Fourth order coefficients can be estimated, using for instance the principle of minimal sensitivity (Stevenson, 
1981) or the simplistic geometric growth rule. The latter method was used in early analyses (Chen, 1998; 
Chen et al, 1999; Davier, 1997; Hocker, 1997). 

Examination of examples of asymptotic series suggests that a reasonable procedure is to truncate the 
expansion where the terms reach their minimum value. The precise prescription — cutting at the minimum, or 
one order before or after, including the full last term or only a fraction of it — is arbitrary and this ambiguity 
must be reflected by a specific uncertainty attached to the procedure. In (Chen et al., 2001; Davier et al, 2001) 
the adopted rule was to keep all terms up to (and including) the minimal one and to assign as a systematic 
uncertainty the value of the last term retained. It follows that the A^o series are summed up to third order 
for k = 0, 1, second order for k = 2 and first order for k = 3, 4. 

Since the culprit is the longitudinal expansion, it could be wise to remove the longitudinal piece beforehand 
from the inclusive data and to apply the QCD analysis only to the well-behaved (1 + 0) part 28 . This method 
was already used in (ALEPH Coll., 1999b), where it was noted that the dominant contribution to the strange 
J = component is given by the kaon pole. Other longitudinal contributions can be identified and estimated, 
giving 

R f.s( K ) = -0.00615 ± 0.00026 , (164) 
R^Js (other 0",0 + ) = -0.0015 ± 0.0015 , (165) 

leading to the result 

6R°°{ J = 0) = 0.147 ±0.031 , (166) 

where the J = contributions in the nonstrange part are assumed to be saturated by the pion pole and have 
been computed to be —0.0078. Similarly, one estimates the J = contributions for all (k,£) moments and 
subtracts them from the measured moments. 

The more recent analyses rely on phenomcnological parameterizations of the additional scalar and pseu- 
doscalar contributions (Jamin et al, 2002; Maltman and Kambor, 2002). The O^ 1 spectral functions con- 
structed in this way rely on not so well-established resonances with poorly known couplings, however they 
probably provide a realistic level for these contributions. Even with their intrinsically large uncertainties they 
are an order of magnitude more precise than the experimental estimate from ALEPH (Eq.(165)). In the 
analysis of (Gamiz et al, 2003) the corresponding contribution to R TiS yields 

5R®°(J = 0, ph) = 0.155 ± 0.005 , (167) 

where "ph" stands for phenomenological parameterizations. The quoted uncertainty in this method is clearly 
much smaller than in Eq. (166). It reflects the fact that pieces of the spin-0 spectral functions can be 
reconstructed using known resonances. 

1. Verifying the stability of the results 

As shown in the analysis of the nonstrange r decays in (ALEPH Coll., 1998c), one can test the validity 
of the QCD analysis by simulating the physics of a hypothetical r lepton of lower mass, ^/sq < m T . This 
is obtained by replacing m 2 by sq everywhere in Eqs. (155, 158), and correcting the latter equation for the 
modified kinematic factor. Under the assumption of quark-hadron duality, the evaluation of the observables as 
function of So represents a test of the OPE approach, since the energy dependence of the theoretical predictions 
is determined once the parameters of the theory are fixed. The results of this exercise are given in Fig. 25, 
showing the variation with so of the first four SR^. moments, and the value for m s (so) derived from each 
moment. The bands indicate the experimental and theoretical uncertainties. The agreement between data and 
theory, perfect at so = by construction, remains acceptable for lower so values, down to 1.6 (2.4) GeV 2 
for k = (4). For the first three moments used in the final determination, the running observed in data 
follows the RGE evolution down to about 2 GeV 2 . The validity of using higher moments thus appears more 
questionable. 



A drawback of this procedure is that it reduces the inclusiveness of the data and hence makes the phenomcnological analysis 
more vulnerable to quark-hadron duality violations. 
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FIG. 25 The observables 5R^°(so) as a function of the "r- mass" -squared so confronted with the QCD predictions (left 
plots), and the derived m s (so) compared with the QCD RGE running (right plots). The bands correspond to the 
experimental and theoretical uncertainties (left plots), and experimental errors only (right plots). By construction, 
data and theory agree at so = m^. The figure is taken from (Chen et ai, 2001). 



2. Results and discussion 

Several analyses of the ALEPH strange spectral function have been performed in order to extract m s (m^). 
More recent analyses take advantage of new experimental input on the strange spectral function from (OPAL 
Coll., 2004) and improvements in branching fractions for some strange channels (see Table I). 

• In (ALEPH Coll., 1999b), a rather conservative road was followed, using only the 1+0 part and a prudent 
estimate of the removed longitudinal part directly from data. A price was paid in the experimental 
sensitivity, leading to large errors. Also, the experimental moments were fit, not only to the strange 
quark mass, but also to the values of the nonperturbative operators (SOq) and (SOs)- The value found 
by ALEPH is m s (m 2 T ) = (176 t%\ exp t 2 2 t\ th ± 14 mcth ) MeV. 

• Using the truncation method with the k = 0,1,2 moments (Chen et ai, 2001), the value m s (m 2 ) = 
(120 ± ll exp ± 8| Vus | ± 19 th ) MeV is found. 

• The analysis of (Korner et ai, 2001b) uses the ALEPH moment and obtains m s (m^) = (130 ± 
27 cxp ± 9th) MeV. It advocates contour-improved resummation and employs an effective charge as well 
as effective masses absorbing the higher perturbative terms. The small quoted theoretical uncertainty is 
not exhaustively motivated in their paper, and apparently no uncertainty is included from \V US \, which 
should be of the order of ±13 MeV. 

• Another analysis (Kambor and Maltman, 2000, 2001) using the ALEPH data makes use of weight 
functions multiplying the correlators in Eq. (158). These weights are designed to improve the convergence 
of the perturbative series, while suppressing the less accurate high-mass part of the strange spectral 
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function. This latter feature is similar to using higher moments in k and it is not entirely clear to what 
extent this procedure provides stable results and hence how reliable the answer is. Nevertheless that 
result is in good agreement with other results, with a smaller theoretical uncertainty, but apparently 
not including the effect from the renormalization scale ambiguity. 

• In (Gamiz et at, 2003), the longitudinal part is subtracted using a phenomenological prescription. The 
result to s (to 2 ) = (107 ± 18)MeV is obtained by taking the weighted average of the m s values from 
k = 0, . . . , 4 and I = moments. This procedure is statistically hardly justified as the experimental 
moments are strongly correlated. The result is more precise than other determinations, however a large 
systematic trend is observed for different k values: the m s central values decrease monotonically from 
151 for k = to 91 MeV for k = 4. No additional systematic uncertainty is included to account for 
this effect, which is larger than the final quoted error, dominated by the experimental error and the 
uncertainty on \V US \. The updated value using the OPAL strange spectral function with world-average 
branching ratios yields (Gamiz et al., 2005) a value m s (m^) = (84 ± 23) MeV. 

• Another analysis (Baikov et al, 2004) uses the full perturbative series for the 1 + correlator up to 
third order and an estimate of the fourth order contributions, while the longitudinal part is subtracted 
using the same phenomenological model as in (Gamiz et al., 2003). The result, based on the k = 3 
I = moment, is m s (m 2 ) = (100 h t{g| r t ) MeV. Here the first (quite small) error shall cover the 
uncertainties on the perturbative treatment, while the remaining error is dominated by the experimental 
and \V US \ uncertainties. The same systematic trend is observed with m s {m%.) values ranging from 
126 MeV for k = 2, to 82 MeV for k = 4. 

The determinations of m s from t decays compare reasonably well to the results of analyses of the divergence 
of the vector and axial- vector current two-point function correlators (Becchi et al, 1981; Chetyrkin et al., 1995, 
1997d; Colangelo et al., 1997; Dominguez and de Rafael, 1987; Domingucz et al., 1991; Gasser and Leutwyler, 
1982; Jamin and Miinz, 1995; Jamin, 1998; Katacv et al., 1983; Narison, 1989). In these approaches the 
phenomenological information on the associated scalar and pseudoscalar spectral functions is reconstructed 
from phase-shift resonance analyses, which are yet incomplete over the considered mass range and need to be 
supplemented by other assumed ingredients, in particular the description of the continuum. 

Another approach (Narison, 1995) considers the difference between isovector and hypercharge vector current 
correlators as related to the 1=1 and 1 = spectral functions accessible in e + e~ annihilation into hadrons 
at low energy. The author obtains m s (l GeV 2 ) = (198 ± 29) MeV. This approach was criticized in (Maltman, 
1998b), pointing out the possibility of large isospin breaking leading to a significant deviation for the extracted 
m s value. Subsequently in (Narison, 1999), new SU(3)-breaking sum rules much less affected by SU(2) breaking 
were studied, with the result m s (l GeV 2 ) = (178 ± 33) MeV. 

Finally, full lattice QCD calculations of m s including dynamical quark actions have become available re- 
cently. A pioneering analysis using the so-called staggered fermion approach (Aubin et al., 2004) finds in the 
MS scheme m s (4 GeV 2 ) = (76 ±8) MeV. The quoted error is dominated by the uncertainty in the perturbative 
theory, needed to obtain the MS mass from the pole mass, itself calculated from the lattice bare mass. The 
included systematic uncertainty from the simulation on the lattice is quoted to be 3 MeV, while the statistical 
error is negligible. Other computations using Nf = 2 Wilson fcrmions find the significantly larger values of 
m s (4GcV 2 ) = (119±10)MeV (Gockeler et al, 2004) and m s (4GeV 2 ) = (116±6)MeV (Becirevic et al., 2004). 

Available results on m s are displayed in Fig. 26. They are scaled for convenience at a mass of 2 GeV, using 
the four-loop RGE 7- function (Larin et al, 1997b). One observes a scatter among the values that exceeds the 
quoted uncertainties in some cases. 

D. Sensitivity to \V US \ and concluding remarks on the m s determinations 

Before concluding the discussion of strange r decays, it is useful to combine the available experimental 
information. Taking detector-corrected experimental mass spectra from (ALEPH Coll., 1999b), the world 
average branching fractions for the strange final states (Table I) and the value (16) for the t~ — ► K~v T 
branching ratio, an improved strange spectral function can be constructed. We shall use it in this last section. 

Besides the strange quark mass, the detailed analysis of the strange r width can bring independent infor- 
mation on the CKM matrix element \V US \- In fact the present uncertainty on \ V US \ is a major contribution to 
the final error on m s . One can turn the argument around and consider the possibility that data on strange 
r decays could ultimately provide a powerful determination of \V US \ (Gamiz et al., 2003, 2005). Figure 27 
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FIG. 26 Determinations of the MS mass m s (4GeV 2 ) based on r spectral functions and on other approaches: K* and 
e + e~ sum rules, lattice calculations. 

illustrates the method: the bands show the regions in the \V US |-vs. -m s (m^) plane which are allowed by the 
different moments <5i?*° for J = 1 + (i.e., longitudinal part excluded) within their total (experimental and 
theoretical) errors. It is found that the various moments are sensitive to |V^ S | and m s in a different way. 
While higher-fc moments exhibit a strong m s dependence, <5i?°° is particularly sensitive to \V US \, almost in- 
dependently of m s . This property can be exploited to achieve a determination of \V US \ from r decays. Since 
the t SU(3) breaking relation (156) involves both \V u d\ and |V^ S |, unitarity of the CKM matrix is used (\V u b\ 
being negligible) and only \V US \ is kept as a variable. 

Using as an independent determination the m s value from recent lattice calculations (Aubin et al, 2004) 
evolved to the r scale, m s (m 2 r ) — (79 ± 8) MeV, one gets from <5i?°° alone 



where the quoted errors are from the t data (dominantly from the strange channels), the theory- at-large 
(Sew, SO4, f v , Jk, a S} m T , the phenomenological longitudinal part), and m s . The theoretical error is 
strongly reduced because the (1 + 0) expansion of the (0,0) moment is well behaved, in contrast to the higher 
moments. Thus the uncertainty on the result is dominated by the experimental error on the branching ratios 
for strange hadronic final states. 

The present result is not yet competitive with the best determinations of \V US \, but it is approaching these. 
Although it could be a statistical fluctuation, the fact that the result (168) is on the low side of the recent 
results from and decays may be worrisome. It could indicate that not all r strange final states have 
been identified. In this context it is important to revisit in particular the branching ratio for the v T K~it + it~ 
mode where the agreement is marginal between ALEPH on the low side, and CLEO and OPAL on the high 
side (Section III.C). 

One could forecast a significant improvement with the analyses of high statistics r samples available at 
the B Factories. The detectors at these facilities (BABAR and Belle) arc equipped with powerful Cherenkov 
counters, well suited to separate kaons and pions in a large momentum domain, and one should expect 
significant improvement in the measurement of the small, statistics-limited strange branching fractions. Then 
it is conceivable that a competitive determination of \V US \ can be obtained. 

We now turn back to the m s determination. In fact, owing to the availability of the different moments, it 
is possible to fit them simultaneously and to derive both \V US \ and m s . In doing so it is mandatory to take 
into account the strong correlations between the moments. With the present best data it seems not possible 
to achieve a consistent description of all five moments (k = — 4) with a unique solution for \V US \ and m s . 



\V US \ = 0.2204 ± 0.0028 cxp ± 0.0003 th ±0.0001 



(168) 
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FIG. 27 The values of \V US \ and m s (m?) derived from the experimental strange and nonstrange moments obtained 
from the ALEPH spectral functions normalized to the world average branching ratios. The various curves correspond 
to the different k, i = moments with k from to 4. The width of the bands indicates the present experimental 
accuracy. The longitudinal components are subtracted using the K and -n branching fractions together with r-ju 
universality, and the phenomenological contributions estimated in (Gamiz et ai, 2003). The data points indicate the 
\V US \ and m s (m^) results from outside r physics. The two |V^ S | values shown correspond to the average of the 
and JiT M 2 results (full circle), and to the y/l — \V u d\ 2 unitarity-based number (open circle) (Isidori, 2005). 

The situation is best illustrated in Fig. 28 showing the 68% CL contours when only pairs of moments are 
considered in the fit. Typical examples are given with the following "adjacent" fc-pairs: going from k = 0, 1 
through k = 3, 4 pairs clearly selects different regions in the plane. The situation is worse than it looks at first 
sight since the moments are highly correlated. Indeed, the correlation coefficient ranges from 0.68 (k = — 4) 
to 0.99 (k — 2 — 3, 3 — 4). A fit of the three k — 0,2,4 moments (adding more does not extract any new 
information from the data) yields a \ 2 of 6.4 for 1 DF, corresponding to a goodness-of-fit probability of only 



We are afraid that the pattern emerging from Fig. 28 undermines the reliability in the previous findings of 
m 8 from r decays, at least within their quoted errors. The more recent determinations (Baikov et ai, 2004; 
Gamiz et at, 2005) have focused on the higher k — 2, 3, 4 moments (only the k = 3 moment is used in (Baikov 
et ai, 2004)) since they are more sensitive to m s and they used a world average value for \V US \. However, 
Fig. 28 shows that the data prefer a lower value, not consistent with the world average. This might be the 
result of the worse behavior of the perturbative expansion of the higher moments as detailed in Eq. (163), even 
for the J — 1 + part. Another sign of potential problems is found in the m s running study as function of Sq 
in Section VIILC.l, where the stability of the results for the higher moments could not be safely established. 

It may be that some of these alarming features will disappear with better experimental data, especially 
if systematic deviations are found. More precise data from the B Factories on the multi-hadronic strange 
final states will certainly add valuable information to this problem. As for the moment, whereas the \V US \ 
measurement from the k — SU(3) breaking moment appears reasonably safe, the m s determination is 
unfortunately plagued with systematic effects and does not appear to be sufficiently robust for a precise and 
model-independent measurement. 

IX. CHIRAL SYMMETRY AND QCD SUM RULES 

In the limit of n/ massless quarks (m., = 0, i = u, d, . . .), the QCD Lagrangian possesses a global SUl(ti/) x 
SU n{nf) chiral symmetry between the left- and right-handed quarks in flavor space. The associated conserved 
Noether currents are the vector (V^) and axial- vector (A^) quark currents. The vector and axial- vector charges 
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FIG. 28 The 68% CL contours in the \V U s |-vs.— m s {m%.) plane from the fit of different pairs of moments derived from 
the experimental strange and nonstrange ALEPH spectral functions normalized to the world average branching ratios. 
The longitudinal components are subtracted using the t-\jl universality-improved K and n branching fractions, and 
the phenomenological contributions estimated in (Gamiz el al., 2003). The results exhibit a systematic trend with the 
orders of the pair chosen. The data points show the best present knowledge of |V« S | (the open circle determined using 
unitarity) and m 3 (m^) from outside r physics. The two \V US \ values shown correspond to the average of the Ktz and 
results (full circle), and to the \/l — \ V u d\ 2 unitarity-based number (open circle) (Isidori, 2005). 




are the generators of the symmetry group. For a state \(f>), which is symmetric under SU(3)l x SU(3)^, one 
then must have (Gasser and Leutwyler, 1982) 29 

MAlWAkvM = (<t>K(x)VZ(y)\<f>) (a,b =,!,..., 8) , (170) 

which requires that for every contribution on the r.h.s. of Eq. (170) (J p = + or 1~) there exists a mass 
degenerate partner on the l.h.s. ( J p = 0~ or 1 + ). However, chiral symmetry, which should be a good symmetry 
for the light u, d, s quarks, is not observed in the low energy hadronic spectra of e + e~ annihilation or r spectral 
functions. In order to be consistent with this experimental fact, the chiral flavor group SU(3)l X SU(3)i{ is 
assumed to be spontaneously broken down to SU(3)l+_r where the vacuum expectation values are symmetric 

(uu) = (dd) = (ss) ^ . (171) 

The spontaneously breaking of the axial charge symmetry generates an octet of massless pseudoscalar 
Goldstone mesons (Goldstone, 1961; Goldstone et al., 1961; Nambu and Jona-Lasinio, 1961), which is identified 
with the 8 lightest hadronic states: 7r + , 7r°, 7r~, rj, K + , K~ , K a and K°. For non-zero quark masses the 
axial-vector current is not conserved and its divergence reads 

d„A^ = (m, + m^q^ij^qj , (172) 

which is associated with the decay constant fp of a pseudoscalar meson P with four momentum g M , defined 

as 

(Old^AplP) = V2f P m 2 P or: {0\A P \P) = V2f P q". (173) 

The non- vanishing physical masses of the pseudoscalars reflect the fact that chiral symmetry is not exact in 
QCD. True zero it, d, s quark masses would generate massless pseudoscalar mesons. Explicit mass relations 



The conserved right- (R = V + A) and left-handed (L = V — A) currents under SU(3)l x SU(3)^j transform like (8, 1) and (1 
8), respectively. If \<f>) is invariant then <^|i?L|0) = =S> ((p\V 2 - A 2 \4>) = 0, which is Eq. (170). 
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between the pscudoscalars can be deduced from Chiral Perturbation Theory (c/. Section IX. D). There is no 
indication of spontaneously breaking of the vector charge symmetry (the vector charge operators annihilate the 
vacuum): no light scalars have been found experimentally, and vector mesons can be classified in degenerate 
multiplets, which are SU(3)y representations. Hence, the longitudinal part of the vector correlator in Eq. (38) 

vanishes in the chiral limit, Il^j v (q 2 ) = 0. 

Reducing chiral symmetry to the nonstrange sector, i.e., to SU(2)l x SU(2)#, one has the Lorentz- 
decomposed two-point correlation function 

KIlr(i) = * J d 4 xe^(0|T(^(x)^(0)t)|0) 

= K^Y + n£Ufo 2 ) + ^Ilr(i 2 ) > (174) 
where IJ 1 and W are the left- and right-handed quark currents 

L M = U7 M (1 — j 5 )d , =B7**(1 + 75 )d . (175) 

With Eqs. (174) and (175) and the correlators for vector and axial-vector currents (38) one finds 

nSU9) = n^(?)-n^(«). (176) 

In the chiral limit and for q 2 — > oo, the correlator 11^ lr(<i) vanishes (which again implies the degeneracy of 
Eq. (170)). With the use of a dispersion relation and from the comparison of the q^q" and the q 2 terms in 
Eq. (174), one obtains the two Weinberg sum rules (WSR) for u, d quark correlators (Weinberg, 1967) 

so^oo 

1. WSR : J dslm [u^ v (s) n£>» + U^ v (s) H^»] = , (177) 







2. WSR: J dssIm[n£l v (s)-n£l A (s) 



(178) 



The upper integration bound so is finite when integrating over experimental spectral functions, and the missing 
saturation of the integral at the cut-off is a source of systematic uncertainty. The first WSR can be simplified 
using the pion decay constant f n defined in Eq. (173), which fixes the integral Jds ImlL^ A (s) = f 2 , and the 
fact the longitudinal vector correlator vanishes. Defining 

Pv-a(s) = vi(s) - ai(s) , (179) 

the two Weinberg sum rules read 

so — >oo 

1. WSR : ^ J d S p V - A (s) = f 2 , (180) 

o 

so — >oo 

2. WSR: J dssp V - A (s) = . (181) 

o 

In the presence of non-zero quark masses, only the first WSR remains valid while the second WSR breaks 
down due to contributions from the difference of non-conserved currents of order m 2 /s, leading to a quadratic 
divergence of the integral (Floratos et al, 1979). However, convergence can be recovered by considering a 
Borel-transformed (also denoted Laplace-transformed) version of the second WSR (see Section IX. C), where 
the result is expressed as a function of the Borel parameter M 2 and the u, d running quark masses at scale M 2 
and quark condensates in powers of 1/M°, 1/M 2 , . . . (Floratos et al, 1979; Peccei and Sola, 1987; Shifman 
et al, 1979). 

Using a coarse narrow-width approximation for the resonances and a saturation of the corresponding vector, 
axial-vector spectral functions by the ir, p(770) and ai(1260), Weinberg deduced from these sum rules the 
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mass formula (Weinberg, 1967) m ai /m p ~ 1.41, which compares not too badly with the experimental value 
of 1.60 ± 0.05 (Particle Data Group, 2004). 

Das, Mathur and Okubo (DMO) (Das et al, 1967b) showed that the derivative of the vector minus axial- 
vector correlator Hud,v-A{Q 2 ) taken at q 2 — is connected with the radiative ir~ — > i~v^ decay axial- vector 
form factor Fa via 30 



d,V-A 



dq 2 



1 

4^2 



ds- p V -a(s) 
s 



r 



-Fa 



(183) 



with the pion charge radius-squared (r 2 ) given in Eq. (66). This relation has been used in the past to perform 
a finite energy sum-rule determination of the electric polarizability 31 of the pion (Hocker, 1997; Kartvelishvili 
et al., 1997; Margvelashvili, 1988; Menke, 1999). Equation (183) is an example for an inverse moment spectral 
sum rule. It is special because the integration kernel is singular at s = so that its residue is sensitive to the 
chiral properties of the spectral functions (cf. Section IX. D). 

The electromagnetic mass splitting between the charged and the neutral pion has been calculated using 
current algebra techniques in chiral symmetry by Das et al. (Das et al., 1967a). They derived the Das- 
Guralnik-Low-Mathur- Young (DGLMY) sum rule in the chiral limit 
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(184) 



where f = f w (l-13ml/(192TT 2 f 2 )-ml(rl)/6 + 0(ml)) = (87.6 ± 0.3) MeV corresponds to = (92.4±0.1± 
0.3) MeV (Particle Data Group, 2004), expressed in the chiral limit (m„ = = 0, m s ^ 0) (Moussallam, 
1999). Making use of the KSFR relation (Kawarabayashi and Suzuki, 1966; Riazzuddin and Fayazuddin, 
1966) and lowest resonance saturation, they obtained the relation 



mlo 



3am 2 p ln2 



5.1 MeV 



(185) 



in approximate agreement with the experimental value of 4.59 MeV (Particle Data Group, 2004). Although 
apparently containing an arbitrary energy scale A, the sum rule (184) is actually independent of A by virtue 
of the second WSR. However, for finite s a residual A dependence is maintained. 



A. Chiral sum rules with r spectral functions 

The spectral integrals (177) through (184) are computed with variable upper integration bounds so < m r 
using the experimental r hadronic spectral functions and their respective covariance matrices. The availability 
of the covariance matrices allows one to perform a straightforward error propagation taking into account the 



30 In cases where the radiated photon is real the differential decay rate can be written as 

d 2 r^ e ^ f = d 2 (r IB +r gD + r INT ) 

dE 1 dE l dE 1 dE l 



(182) 



where Tib, Tgrj and Tjnx are the contributions from inner bremsstrahlung, structure-dependent radiation and their interfer- 
ence. The structure-dependent term is parameterized by using vector and axial- vector form factors Fy/A (f° r t ne complete 
expression see Ref. (Particle Data Group, 2004)). 

The electric polarizability of a physical system can be classically understood as the proportionality constant that governs 
the induction of the dipole moment p of a system in presence of an external electric field vector E: p = agE. The polarizability 
is an important quantity to characterize a particle, i.e., in probing its inner structure. In chiral perturbation theory, it can be 
shown (Holstcin, 1990; Terent'ev, 1973) that the polarizability of the pion is given by = ciFa / {m n f 2 ) , which gives (Hocker, 
1997) a E = (2.86 ±0.33) X 10~ 4 fm 3 . Using the value (66) for the pion charge radius-squared, the analysis of a Borcl-improved 
DMO sum rule yields (Hocker, 1997) a E = (2.96 ± 0.32) X 10~ 4 fm 3 . 
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FIG. 29 The first and second Weinberg sum rules (upper left and right), and the DMO and DGLMY sum rules (lower 
left and right) versus the upper integration cut-off (so). The horizontal bands indicate the chiral predictions within 
errors, assuming full saturation of the integrals. 



strong bin-to-bin correlations of the spectral functions. The vector minus axial- vector r spectral functions 
measured by ALEPH (ALEPH Coll., 2005) and OPAL (OPAL Coll., 1999b) have been shown in the right hand 
plot in Fig. 5. Local saturation of the spectral functions at the integral end point s would require convergence 
of central values and all derivatives at zero. This is not yet achieved at sq = m%- We hence expect that the 
sum rules are not entirely satisfied for the r spectral functions. The sum rules (177) through (184) versus the 
upper integration cut-off So < rr? T are plotted in Fig. 29 for the ALEPH spectral functions. The horizontal 
band depicts the chiral predictions of the integrals, neglecting OPE power corrections. One observes that only 
for the DMO sum rule (lower left plot), where contributions from higher masses are suppressed, saturation 
at the t mass is approximately satisfied within errors. The other sum rules are not saturated at m^, and no 
further conclusions can be drawn. 

Many V — A sum-rule analyses have been performed since the experimental spectral functions became 
available (Bijnens, 2001; Cirigliano et at, 2003a, b; Chilli et at, 2003; Davier et at, 1998; Dominguez and 
Schilcher, 2004; Ioffe and Zyablyuk, 2001; Peris et at, 2001; Rojo and Latorre, 2004; Zyablyuk, 2004). Some 
results of these will be discussed in the following. 
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B. Operator product expansion 

Nonperturbative effects give rise to nonperturbative OPE contributions to the V — A correlation function 

ngyW-ng A ( s )= J2 r4m E c D , v _ A (s,^{o D ^)) v _ A , (i86) 

£>=2,4,6,... ^ ' ' dimO=D 

where we shall omit the V — A subscripts in the following. The dimension D = contribution is of pure 
perturbative origin and is degenerate in all orders of perturbation theory for vector and axial- vector currents. 
To NLO in QCD one has 



Cd(s, h){Gd{h)) = a D (fi) + feoln ( — 



(187) 



The explicit logarithmic dependence of the Wilson coefficients (the bo terms) is known for D = 2,4,6, but 
not for higher dimensions. Useful relations between the OPE terms and the spectral functions are obtained 
from projective sum rules (Kartvelishvili and Margvelashvili, 1992). When neglecting the logarithmic power 
terms in Eq. (187), one finds (Cirigliano et al., 2003a) 



hl ds u) 



pv - a{s) = ~h 
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a D 



|s|=s 



(--■) i;/2 Co 



(-1) -F d ^(I3/2)-l , 



(188) 



where S m ^ n has the value 1 if m = n, and elsewhere. 

The validity of the OPE depends on the size of the cut-off scale s . Apart from the saturation of the 
spectral integral, this represents another source of systematic uncertainty related to the cut-off. The insertion 
of appropriate weight functions w(s) under the spectral integral in (188) may help to reduce the impact of these 
systematics (see also the discussion in (Cirigliano et al, 2003a, b)). While the impact of missing saturation 
can be reduced by suppressing the high-energy tail of the V — A spectral function, the breakdown of the OPE 
due to quark-hadron duality violations is unavoidable to some extent. In practice, this duality violation is 
assumed to be absorbed by the OPE power operators that are fit phenomenologically to the spectral sum 
rules. The fit protects the physical observable (for instance the strong coupling constant as discussed in 
Section VII) that one is aiming to extract. If the nonperturbative corrections are found to be large, this 
protection mechanism may not work sufficiently well anymore, and duality violations can entail systematic 
errors. The optimization between saturation of the integral on one hand, and small duality violation on the 
other hand requires a compromise: if one aims at a minimum saturation loss, i.e., small residual modulations 
of the spectral function in the vicinity of the cut-off s , one has to apply a weight function that suppresses 
the high-energy tail 32 . However, the stronger the suppression the larger the quark-hadron duality violation 
will be. 

Using the formulae of (Braaten et al., 1992) and (Chetyrkin and Kwiatkowski, 1993) for the nonperturbative 
power expansion of the correlators, one obtains for the V — A correlators the OPE (Davier et al., 1998) 



n (0+1) (s) 



a s (s) m 2 (s) 



-87r 2 a s (M 2 ) 
2ir 2 



3 + 4L(s) 



8 . . 59 9 . ,\m(uu + dd) 16 m 4 (s) 



1 



1 ' ^ ■ 



(0 2 ( M 2 )) (Qg) 



(189) 



The suppression of the high-energy spectral function also reduces the experimental uncertainty on the integral since the 
statistical errors are large at the kinematic endpoint of the t phase space. 
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2a s " 1 ( S )-5 



21373 75 



m(uu + dd) 



7tt 2 



2448 17 



C(3) a s (s) 



TO 4 (s) 



(190) 



with a s (s) = a s (s)/ir, L(s) = ln(s/y(i 2 ) and the dimension D = 6 operators 

0\ = u^ 5 T a dd^ 5 T a u - u lf ,T a dd^T a u , 
0\ = ujfjddj^u - u-f^dd-f^-f^ , 



(191) 



and where the SU(3) generators T a are normalized so that tr(T a T 6 ) = S ab /2. The average mass m stands for 
(m u +m<i)/2, i.e., SU(2) symmetry is assumed. The constant C{y?) depends on the renormalization procedure 
and should not affect physical observables. Dimension D = 2 mass terms are calculated perturbatively to 
order a|, which suffices for the light u, d, quarks. The coefficient functions of the dimension D = 4 operators 
for vector and axial-vector currents have been calculated to subleading order in (Chetyrkin et al, 1985; 
Gencralis, 1989). Their vacuum expectation values are expressed in terms of the scale invariant gluon and 
quark condensates. Since the Wilson coefficients of the gluon condensate are symmetric for vector and axial- 
vector currents, they vanish in the difference. The expectation values of the dimension D = 6 operators (191) 
obey the inequalities {0\(p?)) > and (0 2 (^i 2 )) < 0, which can be derived from first principles. The 
corresponding coefficient functions were calculated in the chiral limit (Lanin et al., 1986). Since (C 2 (^ 2 )) 
is subleading in a s (^ 2 ), it is neglected in most phenomenological analyses. In leading order, one can use 
factorization of the matrix elements (that is, vacuum insertion) to rewrite {0\) in the form 33 



64 - 

{0\) oc -— 7r 2 a s (0|uii|0)(0|dd|0) 



(192) 



For the dimension D — 8 operators no higher order calculations are available in the literature, so that it 
is customary in the QCD sum-rule analyses to neglect their logarithmic s dependence. Again, the J = 
contribution vanishes in the chiral limit. A factorization result similar to Eq. (192) can be derived, which now 
involves a five-dimensional quark-gluon mixed condensate (Grozin and Pinclis, 1986; Ioffc and Zyablyuk, 2001; 
Zyablyuk, 2004). Higher OPE dimensions are neglected (see (Zyablyuk, 2004) for an analysis that includes 
the relevant D = 10 operators). 



C. Borel and pinched weight sum rules 

Many ways have been explored to improve the saturation of the integrals at finite cut-off scale sq, and to 
better control the validity of the OPE. In general, the approaches distinguish themselves by the definition of 
the weight w(s) used under the integral 

so 

-^-2 J dsw{s)p V - A {s) . (193) 



Since duality violations are expected to be largest in the vicinity of the timclike real axis (Poggio et al., 
1976), some authors suggest to use weights of the form w n (s) oc (1 — s/s )™ (Le Diberder and Pich, 1992b) (cf. 
Section VII. C). The class of functions that satisfy w(sq) = is introduced as "pinched weights" in (Cirigliano 
et al, 2003a, b). As mentioned before, pinched weights improve the experimental accuracy of the sum rule, but 
also reduce the effective inclusiveness of the spectrum, and hence make it more sensitive to duality violations. 

Another way to improve the saturation of finite energy sum rules is to apply the Borel weight w(s) = 
cxp(— s/M 2 ), emphasizing the lower resonances of the spectrum. Here, the OPE operators of dimension 



Following (Zyablyuk, 2004), factorization for the D = 6 operators can be written as ((?7Ari(i)(a!Ar2w)} = — (16/9)Tr[(tt (g) 
u)Tl (d ® d)T2], where I\ are Dirac matrices (acting on color, flavor and Lorentz indices), the factor 16/8 includes the SU(3) 
color number, and the 4x4 matrix (q <g> q) is proportional to the quark condensate (q (g> q) = — (l/4)(0|gg|0). 
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FIG. 30 The Borel sum rule (194) (left) and the real part of the complex Borel sum rule (196) for (j> — 97r/10 (right), 
at M = 1 GeV, and as a function of the integral cut-off so. The chiral prediction shown in the right hand plot does 
not include OPE power corrections. 



D = n are suppressed as (02n+2)/«!, and all leading terms of the OPE vanish. The Borel-transformed 
sum rules give rise to nonperturbative corrections from quark condensates and non-vanishing quark masses 
that scale with the Borel parameter M 2 and with a s (M 2 ). For example, the Borel-transformed first WSR 
reads (Floratos et at, 1979; Narison, 1982; Pascual and de Rafael, 1982; Peccei and Sola, 1987) 

so — >oo 

f dse-^ 2 p v ^ A (s)=f 2 





a s (M 2 ) 
2tt 



m u (M 2 )m d (M 2 ) 8 f 2 m\ 32tt ((0|gg|0)) 2 



(194) 



7T 2 3 M 2 9 M 4 

where the lowest-order current algebra Gell-Mann-Oakes-Renner relation 

f 2 ml = -{m u + m d )(0|uu|0) , (195) 

obtained from Eqs. (172) and (173) has been used (Gell-Mann et al, 1978), and where m u {M 2 ) and m d (M 2 ) 
denote the running light quark masses (see (Jamin, 2002) for a derivation of chiral corrections to Eq. (195), 
which are found to be of the order of —5%). The Borel sum rule is plotted in Fig. 30 (left hand) for 
M = 1 GeV. The comparison with Fig. 29 (upper left hand plot) exhibits the gain in experimental precision 
and in the approach to the asymptotic regime. Inserting the experimental values for the parameters on the 
r.h.s. of Eq. (194), and neglecting higher order nonperturbative contributions (D(M~ 6 ) and contributions 
other than quark condensates, one obtains (8.44 ± 0.17) x 10 -3 GeV 2 . From Fig. 30 we conclude that further 
nonperturbative power terms need to be added to satisfy the sum rule. 

To reduce the correlation between the different operators, one can consider a Borel transformation in the 
complex plane of the Borel scale parameter, that is M 2 — ► M 2 exp[i(7r — <f>)]. Keeping only the leading order 
coefficients ao in Eq. (187), one finds for the real and imaginary parts of the Borel transformation (Zyablyuk, 
2004) 



SO — ►OQ 



47T 



— J e^smis^pv-Ais) - Z> X ) fc!M 2 fc+2 ■ ( 197 ) 



with Stj, = (s/M 2 ) sin <fi, c$ = (s/M 2 ) cos 4>. If the sum rule is applied to the r spectral functions, the argument 
of the cosine function must be negative to suppress the unknown sq > m 2 tail of the spectrum. In the case of 
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FIG. 31 Real part (196) of the Borel transformation for <j> = 97r/10 (left), and its imaginary part (196) for <j> = 4n/5 
(right), as a function of the modulus M 2 of the Borel parameter. The leading order D — 12 operator an vanishes 
in these sum rules. The shaded band shows the ALEPH data (ALEPH Coll., 1997b, 1998c), and the lines display 
the leading order operator series, i.e., neglecting all logarithmic terms. The numbers indicate the highest dimension 
considered in the operator sum. The figures are taken from (Zyablyuk, 2004). 



4> = 7r(2fc - 1) /(2k) (fc = 2,3, . . .) for the real part (196), and = n(k - l)/k (k = 3, 4, . . .) for the imaginary 
part (197), the sum rules become projective with vanishing leading operators of dimension D = 2k + 2. The 
real part for a particular working point chosen in (Zyablyuk, 2004) as a function of the cut-off sq is shown 
in the right hand plot of Fig. 30. The chiral prediction given contains the pion pole, but no OPE power 
corrections. The sum rule is well saturated at s = m r- The remaining discrepancy with the data is due to 
nonperturbative contributions. An update (Zyablyuk, 2004) of the numerical analysis performed in (Ioffe and 
Zyablyuk, 2001) finds for the leading V - A operators using the ALEPH data (ALEPH Coll., 1998c) 

a 6 = -(7.2 ±1.2) x 10~ 3 GeV 6 , (198) 
a 8 = (7.8 ±2.5) x 10~ 3 GeV 8 . (199) 

The author also estimates the leading D — 10 operator to be a w = —(4.4 ± 2.8) x 10~ 3 GeV 10 , without 
however reporting the statistical correlations with ciq and as. The imaginary part (197) of the Borel sum rule 
were used as input, as well as the angles <fr — 2tt/3, 37r/4, 47r/5 for which the operators as,aio,&i2 vanish, 
respectively. For the Borel parameter M 2 a range between 0.4 and 1.0 GeV 2 is used, which is integrated over 
to determine an average \ 2 estimator (Zyablyuk, 2004) that is minimized in the fit. The fit results for the 
real and imaginary parts of the sum rules (196), (197), and for two working points in the parameter <j> are 
shown in Fig. 31. The contribution from D = 10 is significant here. 

Recent theoretical developments match the D = 6 operator, derived from the V — A spectral sum rules, 
with the vacuum matrix elements and (Os)^, which are related to the electroweak penguin operator 

Qs that contributes to K — > 7T7t decays (see, e.g., Refs. (Bijncns, 2001; Cirigliano et at, 2003b)). Together 
with the QCD penguin Qq, these operators generate direct CP violation in neutral kaon decays to two pions, 
the computation of which suffers from large uncertainties. 



D. Inverse moments sum rules and chiral perturbation theory 

We have seen that at sufficiently high energies nonperturbative effects are parameterized by vacuum con- 
densates, following the OPE rules. A particular role is played by the condensates that are order parameters 
of the spontaneous breakdown of chiral symmetry (SB^S). They vanish at all orders of perturbation theory 
and control the high energy behavior of chiral correlation functions, such as the difference of vector and axial- 
vector current two-point functions. On the other hand, at low energies, SB%S makes it possible to construct 
an effective theory of QCD, the Chiral Perturbation Theory (x?T) (Gasser and Leutwyler, 1984; Weinberg, 
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1979), which uses the Goldstone bosons as fundamental fields and provides a systematic expansion of QCD 
correlation functions in powers of momenta and quark masses. Any missing information is then parameterized 
by low-energy coupling constants, which can be determined phenomcnologically in low-energy experiments 
involving pions and kaons. The fundamental parameters describing chiral symmetry breaking, the running 
quark masses and the quark condensates appear both in low-energy x?T and the high energy OPE expan- 
sion. For this reason it is useful to combine the two expansions in order to get a truly systematic approach 
to the chiral sum rules (Donoghue and Golowich, 1994). Such a combined approach is illustrated in (Davier 
et al., 1998) through the determination of the Lio parameter of the chiral Lagrangian, including high-energy 
corrections coming from the OPE. The connection between the two domains is provided by the experimental 
t spectral function pv-A- 

At leading order %PT, L w is obtained by the DMO sum rule (183) 

so — >oo 

J ds^ PV -a(s) -4Li • (200) 



As it stands the DMO sum rule (183) is subject to chiral corrections due to non-vanishing quark 
masses (Golowich and Kambor, 1997). On the other hand, the truncation of the integral at so < m 2 in- 
troduces an error that competes with the low-energy chiral corrections. Both types of corrections can be 
systematically included through (i) the high-energy expansion in a s (so) and in inverse powers of s , and 
(ii) the low-energy expansion in powers of quark masses and of their logarithms. 

In Section VII. C we have reviewed the simultaneous determination of the strong coupling constant and the 
OPE operators using the spectral weights 

•" ,w,w =K)"fe)'- (201 » 

for k,£ > 0. The extension of the spectral moment analysis to negative integer values of I ("inverse moment 
sum rules", (IMSR) (Golowich and Kambor, 1996)) requires, due to the pole at s = 0, a modified contour 
of integration in the complex plane, as illustrated in Fig. 32. This is where xPT comes into play: along the 
small circle Ci placed at the production threshold s t h = 4mJ one can use xPT predictions for the two-point 
correlators. Generalizing the spectral moments (129) to IMSR's, one obtains for the Cauchy integral 



Rl k f /A = 6ni\V ud \ 2 S EW f^w^(s) 



(202) 



where C = C\ + C2 for the inverse moments and C = C\ for the positive moments (see Fig. 32). 
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The nonstrange correlators (38) are known up to two loops (Gasser and Leutwyler, 1984, 1985; Golowich 
and Kambor, 1995, 1998) in standard xPT 34 

The analysis described in (Davier et al., 1998) performs a fit of the OPE to the IMSR R^~y* A and the 

spectral moments (129). The inner integration circle C 2 occurring in R^y* A introduces the sensitivity to the 

effective xPT parameter Z^o > which is determined by the fit. This quantity is a well defined observable that 
can be related to L"[ (ii x p T ) at each loop in xPT. 

The spectral information is used to simultaneously determine L°q and the nonperturbative phenomeno- 
logical operators. For dimension D = 6 the contribution from 0\ is neglected since it is large-Ac and a 2 s 
suppressed. The results of the fit are reported as (Davier et al., 1998) 

Lf = -(6.36 ± 0.09 cxp ± 0.14th ± 0.07 nt ± 0.06 OPE ) x 10~ 3 , (203) 
(O e ) = (5.0 ± 0.5 cxp ± 0.4th ± 0.2 fit ± I.Iope) x 10~ 4 GeV 6 , (204) 
(O s ) = (8-7 ± 1.0 exp ± 0.1 t h ± 0.6fi t ± 2.1 OPE ) x lO" 3 GeV 8 , (205) 

where Oq = 0\(m 2 ). The fit errors are separated in experimental (first errors) and theoretical (second errors) 
parts, and systematic fit uncertainties (third errors) are added. The last errors estimate the uncertainties due 
to the OPE breakdown (ALEPH Coll., 1998c). They are small for L°q and dominant for the nonperturbative 
operators. Correcting for the factor of about —13 between the definition of (Oq) and a§ of Eq. (198), the two 
estimates are in agreement. The results (204) and (205) have been confirmed in (OPAL Coll., 1999b) using 
OPAL data. In general it is observed that while all available sum-rule analyses 35 agree on the value for the 
D = 6 contribution to the V — A correlator, the D = 8 contribution, due to its smallness, is less stable. It 
exhibits a model dependence that exceeds the experimental error in (205). 

A stability test of the results (203) through (205) at m 2 has been performed in (Davier et al, 1998) by fitting 
with variable "r masses" s < m 2 (Davier et at, 1998). The resulting convergence was found satisfactory. 

Expressing of Eq. (203) at the renormalization scale /j x pt = m-p — 770 MeV, one obtains (Davier et 
al, 1998) L r w (m p ) = —(5.13 ± 0.19) x 10~ 3 , where experimental and theoretical errors of different sources 
have been added in quadrature. This result can be compared with the value of L\ that is obtained from the 
one-loop expression of the axial form factor, F A , of ir~ — ► e~T7 e ~( decays (see (Bijnens and Talavera, 1997) for 
notations) 

F A = ^(L 9 +I 10 ) . (206) 
fa 

Using L r g (m p ) = (6.78 ±0.15) x 10~ 3 (Bijnens et al, 1995a; Davier et al, 1998; Gasser and Leutwyler, 1985) 
and F A = 0.0116 ± 0.0016 (Particle Data Group, 2004), one obtains L 9 + L 10 = (1.36 ± 0.19) x 10~ 3 . This 
gives L r w (m p ) — —(5.42 ± 0.24) x 10~ 3 , in agreement with the r result. 



X. CONCLUSIONS AND PERSPECTIVES 

The physics of hadronic r decays has been the subject of much progress in the last decade, both at the 
experimental and the theoretical level. Somewhat unexpectedly, hadronic r decays provide one of the most 
powerful testing grounds for QCD. This situation results from a number of favorable conditions: 

• The t lepton is heavy enough to decay into a variety of hadrons, with net strangeness and ±1. 

• t leptons are copiously produced in pairs at e + e~ colliders, leading to simple event topologies without 
extra background particles. 



The authors of (Davier et al, 1998) use generalized xPT (G*PT (Knecht and Stern, 1995)) to 0(p 4 ) one-loop order only. 
Standard xPT assumes (Weinberg, 1977) the validity of Eq. (195) and r ~ 2m 2 K /m'^ — 1 ss 25.9, whereas GxPT admits lower 
values of these two quantities (Fuchs et al., 1990, 1991; Knecht and Stern, 1995). The alterations of the standard 0(p 4 ) results 
for nonstrange correlators introduced by GxPT are marginal. They merely concern the symmetry breaking J = component 
of the spectral functions and most of them are absorbed into the renormalization of f^. 

See also the detailed recent analysis (Cata et al, 2005) on duality violations in V — A sum rules, which came to late to be 
included in this review. 
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• The experimental study of r decays could be done with large data samples: at LEP with high efficiency 
and small background; at the Y(4S) energy with much increased statistics, but irreducible backgrounds. 
The two experimental setups offer complementary advantages. 

• As a consequence, r decays are experimentally known to great details and their rates measured with 
high precision. 

• The theoretical description of hadronic r decays is based on solid ground and found to be dominated by 
perturbative QCD, as the result of several lucky circumstances, which shall be briefly recalled below. 

The t decay rates into hadrons are expressed through spectral functions of different final states with 
specific quantum numbers. The spectral functions are the basic ingredients to the theoretical description of 
these decays, since they represent the probability to produce a given hadronic system from the vacuum, as a 
function of its invariant mass-squared s. The spectral functions are dominated by known resonances at low 
mass, and tend to approach the quark-level asymptotic regime towards m 2 .. They have been determined for 
individual hadronic modes, however they are conveniently summed and separated into their inclusive vector 
and axial- vector components for the nonstrange part, while the Cabibbo-suppressed strange part cannot be 
fully separated at present, due to the lack of necessary experimental information. 

The nonstrange r vector spectral functions can be compared to the corresponding quantities obtained in 
e + e~ annihilation by the virtue of isospin symmetry. The 0(%) precision reached makes it necessary to correct 
for isospin-symmetry breaking. Beyond a superficial agreement, the detailed comparison unveils discrepancies 
with some e + e~ data sets. Since the vector spectral functions are the necessary ingredients to compute vacuum 
polarization integrals, required for the evaluation of the running of a(s) or the anomalous magnetic moment 
of the muon a M , this disagreement leads to different results when using r or e+e~ spectral functions. While 
the e + e~-based theoretical a M value disagrees with the measurement by 2.7 "standard deviations", possibly 
indicating a contribution from physics beyond the Standard Model, the r-based calculation is consistent with 
experiment. Although the use of e + e~ data is a priori more direct, the physics at stake is such that the 
present situation must evolve, requiring more experimental and theoretical cross-checks. Newest e + e~ results 
already indicate better agreement with the r data. 

The observation that hadronic r spectral functions can provide precision information on perturbative QCD 
is surprising, considering the moderate energy scale involved in these decays. Hence it is useful to recall briefly 
the reasons behind this success: 

• The total semileptonic decay rate normalized to the leptonic width, i? T , is obtained by integrating over 
s the total spectral function weighed by a known kinematic function originating from the V — A leptonic 
tensor. This integral can be transformed to a complex contour integral over a circle at |s| = m 2 ., hence 
involving only large complex s values where perturbative QCD can be applied. 

• One factor in the weight function, (1 — s/to 2 ) 2 , chops off the integrand in the vicinity of the real axis 
where poles of the current correlator are located. 

• The invariant-mass spectra in the hadronic r decays, dominated by resonances, are thus related to the 
quark-level contributions expected at large mass. This global quark-hadron duality is expected to work 
best in more inclusive situations. In the case of hadronic r decays, this condition is particularly well met, 
with several resonance contributions in each of equal-importance vector and axial- vector components. 

• The perturbative expansion of the spectral function is known to third order in a s , with some estimates 
of the fourth order coefficient. 

• The perturbative and nonperturbative contributions can be treated systematically using the operator 
product expansion (OPE), giving corrections organized in inverse powers D of m T . 

• The a priori dominant nonperturbative term from the gluon condensate (D = 4) is strongly suppressed 
by the second factor of the weight function, (l + 2s/m 2 ). This accidental fact renders hadronic t decays 
particularly favorable since nonperturbative effects are expected to be small. 



The next D — 6 term is subject to a partial cancellation between the vector V and axial- vector A 
contributions, due to the (1 — 75) relative factor. 
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These fortunate circumstances explain why hadronic r decays can be globally described by perturbative 
QCD with very small nonperturbative corrections. This behavior has been verified model-independently using 
the measured spectral functions. From a simultaneous fit of the total hadronic rate, obtained with precision 
from the leptonic branching ratio, and spectral moments that are weighted integrals of the mass-squared 
distributions, the total nonperturbative contribution is found to be 5np,v+a = (— 4.3 ± 1.9) x 10~ 3 relative 
to the V + A hadronic width. This is almost a factor of 50 smaller than the a s -dependent perturbative 
correction to R T . The predicted near-cancellation between the V and A contributions has been also observed, 
as 5 NP y = (2.0 ± 0.3) x 10~ 2 and S NP , A = (-2.8 ± 0.3) x 10~ 2 . 

Several approaches are available to treat the perturbative expansion of the r hadronic width. In fixed-order 
perturbation theory (FOPT), all contributions beyond the known orders of the Adler function are neglected, 
even if some parts resulting from the running of a s are calculable (and calculated). On the contrary, the 
contour-improved FOPT variant (CIPT) keeps all known terms. It has been shown that a faster convergence 
is achieved with CIPT and, since the FOPT expansion represents an unnecessary approximation, the CIPT 
approach should be preferred. Analyses using ALEPH and OPAL spectral functions give a s (m 2 ) values that 
are in agreement. The robustness of the result can be established by comparing results from separate fits based 
on the V or A spectral functions, which are less inclusive and include larger, yet still small, nonperturbative 
contributions: their results are in reasonable agreement. A further test is achieved by partial integration of 
the spectral function up to a variable hypothetical r mass so less than m 2 : the derived a s (so) values are 
in agreement with the expected behavior from the QCD renormalization group down to sq values as low as 
1 GeV 2 , where the OPE begins to break down. Therefore the moderately high mass of the r lepton is shown 
to be fully adequate for a precision measurement of a s . 

Using the final spectral functions from ALEPH, supplemented by the world average leptonic and total 
strange hadronic fractions, the most precise determination of the strong coupling at the r mass scale gives 

a s (m 2 ) = 0.345 ± 0.004 cxp ± 0.009 th ■ (207) 

Evolving this value to the Z mass, and adding experimental and theoretical errors in quadrature, gives 
a s (M§ ) T = 0.1215±0.0012. Taken together with the measurement from the Z width at LEP, a s {M\)z width = 
0.1186±0.0027, the a s (m 2 ) determination provides a powerful test of the running of a s as predicted by QCD. 
The property of asymptotic freedom can be conveniently tested using the estimator 

Aa- 1 

r(si, s 2 ) = A1 s r- = 1.405 ± 0.053 , (208) 
Amy's 

with si = to 2 and s 2 = M§, while QCD predicts 1.353 ± 0.006. 

Further progress in the QCD analyses of r decays will have to wait for a full calculation of the fourth- 
order coefficient K4 in the Adler function expansion. Such a program is underway and will succeed with the 
availability of improved computing power. When K 4 is known, we expect the perturbative uncertainty to 
become really small: the errors on a s (m 2 ) from K 5l estimated assuming a geometric growth of the perturbative 
coefficients with 100% uncertainty, and from the renormalization scale will be 0.002 and 0.002, respectively. 
Using the existing data, the total error on a s (m 2 ) will be reduced by a factor of 2 and reach a relative 
precision of 1.4%. Hence the calculation of K4 will be highly rewarding. 

Further progress can be considered with better quality data. Whereas ameliorating the precision on the 
branching ratios appears to be difficult, several areas could be improved with facilities such as the B Factories, 
affording very large statistics samples. Examples include the study of the high mass part of the spectrum near 
the r-mass endpoint, limited at the moment by statistics. This will improve the knowledge of the spectral 
function in an interesting region and permit the study of the approach to the asymptotic regime. Such topics 
can also be covered at a r-charm factory operating near the r threshold. Detailed physics simulations are 
required to quantitatively assess the potential of these measurements. 

The strange component of the r hadronic width gives access to the study of flavor-symmetry breaking 
induced by the massive s quark as compared to the nearly massless u and d quarks. Although the expected shift 
in the strange hadronic width is indeed found experimentally, it turns out that a competitive determination 
of the strange quark mass m s is difficult because of poor convergence of the mass-dependent perturbative 
series and lack of sufficient inclusiveness in the higher-order spectral moments that are most sensitive to m s . 
Thus various analyses have given somewhat different results, depending on the moments used. The situation 
might improve with a better determined strange spectral function. Contrary to the situation for m s , the 
determination of the CKM matrix element \V US \ from the strange r decays appears on safe grounds, since it 
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relies directly on the non-weighted spectral function with maximum inclusiveness and a better control of the 
perturbative expansion. Present data yield 

| K s | = 0.2204 ±0.0028 exp ± 0.0003th ±0.0001 TOs , (209) 

where m s has been taken from lattice calculations to be m s (m^) = (79 ± 8) MeV. 

The study of strange r decays can benefit significantly from the large statistical samples provided by the B 
Factories, complemented by the excellent K — ir separation available with the BABAR and Belle detectors. The 
prospects are excellent for improving over the present situation. As discussed above, the determination of \ V US \ 
is limited by the experimental error, while the theoretical uncertainty would permit the best measurement 
among all of this quantity. We express our wish that these experiments measure the various strange decay 
channels in a systematic way towards achieving this goal. It could also occur that the present unsatisfactory 
situation on the determination of m a will improve during this process. In any case, better strange spectral 
functions are needed. 

Another active field of research is the study of finite energy chiral sum rules, involving the nonstrange V — A 
spectral function, which are of intrinsically nonperturbative nature. The analyses are divided into spectral 
moments that are weighted integrals over the V — A spectral function with some high-mass damping function 
to reduce the systematic uncertainty due to the finite energy integral cut-off, and inverse moment sum rules, 
which, due to the singularity at s = 0, give direct access to the chiral structure of the QCD Lagrangian. 
The former sum rules are useful to analyze the nonperturbative structure of the OPE , since it allows one to 
measure the higher dimensional operators of the expansion independently of the perturbative contributions. 
Several analyses give agreeing results on the D — 6 contribution, while the much smaller D = 8 contribution 
exhibits some model dependence. Even D = 10 has been considered and found to improve the fit, though 
not significantly. No sign of nonconvergence of the OPE has been detected at the r-mass scale. Inverse 
moment sum rules have been traditionally used through the Das-Mathur-Okubo sum rule to determine the 
polarizability of the pion. In a more evolved later analysis the parameter L w of the chiral Lagrangian has 
been determined, which at one-loop and expressed at the Chiral Perturbation Theory (xPT) renormalization 
scale h x pt = m p ~ 770 MeV, is found to be 

L w( m p) = -(5.13 ±0.19) x 10~ 3 , (210) 

where experimental and theoretical errors of different sources have been added in quadrature. This result 
is in agreement with, but more precise than corresponding determinations from the axial form factor of 
7r~ — » e~~F e j decays. 

While the study of hadronic r decays has achieved a high status thanks to precise and complete sets of 
measurements and the development of the relevant theoretical framework, it is satisfying that advances are 
at hand both at the experimental and theoretical levels. 
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